Polarizable continuum model with the fragment molecular orbital-based time-dependent density functional theory

The polarizable continuum model (PCM) for describing the solvent effect was combined with the fragment molecular orbital-based time-dependent density functional theory (TDDFT). Several levels of the many-body expansion were implemented, and the importance of the many-body contributions to the singlet-excited states was discussed. To calibrate the accuracy, we performed a number of the model calculations using our method and the regular TDDFT in solution, applying them to phenol and polypeptides at the long-range corrected BLYP/6-31G* level. It was found that for systems up to 192 atoms the largest error in the excitation energy was 0.006 eV (vs. the regular TDDFT/PCM of the full system). The solvent shifts and the conformer effects were discussed, and the scaling was found to be nearly linear. Finally, we applied our method to the lowest singlet excitation of the photoactive yellow protein (PYP) in aqueous solution and determined the excitation energy to be in reasonable agreement with experiment. The excitation energy analysis provided the contributions of individual residues, and the main factors as well as their solvent shifts were determined.     

Reaction dynamics following electron capture of chlorofluorocarbon adsorbed on water cluster: a direct density functional theory molecular dynamics study

The electron capture dynamics of halocarbon and its water complex have been investigated by means of the full dimensional direct density functional theory molecular dynamics method in order to shed light on the mechanism of electron capture of a halocarbon adsorbed on the ice surface. The CF(2)Cl(2) molecule and a cyclic water trimer (H(2)O)(3) were used as halocarbon and water cluster, respectively. The dynamics calculation of CF(2)Cl(2) showed that both C-Cl bonds are largely elongated after the electron capture, while one of the Cl atoms is dissociated from CF(2)Cl(2) (-) as a Cl(-) ion. Almost all total available energy was transferred into the internal modes of the parent CF(2)Cl radical on the product state, while the relative translational energy of Cl(-) was significantly low due to the elongation of two C-Cl bonds. In the case of a halocarbon-water cluster system, the geometry optimization of neutral complex CF(2)Cl(2)(H(2)O)(3) showed that one of the Cl atoms interacts with n orbital of water molecules of trimer and the other Cl atom existed as a dangling Cl atom. After the electron capture, only one C-Cl bond (dangling Cl atom) was rapidly elongated, whereas the other C-Cl bond is silent during the reaction. The dangling Cl atom was directly dissociated from CF(2)Cl(2) (-)(H(2)O)(3) as Cl(-). The fast Cl(-) ion was generated from CF(2)Cl(2) (-)(H(2)O)(3) on the water cluster. The mechanism of the electron capture of halocarbon on water ice was discussed on the basis of the theoretical results.     

Electronic properties of liquid water by sequential molecular dynamics/density functional theory

Electronic properties of liquid water were analysed by a sequential molecular dynamics (MD)/density functional theory approach. MD simulations are based on a polarisable model for water. Emphasis was placed on the prediction of the water dipole moment, liquid state polarisability, ionisation potential (IP), and vertical electron affinity. The dipole moment of the water molecule in liquid water is not dependent on the number of molecules included in the quantum mechanical calculations. The polarisability of the water molecule in liquid water is 4% lower than its gas phase value. The IP of liquid water (9.7 ± 0.06 eV) is in good agreement with recent experimental data.     

Modeling ultrafast solvated electronic dynamics using time-dependent density functional theory and polarizable continuum model

A first-principles solvated electronic dynamics method is introduced. Solvent electronic degrees of freedom are coupled to the time-dependent electronic density of a solute molecule by means of the implicit reaction field method, and the entire electronic system is propagated in time. This real-time time-dependent approach, incorporating the polarizable continuum solvation model, is shown to be very effective in describing the dynamical solvation effect in the charge transfer process and yields a consistent absorption spectrum in comparison to the conventional linear response results in solution.     

Bubble nucleation in micellar solution: a density functional study

