Relativistic and electron‐correlation effects on magnetizabilities investigated by the Douglas‐Kroll‐Hess method and the second‐order Møller‐Plesset perturbation theory

Isotropic and anisotropic magnetizabilities for noble gas atoms and a series of singlet and triplet molecules were calculated using the second‐order Douglas‐Kroll‐Hess (DKH2) Hamiltonian containing the vector potential A and in part using second‐order generalized unrestricted Møller‐Plesset (GUMP2) theory. The DKH2 Hamiltonian was resolved into three parts (spin‐free terms, spin‐dependent terms, and magnetic perturbation terms), and the magnetizabilities were decomposed into diamagnetic and paramagnetic terms to investigate the relativistic and electron‐correlation effects in detail. For Ne, Kr, and Xe, the calculated magnetizabilities approached the experimental values, once relativistic and electron‐correlation effects were included. For the IF molecule, the magnetizability was strongly affected by the spin‐orbit interaction, and the total relativistic contribution amounted to 22%. For group 17, 16, 15, and 14 hydrides, the calculated relativistic effects were small (less than 3%), and trends were observed in relativistic and electron‐correlation effects across groups and periods. The magnetizability anisotropies of triplet molecules were generally larger than those of similar singlet molecules. The so‐called relativistic‐correlation interference for the magnetizabilities computed using the relativistic GUMP2 method can be neglected for the molecules evaluated, with exception of triplet SbH. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Description of core excitations by time-dependent density functional theory with local density approximation, generalized gradient approximation, meta-generalized gradient approximation, and hybrid functionals

Time-dependent density functional theory (TDDFT) is employed to investigate exchange-correlation-functional dependence of the vertical core-excitation energies of several molecules including H, C, N, O, and F atoms. For the local density approximation (LDA), generalized gradient approximation (GGA), and meta-GGA, the calculated X1s-->pi* excitation energies (X = C, N, O, and F) are severely underestimated by more than 13 eV. On the other hand, time-dependent Hartree-Fock (TDHF) overestimates the excitation energies by more than 6 eV. The hybrid functionals perform better than pure TDDFT because HF exchange remedies the underestimation of pure TDDFT. Among these hybrid functionals, the Becke-Half-and-Half-Lee-Yang-Parr (BHHLYP) functional including 50% HF exchange provides the smallest error for core excitations. We have also discovered the systematic trend that the deviations of TDHF and TDDFT with the LDA, GGA, and meta-GGA functionals show a strong atom-dependence. Namely, their deviations become larger for heavier atoms, while the hybrid functionals are significantly less atom-dependent.     

Time-dependent approach to electronically excited states of molecules with the multiconfiguration time-dependent Hartree-Fock method

In this paper the authors show how the multiconfiguration time-dependent Hartree-Fock (MCTDHF) method can be used for the calculation of electronic properties of molecules associated with the population of excited states. In contrast to other methods for correlated electron dynamics, such as configuration interaction, MCTDHF does not rely on a solution of the electronic Schrodinger equation prior to the propagation. The authors apply this approach to the calculation of vertical excitation energies, transition dipole moments, and oscillator strengths for two test molecules, lithium hydride and methane.     

van der Waals interactions are critical in Car–Parrinello molecular dynamics simulations of porphyrin–fullerene dyads

The interplay between electrostatic and van der Waals (vdW) interactions in porphyrin‐C₆₀ dyads is still under debate despite its importance in influencing the structural characteristics of such complexes considered for various applications in molecular photovoltaics. In this article, we sample the conformational space of a porphyrin‐C₆₀ dyad using Car–Parrinello molecular dynamics simulations with and without empirical vdW corrections. Long‐range vdW interactions, which are poorly described by the commonly used density functional theory functionals, prove to be essential for a proper dynamics of the dyad moieties. Inclusion of vdW corrections brings porphyrin and C₆₀ close together in an orientation that is in agreement with experimental observations. The structural differences arising from the vdW corrections are shown to be significant for several properties and potentially less important for others. Additionally, our Mulliken population analysis reveals that contrary to the common belief, porphyrin is not the primary electron donating moiety for C₆₀. In the considered dyad, fullerene's affinity for electrons is primarily satisfied by charge transfer from the amide group of the linker. However, we show that in the absence of another suitable bound donor, C₆₀ can withdraw electrons from porphyrin if it is sufficiently close. © 2015 Wiley Periodicals, Inc.     

Implementation of divide-and-conquer method including Hartree-Fock exchange interaction

The divide-and-conquer (DC) method, which is one of the linear-scaling methods avoiding explicit diagonalization of the Fock matrix, has been applied mainly to pure density functional theory (DFT) or semiempirical molecular orbital calculations so far. The present study applies the DC method to such calculations including the Hartree-Fock (HF) exchange terms as the HF and hybrid HF/DFT. Reliability of the DC-HF and DC-hybrid HF/DFT is found to be strongly dependent on the cut-off radius, which defines the localization region in the DC formalism. This dependence on the cut-off radius is assessed from various points of view: that is, total energy, energy components, local energies, and density of states. Additionally, to accelerate the self-consistent field convergence in DC calculations, a new convergence technique is proposed.     

