Molecular dynamics simulations of matrix deposition. III. Site structure analysis for porphycene in argon and xenon

Porphycene (1) and porphyrin (2), two constitutional isomers, reveal completely different electronic spectral patterns in argon and xenon matrices. For the former the spectra recorded in the two solidified gases resemble each other, whereas for the latter they are completely different. This difference can be rationalized by molecular-dynamics simulations of the structure of the microenvironment carried out for the two chromophores embedded in argon and xenon hosts. For 1, the structure of the main substitutional site is the same for Ar and Xe and consists of a hexagonal cavity obtained by removing seven host atoms from the [111] crystallographic plane. An analogous structure is obtained for 2 in xenon. However, in argon the porphyrin chromophore environment is shared between several different sites, with the number of replaced host atoms ranging from seven to ten. These results demonstrate that a relatively minor structural alternation may lead to major changes in the spectral pattern of molecules embedded in rare-gas cryogenic matrices.     

Molecular dynamics and density functional theory simulations of cesium and strontium adsorption on illite/ smectite

Journal of Radioanalytical and Nuclear Chemistry - Based on molecular dynamics (MD) and density function theory (DFT) simulation, the adsorption mechanisms of Cs+ and Sr2+ on the...     

Addressing an instability in unrestricted density functional theory direct dynamics simulations

In Density Functional Theory (DFT) direct dynamics simulations with Unrestricted Hartree Fock (UHF) theory, triplet instability often emerges when numerically integrating a classical trajectory. A broken symmetry initial guess for the wave function is often used to obtain the unrestricted DFT potential energy surface (PES), but this is found to be often insufficient for direct dynamics simulations. An algorithm is described for obtaining smooth transitions between the open-shell and the closed-shell regions of the unrestricted PES, and thus stable trajectories, for direct dynamics simulations of dioxetane and its •O CH₂-CH₂ O• singlet diradical. © 2018 Wiley Periodicals, Inc.     

Real-time, local basis-set implementation of time-dependent density functional theory for excited state dynamics simulations

We present a method suitable for large-scale accurate simulations of excited state dynamics within the framework of time-dependent density functional theory (DFT). This is achieved by employing a local atomic basis-set representation and real-time propagation of excited state wave functions. We implement the method within SIESTA, a standard ground-state DFT package with local atomic basis, and demonstrate its potential for realistic and accurate excited state dynamics simulations using small and medium-sized molecules as examples (H(2), CO, O(3), and indolequinone). The method can be readily applied to problems involving nanostructures and large biomolecules.     

The structure of the second-order reduced density matrix in density matrix functional theory and its construction from formal criteria

Some formal requirements for the second-order reduced density matrix are discussed in the context of density matrix functional theory. They serve as a basis for the ad hoc construction of the second-order reduced density matrix in terms of the first-order reduced density matrix and lead to implicit functionals where the occupation numbers of the natural orbitals are obtained as diagonal elements of an idempotent matrix the elements of which represent the variational parameters to be optimized. The numerical results obtained from a first realization of such an implicit density matrix functional give excellent agreement with the results of full configuration interaction calculations for four-electron systems like LiH and Be. Results for H2O taken as an example for a somewhat larger molecule are numerically less satisfactory but still give reasonable occupation numbers of the natural orbitals and indicate the capability of density matrix functional theory to cope with static electron correlation.     

Towards an assessment of the accuracy of density functional theory for first principles simulations of water. II

A series of 20 ps ab initio molecular dynamics simulations of water at ambient density and temperatures ranging from 300 to 450 K are presented. Car-Parrinello (CP) and Born-Oppenheimer (BO) molecular dynamics techniques are compared for systems containing 54 and 64 water molecules. At 300 K, an excellent agreement is found between radial distribution functions (RDFs) obtained with BO and CP dynamics, provided an appropriately small value of the fictitious mass parameter is used in the CP simulation. However, we find that the diffusion coefficients computed from CP dynamics are approximately two times larger than those obtained with BO simulations for T>400 K, where statistically meaningful comparisons can be made. Overall, both BO and CP dynamics at 300 K yield overstructured RDFs and slow diffusion as compared to experiment. In order to understand these discrepancies, the effect of proton quantum motion is investigated with the use of empirical interaction potentials. We find that proton quantum effects may have a larger impact than previously thought on structure and diffusion of the liquid.     

