Bond breaking in epoxy systems: A combined QM/MM approach
- Air Force Research Laboratory, Materials and Manufacturing Directorate, WPAFB, Ohio 45433 (United States)
- Engility Corporation, Air Force Research Laboratory, WPAFB, Ohio 45433 (United States)
- Miami University, Oxford, Ohio 45056 (United States)
A novel method to combine quantum mechanics (QM) and molecular mechanics has been developed to accurately and efficiently account for covalent bond breaking in polymer systems under high strain without the use of predetermined break locations. Use of this method will provide a better fundamental understanding of the mechano-chemical origins of fracture in thermosets. Since classical force fields cannot accurately account for bond breaking, and QM is too demanding to simulate large systems, a hybrid approach is required. In the method presented here, strain is applied to the system using a classical force field, and all bond lengths are monitored. When a bond is stretched past a threshold value, a zone surrounding the bond is used in a QM energy minimization to determine which, if any, bonds break. The QM results are then used to reconstitute the system to continue the classical simulation at progressively larger strain until another QM calculation is triggered. In this way, a QM calculation is only computed when and where needed, allowing for efficient simulations. A robust QM method for energy minimization has been determined, as well as appropriate values for the QM zone size and the threshold bond length. Compute times do not differ dramatically from classical molecular mechanical simulations.
- OSTI ID:
- 22676023
- Journal Information:
- Journal of Chemical Physics, Vol. 144, Issue 24; Other Information: (c) 2016 U.S. Government; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Development and application of QM/MM methods to study the solvation effects and surfaces
Photodissociation dynamics of CH{sub 3}C(O)SH in argon matrix: A QM/MM nonadiabatic dynamics simulation