Development of spin‐dependent relativistic open‐shell Hartree–Fock theory with time‐reversal symmetry (I): The unrestricted approach

An open‐shell Hartree–Fock (HF) theory for spin‐dependent, two‐component relativistic calculations, termed the Kramers‐unrestricted HF (KUHF) method, is developed. The present KUHF method, which is formulated as a relativistic counterpart of nonrelativistic UHF, is based on quaternion algebra and partly uses time‐reversal symmetry. The fundamental characteristics of KUHF are discussed in this study. From numerical assessments, it was revealed that KUHF gives a corresponding solution to nonrelativistic UHF; furthermore, KUHF properly describes spin‐orbit interactions. In addition, KUHF can improve the self‐consistent field convergence behavior in spin‐dependent calculations, for example, for f‐block elements.     

New analytic approximation to the standard molecular volume definition and its application to generalized Born calculations

In a recent article (Lee, M. S.; Salsbury, F. R. Jr.; Brooks, C. L., III. J Chem Phys 2002, 116, 10606), we demonstrated that generalized Born (GB) theory provides a good approximation to Poisson electrostatic solvation energy calculations if one uses the same definitions of molecular volume for each. In this work, we present a new and improved analytic method for reproducing the Lee-Richards molecular volume, which is the most common volume definition for Poisson calculations. Overall, 1% errors are achieved for absolute solvation energies of a large set of proteins and relative solvation energies of protein conformations. We also introduce an accurate SASA approximation that uses the same machinery employed by our GB method and requires a small addition of computational cost. The combined methodology is shown to yield an efficient and accurate implicit solvent representation for simulations of biopolymers.     

Development of minimized mixing molecular orbital method for designing organic ferromagnets

Predicting the high spin stability of organic radicals correctly for designing organic ferromagnets remains a significant challenge. We have developed a method with an index (L_min) for predicting the high spin stability of conjugated organic radicals at the restricted open‐shell Hartree–Fock level. Unitary transformations were performed for localizing the coefficients of nonbonding molecular orbitals, and subsequently the localized coefficients were used to calculate L_min that indicates the high spin stability of conjugated organic radicals. This method can be combined with the elongation method to treat huge high spin open‐shell systems. Thus, this method is useful for designing organic ferromagnets. © 2015 Wiley Periodicals, Inc.     

A non-orthogonal Kohn-Sham method using partially fixed molecular orbitals

A density functional theory method using partially fixed molecular orbitals (PFMOs) is presented. The PFMOs, which have some fixed molecular orbital coefficients and are non-orthogonal, are a generalization of the extreme localized orbitals (ELMOs) of Couty, Bayse, and Hall (1997) Theor Chem Acc 97:96. A non-orthogonal Kohn-Sham method with these PFMOs is derived, and is applied to molecular calculations on the ionization potential of pyridine, the energy difference between cis- and trans-butadiene, the reaction barrier height of the cyclobutene-cis-butadiene interconversion, and the potential energy curve of the hydrogen shift reaction of hydroxycarbene to formaldehyde. The PFMO Kohn-Sham method reproduces well the results of the full Kohn-Sham method without having a restriction on the molecular orbital coefficients. The difference is less than 0.1 eV in the ionization potential and about 0.1 kcal/mol in the barrier height and in the potential energy calculations.     

Mixing of generalized parity in molecular orbitals

The concept of generalized parity is introduced. It allows mixing of different symmetries in molecular orbitals in the framework of the Parity Mixing in Orbitals method. An extension of this SCF calculational scheme is also discussed and the relevant secular equations are reported.Alexander von Humboldt Fellow. On leave from the Institute Rudjer Bokovi, Zagreb, Croatia, Yugoslavia.     

Balancing an accurate representation of the molecular surface in generalized born formalisms with integrator stability in molecular dynamics simulations

Different integrator time steps in NVT and NVE simulations of protein and nucleic acid systems are tested with the GBMV (Generalized Born using Molecular Volume) and GBSW (Generalized Born with simple SWitching) methods. The simulation stability and energy conservation is investigated in relation to the agreement with the Poisson theory. It is found that very close agreement between generalized Born methods and the Poisson theory based on the commonly used sharp molecular surface definition results in energy drift and simulation artifacts in molecular dynamics simulation protocols with standard 2-fs time steps. New parameters are proposed for the GBMV method, which maintains very good agreement with the Poisson theory while providing energy conservation and stable simulations at time steps of 1 to 1.5 fs.     

Exponential repulsion improves structural predictability of molecular docking

