1. |
- Abraham, Mark J, et al.
(författare)
-
Optimization of parameters for molecular dynamics simulation using smooth particle-mesh Ewald in GROMACS 4.5
- 2011
-
Ingår i: Journal of Computational Chemistry. - : John Wiley & Sons. - 0192-8651 .- 1096-987X. ; 32:9
-
Tidskriftsartikel (refereegranskat)abstract
- Based on our critique of requirements for performing an efficient molecular dynamics simulation with the particle-mesh Ewald (PME) implementation in GROMACS 4.5, we present a computational tool to enable the discovery of parameters that produce a given accuracy in the PME approximation of the full electrostatics. Calculations on two parallel computers with different processor and communication structures showed that a given accuracy can be attained over a range of parameter space, and that the attributes of the hardware and simulation system control which parameter sets are optimal. This information can be used to find the fastest available PME parameter sets that achieve a given accuracy. We hope that this tool will stimulate future work to assess the impact of the quality of the PME approximation on simulation outcomes, particularly with regard to the trade-off between cost and scientific reliability in biomolecular applications.
|
|
2. |
- Abraham, Mark J
(författare)
-
Performance enhancements for GROMACS nonbonded interactions on BlueGene.
- 2011
-
Ingår i: Journal of Computational Chemistry. - : John Wiley & Sons. - 0192-8651 .- 1096-987X. ; 32:9
-
Tidskriftsartikel (refereegranskat)abstract
- Several improvements to the previously optimized GROMACS BlueGene inner loops that evaluate nonbonded interactions in molecular dynamics simulations are presented. The new improvements yielded an 11% decrease in running time for both PME and other kinds of GROMACS simulations that use nonbonded table look-ups. Some other GROMACS simulations will show a small gain.
|
|
3. |
- Andersson, Rasmus, 1990, et al.
(författare)
-
CHAMPION: Chalmers hierarchical atomic, molecular, polymeric and ionic analysis toolkit
- 2021
-
Ingår i: Journal of Computational Chemistry. - : Wiley. - 0192-8651 .- 1096-987X. ; 42:23, s. 1632-1642
-
Tidskriftsartikel (refereegranskat)abstract
- We present CHAMPION (Chalmers hierarchical atomic, molecular, polymeric, and ionic analysis toolkit): a software developed to automatically detect time-dependent bonds between atoms based on their dynamics, classify the local graph topology around them, and analyze the physicochemical properties of these topologies by statistical physics. In stark contrast to methodologies where bonds are detected based on static conditions such as cut-off distances, CHAMPION considers pairs of atoms to be bound only if they move together and act as a bound pair over time. Furthermore, the time-dependent global bond graph is possible to split into dynamically shifting connected components or subgraphs around a certain chemical motif and thereby allow the physicochemical properties of each such topology to be analyzed by statistical physics. Applicable to condensed matter and liquids in general, and electrolytes in particular, this allows both quantitative and qualitative descriptions of local structure, as well as dynamical processes such as speciation and diffusion. We present here a detailed overview of CHAMPION, including its underlying methodology, implementation, and capabilities.
|
|
4. |
- Andres, Juan, et al.
(författare)
-
Nine questions on energy decomposition analysis
- 2019
-
Ingår i: Journal of Computational Chemistry. - : Wiley. - 0192-8651 .- 1096-987X. ; 40:26, s. 2248-2283
-
Tidskriftsartikel (refereegranskat)abstract
- The paper collects the answers of the authors to the following questions: Is the lack of precision in the definition of many chemical concepts one of the reasons for the coexistence of many partition schemes? Does the adoption of a given partition scheme imply a set of more precise definitions of the underlying chemical concepts? How can one use the results of a partition scheme to improve the clarity of definitions of concepts? Are partition schemes subject to scientific Darwinism? If so, what is the influence of a community's sociological pressure in the “natural selection” process? To what extent does/can/should investigated systems influence the choice of a particular partition scheme? Do we need more focused chemical validation of Energy Decomposition Analysis (EDA) methodology and descriptors/terms in general? Is there any interest in developing common benchmarks and test sets for cross-validation of methods? Is it possible to contemplate a unified partition scheme (let us call it the “standard model” of partitioning), that is proper for all applications in chemistry, in the foreseeable future or even in principle? In the end, science is about experiments and the real world. Can one, therefore, use any experiment or experimental data be used to favor one partition scheme over another?
|
|
5. |
- Aquilante, Francesco, et al.
(författare)
-
Molcas 8 : New capabilities for multiconfigurational quantum chemical calculations across the periodic table
- 2016
-
Ingår i: Journal of Computational Chemistry. - : Wiley. - 0192-8651 .- 1096-987X. ; 37:5, s. 506-541
-
Tidskriftsartikel (refereegranskat)abstract
- In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.
|
|
6. |
- Aquilante, Francesco, et al.
(författare)
-
Software news and update MOLCAS 7 : The Next Generation
- 2010
-
Ingår i: Journal of Computational Chemistry. - : Wiley. - 0192-8651 .- 1096-987X. ; 31:1, s. 224-247
-
Forskningsöversikt (refereegranskat)abstract
- Some of the new unique features of the MOLCAS quantum chemistry package version 7 are presented in this report. In particular, the Cholesky decomposition method applied to some quantum chemical methods is described. This approach is used both in the context of a straight forward approximation of the two-electron integrals and in the generation of so-called auxiliary basis sets. The article describes how the method is implemented for most known wave functions models: self-consistent field, density functional theory, 2nd order perturbation theory, complete-active space self-consistent field multiconfigurational reference 2nd order perturbation theory, and coupled-cluster methods. The report further elaborates on the implementation of a restricted-active space self-consistent field reference function in conjunction with 2nd order perturbation theory. The average atomic natural orbital basis for relativistic calculations, covering the whole periodic table, are described and associated unique properties are demonstrated. Furthermore, the use of the arbitrary order Douglas-Kroll-Hess transformation for one-component relativistic calculations and its implementation are discussed. This section especially focuses on the implementation of the so-called picture-change-free atomic orbital property integrals. Moreover, the ElectroStatic Potential Fitted scheme, a version of a quantum mechanics/molecular mechanics hybrid method implemented in MOLCAS, is described and discussed. Finally, the report discusses the use of the MOLCAS package for advanced studies of photo chemical phenomena and the usefulness of the algorithms for constrained geometry optimization in MOLCAS in association with such studies.
|
|
7. |
- Blomberg, Margareta R. A., et al.
(författare)
-
Improved free energy profile for reduction of NO in cytochrome c dependent nitric oxide reductase (cNOR)
- 2016
-
Ingår i: Journal of Computational Chemistry. - : Wiley. - 0192-8651 .- 1096-987X. ; 37:19, s. 1810-1818
-
Tidskriftsartikel (refereegranskat)abstract
- Quantum chemical calculations play an essential role in the elucidation of reaction mechanisms for redox-active metalloenzymes. For example, the cleavage and the formation of covalent bonds can usually not be described only on the basis of experimental information, but can be followed by the calculations. Conversely, there are properties, like reduction potentials, which cannot be accurately calculated. Therefore, computational and experimental data has to be carefully combined to obtain reliable descriptions of entire catalytic cycles involving electron and proton uptake from donors outside the enzyme. Such a procedure is illustrated here, for the reduction of nitric oxide (NO) to nitrous oxide and water in the membrane enzyme, cytochrome c dependent nitric oxide reductase (cNOR). A surprising experimental observation is that this reaction is nonelectrogenic, which means that no energy is conserved. On the basis of hybrid density functional calculations a free energy profile for the entire catalytic cycle is obtained, which agrees much better with experimental information on the active site reduction potentials than previous ones. Most importantly the energy profile shows that the reduction steps are endergonic and that the entire process is rate-limited by high proton uptake barriers during the reduction steps. This result implies that, if the reaction were electrogenic, it would become too slow when the gradient is present across the membrane. This explains why this enzyme does not conserve any of the free energy released.
|
|
8. |
- Blomgren, Fredrik, 1973, et al.
(författare)
-
Exploring the potential energy surface of retinal, a comparison of the performance of different methods
- 2005
-
Ingår i: Journal of Computational Chemistry. - : Wiley. - 0192-8651 .- 1096-987X. ; 26:7, s. 738-742
-
Tidskriftsartikel (refereegranskat)abstract
- The ground state structure of retinal has been investigated. We found that DFT and CASSCF produce different results for the bond length alternation in a model system of retinal. Quantum mechanics/molecular mechanics calculations including the closest surrounding amino acids have been performed, using DFT and CASSCF to calculate the structure of retinal in the protein cavity. The planarity of the retinal molecule is affected by the surrounding protein. DFT and CASSCF produce different twist angles. The difference between CASSCF and DFT appears to be related to the positively charged nitrogen of the Schiff base, which leads to different π-bond orders produced by the two methods.
|
|
9. |
- Bondesson, Laban, et al.
(författare)
-
Basis set dependence of solute-solvent interaction energy of benzene in water : a HF/DFT study
- 2008
-
Ingår i: Journal of Computational Chemistry. - : Wiley. - 0192-8651 .- 1096-987X. ; 29:11, s. 1725-1732
-
Tidskriftsartikel (refereegranskat)abstract
- Solute-solvent interaction energies for the benzene molecule dissolved in water are computed using Hartree-Fock and B3LYP density functional theories. Explicit solvent molecules up to 14-angstrom away from the dissolved benzene molecule are included in the calculation of interaction energies. Both basis set dependence and basis Set Superposition errors are carefully examined. It is found that the use of a larger basis set for the region near the solute together with a smaller basis set for the outer region gives results very close to what would have been obtained if the larger basis set had been used for the whole system. It is also shown that a correction for the basis Set superposition error is a necessary component in this kind of calculations. With this correction, results obtained with different tested basis sets converge to within 1 kcal/mol.
|
|
10. |
|
|