Vibrational assignment and structure of dibenzoylmethane. A density functional theoretical study

Molecular structure and vibrational frequencies of 1,3-diphenyl-1,3-propanedione, known as dibenzoylmethane (DBM), have been investigated by means of density functional theory (DFT) calculations. The results were compared with those of benzoylacetone (BA) and acetylacetone (AA), the parent molecule. IR and Raman spectra of DBM and its deuterated analogue were clearly assigned.The calculated hydrogen bond energy of DBM is 16.15 kcal/mol, calculated at B3LYP/6-311++G** level of theory, which is 0.28 kcal/mol more than that of AA. This result is in agreement with the vibrational and NMR spectroscopy results. The molecular stability and the hydrogen bond strength were investigated by applying the Natural Bond Orbital analysis (NBO) and geometry calculations. The theoretical calculations indicate that the hydrogen bond in DBM is relatively stronger than that in BA and AA.     

Vibrational assignment and structure of trifluorobenzoylacetone. A density functional theoretical study

Molecular structure and vibrational frequencies of 4,4,4-trifluoro-1-phenyl-1,3-butanedione, known as trifluorobenzoylacetone (TFBA), have been investigated by means of density functional theory (DFT) calculations. The results were compared with those of benzoylacetone (BA), acetylacetone (AA), and trifluoroacetylacetone (TFAA). Comparing the calculated and experimental band frequencies and intensities suggests coexisting of both stable cis-enol conformers in comparable proportions in the sample. The energy difference between the two stable chelated enol forms is negligible, 0.96 kcal/mol, calculated at B3LYP/6-311++G** level of theory. The molecular stability and the hydrogen bond strength were investigated by applying the natural bond orbital (NBO) theory and geometry calculations. The theoretical calculations and spectroscopic results indicate that the hydrogen bond strength of TFBA is between those of TFAA and AA, considerably weaker than that of BA.     

Molecular dynamics simulations of defect formation in thin graphite films using the density functional tight-binding method

Point defect formation within graphene and ultra-thin graphite films is considered by means of molecular dynamics simulations in the framework of the density functional tight-binding approach. The barrier energy for vacancy formation is estimated and two types of defect formation are revealed.     

Molecular dynamics simulations of structure and dynamics of organic molecular crystals

A set of model compounds covering a range of polarity and flexibility have been simulated using GAFF, CHARMM22, OPLS and MM3 force fields to examine how well classical molecular dynamics simulations can reproduce structural and dynamic aspects of organic molecular crystals. Molecular structure, crystal structure and thermal motion, including molecular reorientations and internal rotations, found from the simulations have been compared between force fields and with experimental data. The MM3 force field does not perform well in condensed phase simulations, while GAFF, CHARMM and OPLS perform very similarly. Generally molecular and crystal structure are reproduced well, with a few exceptions. The atomic displacement parameters (ADPs) are mostly underestimated in the simulations with a relative error of up to 70%. Examples of molecular reorientation and internal rotation, observed in the simulations, include in-plane reorientations of benzene, methyl rotations in alanine, decane, isopropylcyclohexane, pyramidal inversion of nitrogen in amino group and rotation of the whole group around the C-N bond. Frequencies of such dynamic processes were calculated, as well as thermodynamic properties for reorientations in benzene and alanine. We conclude that MD simulations can be used for qualitative analysis, while quantitative results should be taken with caution. It is important to compare the outcomes from simulations with as many experimental quantities as available before using them to study or quantify crystal properties not available from experiment.     

Molecular structure and vibrational assignment of α‐chloro acetylacetone: A density functional theory study