We explore the free energetics of bubble nucleation in the micellar solution subjected to a negative pressure using a density functional model of a non-ionic surfactant solution. In this two-component model, the solvent is represented by a single hard-core sphere and the surfactant is represented by two tangent hard-core spheres connected by a rigid bond. The attractive interactions between the particles are modeled by the simple 1/R(6) form. Under all conditions of pressure and interparticle interactions we studied, the free energy barrier of bubble nucleation is found to be lower in the binary surfactant solution than that in a pure solvent and to continue to decrease as the mole fraction of the surfactant in the solution increases. We analyze the free energy surface of the model system under the conditions where both the critical bubble nucleus and the stable micelle exist in equilibrium with the same metastable solution. Our study shows that at moderately low pressures, bubbles can nucleate from the stable micelle and that the resulting free energy barrier of bubble nucleation is expected to be lower than that in the absence of this mechanism. However, as the spinodal is approached at lower pressures, the mechanism of micelle-assisted bubble nucleation becomes less effective. The liquid-liquid miscibility of the model system correlates well with the mechanism of bubble nucleation from the stable micelle.     

Applicability of hybrid density functional theory methods to calculation of molecular hyperpolarizability

The donor/acceptor (D/A) substituted pi-conjugated organic molecules possess extremely fast nonlinear optical (NLO) response time that is purely electronic in origin. This makes them promising candidates for optoelectronic applications. In the present study, we utilized four hybrid density functionals (B3LYP, B97-2, PBE0, BMK), Hartree-Fock, and second order Moller-Plesset correlation energy correction, truncated at second-order (MP2) methods with different basis sets to estimate molecular first hyperpolarizability (beta) of D/A-substituted benzenes and stilbenes (D=OMe, OH, NMe(2), NH(2); A=NO(2), CN). The results of density functional theory (DFT) calculations are compared to those of MP2 method and to the experimental data. We addressed the following questions: (1) the accurate techniques to compare calculated results to each other and to experiment, (2) the choice of the basis set, (3) the effect of molecular planarity, and (4) the choice of the method. Comparison of the absolute values of hyperpolarizabilities obtained computationally and experimentally is complicated by the ambiguities in conventions and reference values used by different experimental groups. A much more tangible way is to compare the ratios of beta's for two (or more) given molecules of interest that were calculated at the same level of theory and measured at the same laboratory using the same conventions and reference values. Coincidentally, it is the relative hyperpolarizabilities rather than absolute ones that are of importance in the rational molecular design of effective NLO materials. This design includes prediction of the most promising candidates from particular homologous series, which are to be synthesized and used for further investigation. In order to accomplish this goal, semiquantitative level of accuracy is usually sufficient. Augmentation of the basis set with polarization and diffuse functions changes beta by 20%; however, further extension of the basis set does not have significant effect. Thus, we recommend 6-31+G(*) basis set. We also show that the use of planar geometry constraints for the molecules, which can somewhat deviate from planarity in the gas phase, leads to sufficient accuracy (with an error less than 10%) of predicted values. For all the molecules studied, MP2 values are in better agreement with experiment, while DFT hybrid methods overestimate beta values. BMK functional gives the best agreement with experiment, with systematic overestimation close to the factor of 1.4. We propose to use the scaled BMK results for prediction of molecular hyperpolarizability at semiquantitative level of accuracy.     

Dynamic behavior of chemical reactivity indices in density functional theory: A Bohn-Oppenheimer quantum molecular dynamics study

Dynamic behaviors of chemical concepts in density functional theory such as frontier orbitals (HOMO/LUMO), chemical potential, hardness, and electrophilicity index have been investigated in this work in the context of Bohn-Oppenheimer quantum molecular dynamics in association with molecular conformation changes. Exemplary molecular systems like CH ₅ ⁺ , Cl⁻ (H₂O)₃₀ and Ca²⁺ (H₂O)1₅ are studied at 300 K in the gas phase, demonstrating that HOMO is more dynamic than LUMO, chemical potential and hardness often fluctuate concurrently. It is argued that DFT concepts and indices may serve as a good framework to understand molecular conformation changes as well as other dynamic phenomena.     

The structure of 1-chlorosilatrane: An ab initio molecular orbital and a density functional theory study

