CN112199909A - Method for accurately calculating absolute free energy of gas molecules - Google Patents

Method for accurately calculating absolute free energy of gas molecules Download PDF

Info

Publication number
CN112199909A
CN112199909A CN202011141320.0A CN202011141320A CN112199909A CN 112199909 A CN112199909 A CN 112199909A CN 202011141320 A CN202011141320 A CN 202011141320A CN 112199909 A CN112199909 A CN 112199909A
Authority
CN
China
Prior art keywords
energy
free energy
gas molecules
potential energy
mol
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011141320.0A
Other languages
Chinese (zh)
Inventor
彭春望
马健
温书豪
杨明俊
孙广旭
赖力鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Jingtai Technology Co Ltd
Original Assignee
Shenzhen Jingtai Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Jingtai Technology Co Ltd filed Critical Shenzhen Jingtai Technology Co Ltd
Priority to CN202011141320.0A priority Critical patent/CN112199909A/en
Publication of CN112199909A publication Critical patent/CN112199909A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method for accurately calculating absolute free energy of gas molecules, which comprises the following steps: s1, reducing the energy barrier of the flexible dihedral angle potential energy surface in the gas molecules; and S2, calculating absolute free energy of gas molecules. The invention realizes the automatic identification and potential energy surface scanning of flexible dihedral angles in molecules, realizes the fitting of the dihedral angle potential energy surface, automatically modifies the dihedral angle parameters to reduce the potential energy surface to a specific value, realizes a new WHAM program, can obtain actual free energy from potential energy reweights of biased and unbiased, designs a set of feasible calculation flow, and realizes the gradual change of potential energy when gas molecules gradually remove various types of energy; all the processes are connected together, and the automatic calculation of the absolute free energy of the gas molecules is realized.

Description