The intramolecular hydrogen bond, molecular structure, and vibrational frequencies of α‐chloro acetylacetone have been investigated. Fourier transform infrared and Fourier transform Raman spectra of this compound and its deuterated analogue were recorded in the regions 400–4,000 cm^?1 and 50–4,000 cm^?1, respectively. Rigorous normal coordinate analysis has been performed at the B3LYP/6‐311++G** level of theory for purposes of comparison. The complete vibrational assignment for TFAA has been made on the basis of the calculated potential energy distribution. We also applied the atoms in molecules theory and natural bond orbital method for the analysis of the hydrogen bond in α‐Chloro acetylacetone and acetylacetone. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Molecular structure and vibrational assignment of melaminium phthalate by density functional theory (DFT) and ab initio Hartree-Fock (HF) calculations

The molecular geometry, the normal mode frequencies and corresponding vibrational assignment of melaminium phthalate (C₃H₇N₆⁺·C₈H₅O₄⁻) in the ground state were performed by HF and B3LYP levels of theory using the 6-31G(d) basis set. The optimized bond length numbers with bond angles are in good agreement with the X-ray data. The vibrational spectra of melaminium phthalate which is calculated by HF and B3LYP methods, reproduces vibrational wave numbers with an accuracy which allows reliable vibrational assignments. The title compound has been studied in the 4000–100 cm⁻¹ region where the theoretical evaluation and assignment of all observed bands were made.     

Projected gradient algorithms for Hartree-Fock and density matrix functional theory calculations

We present projected gradient algorithms designed for optimizing various functionals defined on the set of N-representable one-electron reduced density matrices. We show that projected gradient algorithms are efficient in minimizing the Hartree-Fock or the Muller-Buijse-Baerends functional. On the other hand, they converge very slowly when applied to the recently proposed BBk (k=1,2,3) functionals [O. Gritsenko et al., J. Chem. Phys. 122, 204102 (2005)]. This is due to the fact that the BBk functionals are not proper functionals of the density matrix.     

Density functional theory for efficient ab initio molecular dynamics simulations in solution

We present a density functional for first-principles molecular dynamics simulations that includes the electrostatic effects of a continuous dielectric medium. It allows for numerical simulations of molecules in solution in a model polar solvent. We propose a smooth dielectric model function to model solvation into water and demonstrate its good numerical properties for total energy calculations and constant energy molecular dynamics.     

Structure, spectroscopy, and reactivity properties of porphyrin pincers: a conceptual density functional theory and time-dependent density functional theory study

Porphyrin and pincer complexes are both important categories of compounds in biological and catalytic systems. The idea to combine them is computationally investigated in this work. By employment of density functional theory (DFT), conceptual DFT, and time-dependent DFT approaches, structure, spectroscopy, and reactivity properties of porphyrin pincers are systematically studied for a selection of divalent metal ions. We found that the porphyrin pincers are structurally and spectroscopically different from their precursors and are more reactive in electrophilic and nucleophilic reactions. A few quantitative linear/exponential relationships have been discovered between bonding interactions, charge distributions, and DFT chemical reactivity indices. These results are implicative in chemical modification of hemoproteins and understanding chemical reactivity in heme-containing and other biologically important complexes and cofactors.     

Resonance Raman spectroscopy and density functional theory calculation study of photodecay dynamics of tetra(4-carboxyphenyl) porphyrin

The 397.9 nm, 416.0 nm and 435.7 nm resonance Raman spectra were acquired for meso-tetrakis(4-carboxyphenyl)porphyrin (TCPP) in tetrahydrofuran solution, and the Raman effect of relaxation dynamics was analyzed according to Herzberg-Teller (vibronic coupling) contributions. Density functional calculations were done to help the elucidation of the Soret (B(x) and B(y)-band) electronic transitions and the corresponding photo relaxation dynamics of TCPP. The spectra indicate that the Franck-Condon region photo relaxation dynamics upon S(0) → S(4) electronic transition are predominantly along the totally symmetric C(m)-ph stretch and Porphin ring breath stretch, and simultaneously along the asymmetric ν(C(m)-Phenyl) + δ(N-H) and ν(C(α)-C(m)-C(α))(as) + def (pyr) vibrational relaxation processes. The excited state structural dynamics of TCPP determined from the resonance Raman spectra show that the internal conversion between the B(y) and B(x) electronic states occurs in tens of femtoseconds, and the electronic relaxation dynamics were firstly interpreted taking into account the time-dependent wave packet theory and Herzberg-Teller (vibronic coupling) contributions.     