Molecular docking is a powerful tool for theoretical prediction of the preferred conformation and orientation of small molecules within protein active sites. The obtained poses can be used for estimation of binding energies, which indicate the inhibition effect of designed inhibitors, and therefore might be used for in silico drug design. However, the evaluation of ligand binding affinity critically depends on successful prediction of the native binding mode. Contemporary docking methods are often based on scoring functions derived from molecular mechanical potentials. In such potentials, nonbonded interactions are typically represented by electrostatic interactions between atom‐centered partial charges and standard 6–12 Lennard–Jones potential. Here, we present implementation and testing of a scoring function based on more physically justified exponential repulsion instead of the standard Lennard–Jones potential. We found that this scoring function significantly improved prediction of the native binding modes in proteins bearing narrow active sites such as serine proteases and kinases. © 2016 Wiley Periodicals, Inc.     

Comparison of localization and delocalization indices obtained with Hartree-Fock and conventional correlated methods: effect of Coulomb correlation

Atomic populations and localization [lambda(A)] and delocalization [delta(A,B)] indices (LIs and DIs) are calculated for a large set of molecules at the Hartree-Fock (HF), MP2, MP4(SDQ), CISD, and QCISD levels with the 6-311++G(2d,2p) basis set. The HF method and the conventional correlation methods [MP2, MP4(SDQ), CISD, and QCISD] yield distinct sets of LIs and DIs. Yet, within the four conventional correlation methods the differences in atomic populations and LIs and DIs are small. Relative to HF, the conventional correlation methods [MP2, MP4(SDQ), CISD, QCISD] yield virtually the same LIs and DIs for molecules with large charge separations while LIs and DIs that differ significantly from the HF values--the LIs are increased and DIs decreased--are obtained for bonds with no or small charge separations. Such is the case in the archetypal homopolar molecules HC(triple bond)CH, H2C=CH2, CH3-CH3, and "protonated cyclopropane" C(3)H(7) (+), in which case the bonding may be atypical. Relative to HF, the typical effect of the conventional correlation methods is to decrease the DI between atoms.     

Molecular view of charge exchange in ion–C60 collisions

We present a theoretical study of charge transfer in H⁺+C₆₀ and He²⁺+C₆₀ collisions using an extension of the molecular time‐dependent method of ion–atom collisions. Energy‐correlation diagrams have been evaluated for the corresponding (C₆₀–H)⁺ and (C₆₀–He)²⁺ quasi‐molecules. Single and double charge‐transfer cross sections in C₆₀+He²⁺ collisions are reported for the first time. The results show that double charge‐transfer cross sections are only one order of magnitude smaller than single charge‐transfer cross sections. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Ab initio molecular dynamics of protonated water clusters by integrated multicenter molecular-orbital method

To analyze large-scale cluster systems theoretically, we recently developed an "integrated multicenter molecular-orbital" (IMiC-MO) method. This method calculates the force of an entire system by dividing the system into small regions. We used the method to analyze the effect of cluster size and the process of hydrogen bond network (HBN) growth to form H(+)(H(2)O)(n) (n = 9, 17, and 33) clusters. Our simulations reveal that H(3)O(+) and water molecules in the first solvation shell function take an important role to grow the HBN. In addition, the number of hydrogen donors in each water molecule is strongly related to the shape of the HBN.     

GC–MS with photoionization of cold molecules in supersonic molecular beams—Approaching the softest ionization method

A new type of photoionization ion source was developed for the ionization of cold molecules in supersonic molecular beams (named Cold PI). The system was based on a GC–MS with supersonic molecular beams and its fly‐through EI of cold molecules ion source (Cold EI) plus quadrupole mass analyzer. A continuously operated deuterium VUV photoionization lamp was added and placed above and between the supersonic nozzle and skimmer whereas the Cold EI ion source served only as a portion of the ion transfer ion optics. The supersonic nozzle and skimmer were voltage biased and the VUV light crossed the supersonic expansion about 10 mm from the nozzle. We obtained over three orders of magnitude enhancement in the relative abundance of the molecular ion of squalane in Cold PI versus in photoionization of this compound as a thermal compound. Accordingly, we also proved that standard photoionization is not as soft ionization method as previously perceived for large compounds. We found that Cold PI is as soft as and possibly softer than field ionization; thus, it could be the softest known ionization method. The ionization yield was about 200–300 times weaker than with Cold EI yet our limit of detection was about 200 femtogram in SIM mode for cholesterol and pyrene which is reasonable. Practically, all hydrocarbons gave only molecular ions with rather uniform response whereas alcohols gave some molecular ions plus major fragment ions particularly with a loss of water (similarly to field ionization). We tested Cold PI in the GC–MS analysis of diesel fuels and analyzed the time averaged data for group type information. We also found that we can analyze the diesel fuels by fast under 20‐s flow injection analysis in which the generated averaged mass spectrum of molecular ions only could serve for the characterization of fuels.     

Molecular dynamics method in studies of molecular film growth processes

This review is an attempt to analyze some of the experimental problems arising in the course of growth of orientated molecular films using metal monophthalocyanine (MPc) films as an example and to demonstrate the possibilities of molecular dynamics simulation of these processes for solving experimental problems. Examples of theoretical simulation of adsorption processes are given; formation of a molecular monolayer is considered for copper phthalocyanine films as an example.