The molecular geometries of the 1-chloro-, 1-fluoro-, 1-methyl-, and 1-hydrogenosilatranes were fully optimized by the restricted Hartree-Fock (HF) method supplemented with 3-21G, 3-21G(d), 6-31G(d), and CEP-31G(d) basis sets; by MP2 calculations using 6-31G(d) and CEP-31G(d) basis sets; and by GGA-DFT calculations using 6-31G(d5) basis set with the aim of locating the positions of the local minima on the energy hypersurface. The HF/6-31G(d) calculations predict long (>254 pm) and the MP2/CEP calculations predicted short (∼225 pm) equilibrium Si(SINGLE BOND)N distances. The present GGA-DFT calculations reproduce the available gas phase experimental Si(SINGLE BOND)N distances correctly. The solid phase experimental results predict that the Si(SINGLE BOND)N distance is shorter in 1-chlorosilatrane than in 1-fluorosilatrane. In this respect the HF results show a strong basis set dependence, the MP2/CEP results contradict the experiment, and the GGA-DFT results in electrolytic medium agree with the experiment. The latter calculations predict that 1-chlorosilatrane is more polarizable than 1-fluorosilatrane and also support a general Si(SINGLE BOND)N distance shortening trend for silatranes during the transition from gas phase to polar liquid or solid phase. The calculations predict that the ethoxy links of the silatrane skeleton are flexible. Consequently, it is difficult to measure experimentally the related bond lengths and bond and torsion angles. This is the probable origin of the surprisingly large differences for the experimental structural parameters. On the basis of experimental analogies, ab initio calculations, and density functional theory (DFT) calculations, a gas phase equilibrium (r_e) geometry is predicted for 1-chlorosilatrane. The semiempirical methods predict a so-called exo minimum (at above 310 pm Si(SINGLE BOND)N distance); however, the ab initio and GGA-DFT calculations suggest that this form is nonexistent. The GGA-DFT geometry optima were characterized by frequency analysis. © 1996 by John Wiley & Sons, Inc.     

Nuclear magnetic shielding of the 113Cd(II) ion in aqua solution: a combined molecular dynamics/density functional theory study

We present a combined molecular dynamics simulation and density functional theory investigation of the nuclear magnetic shielding constant of the (113)Cd(II) ion solvated in aqueous solution. Molecular dynamics simulations are carried out for the cadmium-water system in order to produce instantaneous geometries for subsequent determination of the nuclear magnetic shielding constant at the density functional theory level. The nuclear magnetic shielding constant is computed using a perturbation theory formalism, which includes nonrelativistic and leading order relativistic contributions to the nuclear magnetic shielding tensor. Although the NMR shielding constant varies significantly with respect to simulation time, the value averaged over increasing number of snapshots remains almost constant. The paramagnetic nonrelativistic contribution is found to be most sensitive to dynamical changes in the system and is mainly responsible for the thermal and solvent effects in solution. The relativistic correction features very little sensitivity to the chemical environment, and can be disregarded in theoretical calculations when a Cd complex is used as reference compound in (113)Cd NMR experiments, due to the mutual cancelation between individual relativistic corrections.     

Acidity constants from vertical energy gaps: density functional theory based molecular dynamics implementation

The question of how to compute acidity constants (pK(a)) treating solvent and solute at the same level of theory remains of some interest, for example in the case of high or low pH conditions. We have developed a density functional theory based molecular dynamics implementation of such a method. The method is based on a half reaction scheme computing free energies of dissociation from the vertical energy gaps for insertion or removal of protons. Finite system size effects are important, but largely cancel when half reactions are combined to full reactions. We verified the method by investigating a series of organic and inorganic acids and bases spanning a wide range of pK(a) values (20 units). We find that the response of the aqueous solvent to vertical protonation/deprotonation is almost always asymmetric and correlated with the strength of the hydrogen bonding of the deprotonated base. We interpret these observations in analogy with the picture of solvent response to electronic ionization.     

Molecular dynamics effects on luminescence properties of oligothiophene derivatives: a molecular mechanics-response theory study based on the CHARMM force field and density functional theory