Method for accurately calculating absolute free energy of gas molecules
Technical Field
The invention belongs to the technical field of molecular dynamics simulation, and particularly relates to a method for accurately calculating absolute free energy of gas molecules.
Background
In recent years, molecular simulation techniques have been widely used, and there has been an increasing demand for free energy calculations. Currently, because the stability of a crystal structure is better and the sampling space is relatively smaller, some methods can accurately predict the absolute free energy of a solid molecule and can be applied to prediction and screening of a stable experimental crystal form of a drug molecule.
There is no good accurate prediction method for the absolute free energy of liquid and gas molecules. This makes it impossible to predict the properties related to solid-liquid (e.g. melting point, solubility), solid-gas (e.g. enthalpy of sublimation), etc. very effectively.
The calculation of the absolute free energy of gas molecules and liquid molecules currently has no very effective method:
for gas molecules, due to the extremely complex potential energy surface, the gas molecules are easy to fall into the local energy lowest point during simulation, thereby influencing the convergence of the free energy calculation result.
For liquid molecules, the situation is more complicated, and not only the influence of the potential energy surface of the liquid molecules is needed to be overcome, but also the sampling problem caused by the rapid diffusion of the molecules is also overcome. Therefore, it is more difficult to accurately calculate the absolute free energy of the liquid molecules.
Disclosure of Invention
In view of the above technical problems, the present invention provides a method for accurately calculating the absolute free energy of gas molecules, which can accurately calculate the absolute free energy of solid and gas, and can indirectly calculate the absolute free energy of liquid molecules by calculating the solvation free energy.
In order to achieve the purpose, the invention provides the following technical scheme:
a method for accurately calculating the absolute free energy of a gas molecule, comprising the steps of:
s1, reducing the energy barrier of the flexible dihedral angle potential energy surface in the gas molecules;
s1.1, preparing a coordinate file of a molecule to be calculated, and grabbing force field parameters through a related tool;
s1.2, calling yoda program, loading molecular coordinates and a force field, and searching all flexible dihedral angles in molecules;
s1.3, potential energy surface scanning is carried out on all the found flexible dihedral angles;
s1.4, fitting a dihedral angle potential energy surface with an energy barrier larger than 5kcal/mol in the potential energy surface, modifying corresponding dihedral angle parameters, and reducing the energy barrier;
s1.5 repeating steps S1.3 and S1.4 until no energy barrier greater than 5kcal/mol exists in the potential energy planes of all compliant angles;
s1.6, generating a new force field parameter file;
s2, calculating absolute free energy of gas molecules;
s2.1, establishing a simulation box with the size of 20nm multiplied by 20nm by using new force field parameters, carrying out NVT simulation at a series of temperatures, and calculating the electrostatic acting force by using a cut-off mode, wherein the step length is 1fs, and the simulation is 10 ns;
s2.2, selecting two balance structures at reference temperature as initial structures of free energy simulation;
s2.3 stepwise limiting the position of a central heavy atom, the force constant varying from 0 to 10000kJ/mol/nm2The step is difficult to converge by calculating the free energy through MD simulation, but can be derived from theory;
s2.4 stepwise limiting the position of other heavy atoms, the force constant varying from 0 to 10000kJ/mol/nm2
S2.5 stepwise limiting the position of the non-heavy atoms, the force constant varying from 0 to 10000kJ/mol/nm2
S2.6 stepwise increasing the position restriction of all atoms, the force constant varying from 10000 to 200000kJ/mol/nm2
S2.7 step-by-step removal of all atomic intramolecular (electrostatic, van der waals, bond, angle, dihedral) interactions;
s2.8, after the simulation is normal, all using the original force field parameters to rerun the track to obtain the unbound potential energy;
s2.9, loading the binary and unbiased energy of each lambda or T window by using the modified WHAM program to obtain the final free energy after reweigh;
s2.10 uses the reweight Δ f (T) and PSCP free energy and related theoretical derivation values to calculate the final absolute free energy change with temperature.
Wherein, the related tool in step S1.1 is antechamber or CHARMM-GUI.
The series of temperatures in the step S2.1 are within the range of 100-400K.
Step S2.2 selects two reference temperatures of 200K and 300K.
Compared with the prior art, the invention has the beneficial effects that:
(1) and the automatic identification of flexible dihedral angles in molecules and the potential energy surface scanning are realized.
(2) Fitting the dihedral angle potential energy surface is realized, and the dihedral angle parameters are automatically modified so as to reduce the potential energy surface to a specific value.
(3) A new WHAM program is realized, and the actual free energy can be obtained from the potential energies reweigh of biased and unbound.
(4) A set of feasible calculation process is designed to realize the gradual change of potential energy when the gas molecules gradually remove various energies.
(5) All the processes are connected together, and the automatic calculation of the absolute free energy of the gas molecules is realized.
Drawings
FIG. 1 is a schematic representation of the molecular structure of butane in accordance with the examples;
FIG. 2 is a dihedral angle potential profile of an example butane molecule initial C _ C1_ C2_ C3;
FIG. 3 is a potential energy surface obtained by scanning according to an embodiment;
FIG. 4 is a flow chart of the calculation of absolute free energy of gas molecules according to the present invention;
FIG. 5 is an absolute free energy curve obtained using an original force field;
FIG. 6 is a plot of the absolute free energy obtained by the WHAM program reweight using the new force field in this example.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
This example uses the calculation procedure of the absolute free energy of butane (butane) molecule, as shown in FIG. 1:
(1) preparing a coordinate (pdb) file of butane molecules, and grabbing corresponding GAFF2 force field parameters by using antechamber;
(2) calling yoda program, loading molecular coordinates and force field, and searching all flexible dihedral angles in molecules;
(3) potential energy surface scanning is carried out on all the found flexible dihedral angles;
the butane molecule has only one flexible dihedral angle (C _ C1_ C2_ C3) as shown in fig. 2.
(4) Fitting a dihedral angle potential surface with an energy barrier larger than 5kcal/mol in the potential surface, modifying corresponding dihedral angle parameters, and reducing the energy barrier;
fitting function:
f(x)=k1*(1+cosx)+k2*(1-cos2x)+k3*(1+cos3x)+k4*(1+cos4x)+k6*(1+cos6x)+b
(5) repeating steps (3) and (4) until no energy barrier greater than 5kcal/mol exists in the potential energy surfaces of all the flexible angles; after modifying the dihedral angle parameter C _ C1_ C2_ C3, scanning the obtained potential energy surface, as shown in FIG. 3;
(6) generating a new force field parameter file;
the flow of calculation of the absolute free energy of the gas molecules is shown in FIG. 4.
(7) Establishing a 20nm multiplied by 20nm simulated box by using new force field parameters, carrying out NVT simulation under a series of temperatures (100-400K), and calculating the electrostatic acting force by using a cut-off mode, wherein the step size is 1fs and the simulation is 10 ns;
through the potential energy data H (T) obtained by simulation at a series of temperatures, the non-dimensional change of the free energy along with the temperature Δ f (T) can be obtained by integration:
Figure BDA0002738364950000031
(8) selecting a balance structure under two reference temperatures (200K and 300K for convergence verification) as an initial structure of free energy simulation;
(9) gradually limiting the position of a central heavy atom, and changing the force constant from 0 to 10000kJ/mol/nm2The step is difficult to converge by calculating the free energy through MD simulation, but can be derived from theory (Δ A)1);
(10) The positions of other heavy atoms are gradually limited, and the force constant is changed from 0 to 10000kJ/mol/nm2(ΔA2);
(11) The position of non-heavy atoms is gradually limited, and the force constant is changed from 0 to 10000kJ/mol/nm2(ΔA3);
(12) The position limits of all atoms are gradually increased, and the force constant is changed from 10000 to 200000kJ/mol/nm2(ΔA4);
(13) Intramolecular (electrostatic, van der waals, bond, angle, dihedral) interactions (Δ a) that gradually remove all atoms5);
(14) After the simulation is normal, all the original force field parameters are used for rerun on the track to obtain the unbound potential energy;
the total PSCP free energy change (. DELTA.G) can be calculated by the above simulationPSCP(T))
Figure BDA0002738364950000041
(15) Loading the binary and unbiased energy of each lambda or T window by using the modified WHAM program to obtain the final free energy after reweigh;
the principle of the WHAM program for reweight:
calculation formula of standard WHAM:
Figure BDA0002738364950000042
Figure BDA0002738364950000043
Figure BDA0002738364950000044
(16) the change of the final absolute free energy with the temperature is calculated by using the delta f (T) and PSCP free energy after reweight and related theoretical derivation values.
The final absolute free energy is calculated as follows:
Figure BDA0002738364950000045
the results are shown in FIGS. 5 and 6.
From the calculation results, it can be seen that: the result of the new force field is combined with the free energy after reweight, and the convergence is obviously improved. This increase is more pronounced on more complex molecules.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.