The Gaussian and augmented-plane-wave density functional method for ab initio molecular dynamics simulations

A new algorithm for density-functional-theory-based ab initio molecular dynamics simulations is presented. The Kohn–Sham orbitals are expanded in Gaussian-type functions and an augmented-plane-wave-type approach is used to represent the electronic density. This extends previous work of ours where the density was expanded only in plane waves. We describe the total density in a smooth extended part which we represent in plane waves as in our previous work and parts localised close to the nuclei which are expanded in Gaussians. Using this representation of the charge we show how the localised and extended part can be treated separately, achieving a computational cost for the calculation of the Kohn–Sham matrix that scales with the system size N as O(NlogN). Furthermore, we are able to reduce drastically the size of the plane-wave basis. In addition, we introduce a multiple-cutoff method that improves considerably the performance of this approach. Finally, we demonstrate with a series of numerical examples the accuracy and efficiency of the new algorithm, both for electronic structure calculations and for ab initio molecular dynamics simulations. Received: 15 December 1998 /Accepted: 18 February 1999 /Published online: 14 July 1999     

Structure simulation of ultrathin dichloromethane layer on a solid substrate by density functional theory and molecular dynamics simulations

The method for prediction of structural properties of ultrathin liquid layers has been developed on the base of the atomistic molecular dynamics (AMD) and the density functional theory (DFT). A comparative analysis of ultrathin dichloromethane layer density profiles on three types of solid flat substrates showed that these approaches can be effectively used as mutually complementary procedures to describe the structural properties of nanometer scale surface layers. We used AMD calculations to predict the dichloromethane layer density profile on a solid substrate. However, it is difficult and computationally expensive to calculate structural and thermodynamic layers properties. At the same time, DFT can retain the microscopic details of macroscopic systems at the calculative cost significantly lower than that used in AMD. Therefore, in context of DFT, the substrate potential parameters are adjusted to reproduce AMD data. Thus, the obtained potential allows us to compute structural characteristics and, further, can be used to predict other physical properties of ultrathin films within the DFT framework. For instance, we calculated the coefficient of thermal expansion of dichloromethane in the case of three different substrates such as graphite, silicon oxide, and gold.     

Molecular dynamics simulation and density functional theory insight into the cocrystal explosive of hexaazaisowurtzitane/nitroguanidine

Theoretical methods involving molecular dynamics (MD) simulation and density functional theory were performed to investigate the different molecular ratios, mechanical Properties, structure, trigger bond, and intermolecular interaction of hexaazaisowurtzitane (CL‐20)/nitroguanidine (NQ) cocrystal explosive. Results of MD simulation show that CL‐20 and NQ packed in ratios of 1:1 present the larger binding energy and better mechanical properties than any other molecular ratios, which indicates 1:1 cocrystal can form the stable crystal structure. Shorter length and larger dissociation energy of trigger bond in composite structure than in isolated CL‐20 component suggests that the cocrystal may exhibit less sensitive than CL‐20. Analyses of atoms in molecules, reduced density gradient, and natural bond orbital confirm that intermolecular interactions are mainly derived from a series of weak hydrogen bond and strong vdW forces, involving of NH···O, CH···O, CH···N, O···N, and O···O. Additionally, composite structures of 2 and 3 bringing us more attractive performance will act as a key role in constructing of CL‐20/NQ cocrystal explosive. © 2015 Wiley Periodicals, Inc.