Accurate and efficient integration of the equations of motion is indispensable for molecular dynamics(MD)simulations.Despite the massive use of the conventional leapfrog(LF)integrator in modern computational tools wit...Accurate and efficient integration of the equations of motion is indispensable for molecular dynamics(MD)simulations.Despite the massive use of the conventional leapfrog(LF)integrator in modern computational tools within the framework of MD propagation,further development for better performance is still possible.The alternative version of LF in the middle thermostat scheme(LFmiddle)achieves a higher order of accuracy and efficiency and maintains stable dynamics even with the integration time stepsize extended by several folds.In this work,we perform a benchmark test of the two integrators(LF and LF-middle)in extensive conventional and enhanced sampling simulations,aiming at quantifying the time-stepsizeinduced variations of global properties(e.g.,detailed potential energy terms)as well as of local observables(e.g.,free energy changes or bondlengths)in practical simulations of complex systems.The test set is composed of six chemically and biologically relevant systems,including the conformational change of dihedral flipping in the N-methylacetamide and an AT(AdenineThymine)tract,the intra-molecular proton transfer inside malonaldehyde,the binding free energy calculations of benzene and phenol targeting T4 lysozyme L99A,the hydroxyl bond variations in ethaline deep eutectic solvent,and the potential energy of the blue-light using flavin photoreceptor.It is observed that the time-step-induced error is smaller for the LFmiddle scheme.The outperformance of LF-middle over the conventional LF integrator is much more significant for global properties than local observables.Overall,the current work demonstrates that the LF-middle scheme should be preferably applied to obtain accurate thermodynamics in the simulation of practical chemical and biological systems.展开更多
Is it possible to improve CO_(2)solubility in potassium carbonate(K_(2)CO_(3))-based transition temperature mixtures(TTMs)?To assess this possibility,a ternary transition-temperature mixture(TTTM)was prepared by using...Is it possible to improve CO_(2)solubility in potassium carbonate(K_(2)CO_(3))-based transition temperature mixtures(TTMs)?To assess this possibility,a ternary transition-temperature mixture(TTTM)was prepared by using a hindered amine,2-amino-2-methyl-1,3-propanediol(AMPD).Fourier transform infrared spectroscopy(FT-IR)was employed to detect the functional groups including hydroxyl,amine,carbonate ion,and aliphatic functional groups in the prepared solvents.From thermogravimetric analysis(TGA),it was found that the addition of AMPD to the binary mixture can increase the thermal stability of TTTM.The viscosity findings showed that TTTM has a higher viscosity than TTM while their difference was decreased by increasing temperature.In addition,Eyring’s absolute rate theory was used to compute the activation parameters(∆G^(*),∆H^(*),and ∆S^(*)).The CO_(2)solubility in liquids was measured at a temperature of 303.15 K and pressures up to 1.8 MPa.The results disclosed that the CO_(2)solubility of TTTM was improved by the addition of AMPD.At the pressure of about 1.8 MPa,the CO_(2)mole fractions of TTM and TTTM were 0.1697 and 0.2022,respectively.To confirm the experimental data,density functional theory(DFT)was employed.From the DFT analysis,it was found that the TTTM+CO_(2)system has higher interaction energy(|∆E|)than the TTM+CO_(2)system indicating the higher CO_(2)affinity of the former system.This study might help scientists to better understand and to improve CO_(2)solubility in these types of solvents by choosing a suitable amine as HBD and finding the best combination of HBA and HBD.展开更多
基金supported by the National Natural Science Foundation of China(No.21961142017)the Ministry of Science and Technology of China(No.2017YFA0204901)。
文摘Accurate and efficient integration of the equations of motion is indispensable for molecular dynamics(MD)simulations.Despite the massive use of the conventional leapfrog(LF)integrator in modern computational tools within the framework of MD propagation,further development for better performance is still possible.The alternative version of LF in the middle thermostat scheme(LFmiddle)achieves a higher order of accuracy and efficiency and maintains stable dynamics even with the integration time stepsize extended by several folds.In this work,we perform a benchmark test of the two integrators(LF and LF-middle)in extensive conventional and enhanced sampling simulations,aiming at quantifying the time-stepsizeinduced variations of global properties(e.g.,detailed potential energy terms)as well as of local observables(e.g.,free energy changes or bondlengths)in practical simulations of complex systems.The test set is composed of six chemically and biologically relevant systems,including the conformational change of dihedral flipping in the N-methylacetamide and an AT(AdenineThymine)tract,the intra-molecular proton transfer inside malonaldehyde,the binding free energy calculations of benzene and phenol targeting T4 lysozyme L99A,the hydroxyl bond variations in ethaline deep eutectic solvent,and the potential energy of the blue-light using flavin photoreceptor.It is observed that the time-step-induced error is smaller for the LFmiddle scheme.The outperformance of LF-middle over the conventional LF integrator is much more significant for global properties than local observables.Overall,the current work demonstrates that the LF-middle scheme should be preferably applied to obtain accurate thermodynamics in the simulation of practical chemical and biological systems.
基金support of National Key Research and Development Program of China(Grant No.2019YFC1904602)the Key SCI-Tech Innovation 2025 Program of Ningbo,China(Grant No.2018B10025).
文摘Is it possible to improve CO_(2)solubility in potassium carbonate(K_(2)CO_(3))-based transition temperature mixtures(TTMs)?To assess this possibility,a ternary transition-temperature mixture(TTTM)was prepared by using a hindered amine,2-amino-2-methyl-1,3-propanediol(AMPD).Fourier transform infrared spectroscopy(FT-IR)was employed to detect the functional groups including hydroxyl,amine,carbonate ion,and aliphatic functional groups in the prepared solvents.From thermogravimetric analysis(TGA),it was found that the addition of AMPD to the binary mixture can increase the thermal stability of TTTM.The viscosity findings showed that TTTM has a higher viscosity than TTM while their difference was decreased by increasing temperature.In addition,Eyring’s absolute rate theory was used to compute the activation parameters(∆G^(*),∆H^(*),and ∆S^(*)).The CO_(2)solubility in liquids was measured at a temperature of 303.15 K and pressures up to 1.8 MPa.The results disclosed that the CO_(2)solubility of TTTM was improved by the addition of AMPD.At the pressure of about 1.8 MPa,the CO_(2)mole fractions of TTM and TTTM were 0.1697 and 0.2022,respectively.To confirm the experimental data,density functional theory(DFT)was employed.From the DFT analysis,it was found that the TTTM+CO_(2)system has higher interaction energy(|∆E|)than the TTM+CO_(2)system indicating the higher CO_(2)affinity of the former system.This study might help scientists to better understand and to improve CO_(2)solubility in these types of solvents by choosing a suitable amine as HBD and finding the best combination of HBA and HBD.