Claims (6)

1. A method for accurately calculating the absolute free energy of a gas molecule is characterized by comprising the following steps:
s1, reducing the energy barrier of the flexible dihedral angle potential energy surface in the gas molecules;
and S2, calculating absolute free energy of gas molecules.
2. The method for accurately calculating the absolute free energy of gas molecules as claimed in claim 1, wherein the step S1 comprises the following sub-steps:
s1.1, preparing a coordinate file of a molecule to be calculated, and grabbing force field parameters through a related tool;
s1.2, calling yoda program, loading molecular coordinates and a force field, and searching all flexible dihedral angles in molecules;
s1.3, potential energy surface scanning is carried out on all the found flexible dihedral angles;
s1.4, fitting a dihedral angle potential energy surface with an energy barrier larger than 5kcal/mol in the potential energy surface, modifying corresponding dihedral angle parameters, and reducing the energy barrier;
s1.5 repeating steps S1.3 and S1.4 until no energy barrier greater than 5kcal/mol exists in the potential energy planes of all compliant angles;
s1.6, generating a new force field parameter file.
3. A method for accurately calculating the absolute free energy of a gas molecule according to claim 2, wherein the relevant tool in step S1.1 is antechamber or CHARMM-GUI.
4. The method for accurately calculating the absolute free energy of gas molecules as claimed in claim 1, wherein the step S2 comprises the following sub-steps:
s2.1, establishing a simulation box with the size of 20nm multiplied by 20nm by using new force field parameters, carrying out NVT simulation at a series of temperatures, and calculating the electrostatic acting force by using a cut-off mode, wherein the step length is 1fs, and the simulation is 10 ns;
s2.2, selecting two balance structures at reference temperature as initial structures of free energy simulation;
s2.3 stepwise limiting the position of a central heavy atom, the force constant varying from 0 to 10000kJ/mol/nm2Derived from theory;
s2.4 stepwise limiting the position of other heavy atoms, the force constant varying from 0 to 10000kJ/mol/nm2
S2.5 stepwise limiting the position of the non-heavy atoms, the force constant varying from 0 to 10000kJ/mol/nm2
S2.6 stepwise increasing the position restriction of all atoms, the force constant varying from 10000 to 200000kJ/mol/nm2
S2.7 removing the intramolecular interactions of all atoms step by step;
s2.8, after the simulation is normal, all the original force field parameters are used for rerun on the track to obtain the unbound potential energy;
s2.9, loading the binary and unbiased energy of each lambda or T window by using the modified WHAM program to obtain the final free energy after reweigh;
s2.10 uses the reweigh Δ f (T) and PSCP free energies, and theoretical derivation values, to calculate the final absolute free energy change with temperature.
5. The method of claim 4, wherein the series of temperatures in step S2.1 is in the range of 100-400K.
6. The method according to claim 4, wherein the two reference temperatures of 200K and 300K are selected in step S2.2.
CN202011141320.0A 2020-10-22 2020-10-22 Method for accurately calculating absolute free energy of gas molecules Pending CN112199909A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011141320.0A CN112199909A (en) 2020-10-22 2020-10-22 Method for accurately calculating absolute free energy of gas molecules

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011141320.0A CN112199909A (en) 2020-10-22 2020-10-22 Method for accurately calculating absolute free energy of gas molecules