CHARMM force field parameter values for a class of oligothiophene derivatives have been derived with reference to density functional theory/B3LYP potential energy surfaces. The force field parametrization of these luminescent conjugated polyelectrolytes includes the electronic ground state as well as the strongly light absorbing first excited state. In conjunction with quantum chemical response theory calculations of transition state properties, a molecular dynamical model of the Stokes shift is obtained. The theoretical model is benchmarked against experimental data recorded at room temperature which refer to sodium salts of p-HTAA and p-FTAA with distilled water as a solvent. For p-HTAA the theoretically predicted Stokes shift of 112 nm is in good agreement with the experimental result of 124 nm, given the approximations about exciton localization that were introduced to obtain a force field for the excited state.     

Adsorption of copolymers in a selective nanoslit: a hybrid density functional theory

A hybrid density functional theory (DFT) is developed for adsorption of copolymers in a selective nanoslit. The DFT incorporates a single-chain simulation for the ideal-gas free energy functional with two weighted density approximations for the residual free energy functional. The theory is found to be insensitive to the width parameter used in the weighted density. Theoretical predictions are in excellent agreement with simulation results in the segment density profiles and the adsorption configurations including tail, loop, and train for copolymers with various sequences over a wide range of surface affinity. The bridge conformation is also observed in multiblock copolymers. Ordered assembly is facilitated in copolymers with longer chain/block and at stronger attraction between segment B and the slit wall. While diblock copolymer shows the longest tail, alternating copolymer has the shortest. As the attraction between segment B and the slit wall increases, the average size and fraction decrease for tail, but increase for loop and train.     

Lagrangian approach to molecular vibrational Raman intensities using time-dependent hybrid density functional theory

The authors propose a new route to vibrational Raman intensities based on analytical derivatives of a fully variational polarizability Lagrangian. The Lagrangian is constructed to recover the negative frequency-dependent polarizability of time-dependent Hartree-Fock or adiabatic (hybrid) density functional theory at its stationary point. By virtue of the variational principle, first-order polarizability derivatives can be computed without using derivative molecular orbital coefficients. As a result, the intensities of all Raman-active modes within the double harmonic approximation are obtained at approximately the same cost as the frequency-dependent polarizability itself. This corresponds to a reduction of the scaling of computational expense by one power of the system size compared to a force constant calculation and to previous implementations. Since the Raman intensity calculation is independent of the harmonic force constant calculation more, computationally demanding density functionals or basis sets may be used to compute the polarizability gradient without much affecting the total time required to compute a Raman spectrum. As illustrated for fullerene C60, the present approach considerably extends the domain of molecular vibrational Raman calculations at the (hybrid) density functional level. The accuracy of absolute and relative Raman intensities of benzene obtained using the PBE0 hybrid functional is assessed by comparison with experiment.     

An ab initio molecular dynamics and density functional theory study of the formation of phosphate chains from metathiophosphates

The reactions that occur between metathiophosphate (MTP) molecules are identified and examined through ab initio molecular dynamics simulations and static quantum chemical calculations at the density functional level of theory. The simulations show that certain types of MTPs can react to yield phosphate chains, while others only dimerize. These differences are rationalized in terms of reaction energies and the electronic structures of these molecules. In the reaction leading to the formation of phosphate chains, the reactive center, a tri-coordinate phosphorus atom, is continually regenerated. A polymerization mechanism linking MTPs to phosphate chains is developed on the basis of these results. This information sheds light on the underlying processes that may be responsible for the formation of phosphates under high-temperature conditions and may prove useful in the development of protocols for the rational synthesis of complex phosphate structures.     

Exploring odd-even effects of simple oligomer-like DRCNnT series: a study based on density functional theory/time-dependent density functional theory calculations

Odd-even effects of short-circuit current density and power conversion efficiency (PCE) are an interesting phenomenon in some organic solar cells. Although some explanations have been given, why they behave in such a way is still an open question. In the present work, we investigate a set of acceptor-donor-acceptor simple oligomer-like small molecules, named the DRCNnT (n = 5-9) series, to give an insight into this phenomenon because the solar cells based on them have high PCE (up to 10.08%) and show strong odd-even effects in experiments. By modeling the DRCNnT series and using density functional theory, we have studied the ground-state electronic structures of the DRCNnT (n = 5-9) series in condensed phase. The calculated results reproduce the experimental trends well. Furthermore, we find that the exciton-binding energies of the DRCNnT series may be one of the key parameters to explain this phenomenon because they also show odd-even effects. In addition, by studying the effects of alkyl branch and terminal group on odd-even effects of dipole moment, we find that eliminating one or two alkyl branches does not break the odd-even effects of dipole moments, but eliminating one or two terminal groups does. Finally, we conclude that removing one alkyl branch close to the terminal group of DRCN5T can decrease highest occupied molecular orbital (HOMO) energy (thus increasing open circuit voltage) and increase dipole moment (thus enhancing charge separation and short-circuit current). This could be a new and simple method to increase the PCE of DRCN5T-based solar cells.     

Homogeneous nucleation in vapor-liquid phase transition of Lennard-Jones fluids: a density functional theory approach

Density functional theory (DFT) with square gradient approximation for the free energy functional and a model density profile are used to obtain an analytical expression for the size-dependent free energy of formation of a liquid drop from the vapor through the process of homogeneous nucleation, without invoking the approximations used in classical nucleation theory (CNT). The density of the liquid drop in this work is not the same as the bulk liquid density but it corresponds to minimum free energy of formation of the liquid drop. The theory is applied to study the nucleation phenomena from supersaturated vapor of Lennard-Jones fluid. The barrier height predicted by this theory is significantly lower than the same in CNT which is rather high. The density at the center of the small liquid drop as obtained through optimization is less than the bulk density which is in agreement with other earlier works. Also proposed is a sharp interface limit of the proposed DFT of nucleation, which is as simple as CNT but with a modified barrier height and this modified classical nucleation theory, as we call it, is shown to lead to improved results.     

A computational study of microsolvation effect on ethylene glycol by density functional method

This study focuses on the conformational analysis of ethylene glycol-(water)n (n=1-3) complex by using density functional theory method and the basis set 6-311++G*. Different conformers are reported and the basis set superposition error corrected total energy is -306.767 5171, -383.221 3135, and -459.694 1528 for lowest energy conformer with 1, 2, and 3 water molecules, respectively, with corresponding binding energy -7.75, -15.43, and -36.28 kcal/mol. On applying many-body analysis it has been found that relaxation energy, two-body, three-body energy have significant contribution to the binding energy for ethylene glycol-(water)3 complex whereas four-body energies are negligible. The most stable conformers of ethylene glycol-(water)n complex are the cyclic structures in which water molecules bridge between the two hydroxyl group of ethylene glycol.     

Lattice density functional theory of molecular diffusion

A density functional theory of diffusion is developed for lattice fluids with molecular flux as a functional of the density distribution. The formalism coincides exactly with the generalized Ono-Kondo density functional theory when there is no gradient of chemical potential, i.e., at equilibrium. Away from equilibrium, it gives Fick's first law in the absence of a potential energy gradient, and it departs from Fickian behavior consistently with the Maxwell-Stefan formulation. The theory is applied to model a nanopore, predicting nonequilibrium phase transitions and the role of surface diffusion in the transport of capillary condensate.     

Mechanisms of difluoroethylene ozonolysis: a density functional theory study

We present a density functional theory (DFT) study on the mechanisms of gas-phase ozonolysis of three isomers of difluoroethylene, namely, cis-1,2-difluoroethylene, trans-1,2-difluoroethylene, and 1,1-difluoroethylene. MPW1K/cc-pVDZ and BHandHLYP/cc-pVDZ methods are employed to optimize the geometries of stationary points as well as the points on the minimum energy path (MEP). The energies of all the points were further refined at the QCISD(T)/cc-pVDZ and QCISD(T)/6-31+G(df,p) levels of theory with zero-point energy (ZPE) corrections. The ozone-cis-1,2-difluoroethylene reaction is predicted to be slower than the ozone-trans-1,2-difluoroethylene reaction. The enhanced reactivity of trans-1,2-difluoroethylene relative to the cis isomer is similar to the reactions of ozone with cis- and trans-dichloroethylene. The ozone-1,1-difluoroethylene reaction is predicted to be slower than the ozone-trans-1,2-difluoroethylene reaction. These results are in agreement with experimental studies. The calculated mechanisms indicate that in ozone-difluoroethylene reactions the yields of OH might be trivial, which is different from the reactions of ozone with unsaturated hydrocarbons.