Publications (1)

Publication Number Publication Date
CN112199909A true CN112199909A (en) 2021-01-08

Family

ID=74012344

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011141320.0A Pending CN112199909A (en) 2020-10-22 2020-10-22 Method for accurately calculating absolute free energy of gas molecules

Country Status (1)

Country Link
CN (1) CN112199909A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114187971A (en) * 2021-12-10 2022-03-15 上海智药科技有限公司 Molecular free energy calculation and stability analysis method, device, equipment and storage medium
WO2023102870A1 (en) * 2021-12-10 2023-06-15 上海智药科技有限公司 Method and apparatus for computing free energy of molecule, and stability analysis method and apparatus, and device and storage medium

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030195734A1 (en) * 2002-04-16 2003-10-16 Huai Sun Method and expert system of molecular mechanics force fields for computer simulation of molecular systems
US20150317459A1 (en) * 2012-12-11 2015-11-05 Asaf FARHI Method to calculate free energies
JP2016181257A (en) * 2015-03-24 2016-10-13 東レ株式会社 Calculation method of pressures generated by van der waals interactions, and sampling method of high molecule-low molecule interaction energy or arrangement of low molecules in high molecules
CN109859806A (en) * 2019-01-17 2019-06-07 中山大学 A kind of prediction drug-target bond strength absolute freedom energy perturbation method
CN111161810A (en) * 2019-12-31 2020-05-15 中山大学 Free energy perturbation method based on constraint probability distribution function optimization
CN111415710A (en) * 2020-03-06 2020-07-14 深圳晶泰科技有限公司 Potential energy surface scanning method and system for molecular conformation space analysis
CN111653320A (en) * 2020-04-21 2020-09-11 深圳晶泰科技有限公司 Molecular force field fitting method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030195734A1 (en) * 2002-04-16 2003-10-16 Huai Sun Method and expert system of molecular mechanics force fields for computer simulation of molecular systems
US20150317459A1 (en) * 2012-12-11 2015-11-05 Asaf FARHI Method to calculate free energies
JP2016181257A (en) * 2015-03-24 2016-10-13 東レ株式会社 Calculation method of pressures generated by van der waals interactions, and sampling method of high molecule-low molecule interaction energy or arrangement of low molecules in high molecules
CN109859806A (en) * 2019-01-17 2019-06-07 中山大学 A kind of prediction drug-target bond strength absolute freedom energy perturbation method
CN111161810A (en) * 2019-12-31 2020-05-15 中山大学 Free energy perturbation method based on constraint probability distribution function optimization
CN111415710A (en) * 2020-03-06 2020-07-14 深圳晶泰科技有限公司 Potential energy surface scanning method and system for molecular conformation space analysis
CN111653320A (en) * 2020-04-21 2020-09-11 深圳晶泰科技有限公司 Molecular force field fitting method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"2_分子力场方法.pdf", pages 1 - 87, Retrieved from the Internet <URL:https://max.book118.com/html/2019/0814/8122126044002042.shtm> *
J. KALIKKA ETC.: "Density functional and classical simulations of liquid and glassy selenium", PHYSICAL REVIEW B, vol. 102, no. 10, pages 1 - 13 *
LI WENTAO ETC.: "Free energy study of H2O, N2O5, SO2, and O3 gas sorption", THE JOURNAL OF CHEMICAL PHYSICS, vol. 148, no. 16, pages 1 - 8 *
LI, HONGZHI ETC.: "Forging the missing link in free energy estimations λ-WHAM in thermodynamic integration, overlap histogramming, and free energy perturbation", CHEMICAL PHYSICS LETTERS, vol. 440, no. 1, pages 155 - 159, XP022067602, DOI: 10.1016/j.cplett.2007.04.009 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114187971A (en) * 2021-12-10 2022-03-15 上海智药科技有限公司 Molecular free energy calculation and stability analysis method, device, equipment and storage medium
WO2023102870A1 (en) * 2021-12-10 2023-06-15 上海智药科技有限公司 Method and apparatus for computing free energy of molecule, and stability analysis method and apparatus, and device and storage medium

Similar Documents

Publication Publication Date Title
CN112199909A (en) Method for accurately calculating absolute free energy of gas molecules
Delogu Structural and energetic properties of unsupported Cu nanoparticles from room temperature to the melting point: molecular dynamics simulations
CN111666615B (en) Suspension bridge space cable shape finding method based on finite element
CN108237531B (en) Humanoid robot gait self-learning control method
CN106022212B (en) A kind of gyro Temperature Drift Modeling
Shi et al. Atomic diffusion mediated by vacancy defects in pure and transition element (TM)-doped (TM= Ti, Y, Zr or Hf) L12 Al3Sc
Vopilov et al. Effect of supercritical carbon dioxide on ultradispersed polytetrafluoroethylene
CN113223631A (en) Crystal structure analysis method and system
CN107944116B (en) A kind of efficient predicting method of transient energy response for Time variable structure
WO2017152397A1 (en) Method for testing phase transformation point of aluminium alloy
WO2022082598A1 (en) Method for accurately calculating absolute free energy of gas molecules
CN109583030B (en) Mixed wave height and period combined distribution construction method based on mixed model
CN108182340B (en) Rapid screening method of high solid solubility multi-principal element alloy
Karkin et al. SEGREGATION OF Mg TO GENERIC TILT GRAIN BOUNDARIES IN Al: MONTE CARLO MODELING.
CN106777825B (en) Rectangular plate vibration mode calculation method based on spectrum finite element
CN111428190B (en) Real-time hybrid simulation self-adaptive compensation method based on frequency domain evaluation index
CN108277506A (en) A kind of adaptive oversampler method of aluminium electroloysis degree of superheat unbalanced data
CN106503394B (en) A method of establishing the fractional order circuit model of tantalum capacitor
CN113392372A (en) Gas well full-life-cycle unobstructed flow determination method and system
CN105910886B (en) The application of the Browman of aluminium alloy stress-strain relation this structure Optimized models
CN112084455B (en) Least square interpolation method for S parameter
Krūkle‐Bērziņa et al. Powder X‐Ray Diffraction Investigation of Xylazine Hydrochloride Solid Phase Transformation Kinetics
Zhang et al. Modified Lennard-Jones potentials for Cu and Ag based on the dense gaslike model of viscosity for liquid metals
Rincent et al. On the exploration of the melting behavior of metallic compounds and solid solutions via multiple classical molecular dynamics approaches: application to Al-based systems
CN106980731B (en) Method for calculating mass fraction of structural component of iron-carbon alloy

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination