Ground State Minimum Research Articles

Towards experimental determination of conical intersection properties: a twin state based comparison with bound excited states

The energy and approximate structure of certain S(0)/S(1) conical intersections (CI) are shown computationally to be deducible from those of two bound states: the first triplet (T(1)), which is iso-energetic with the CI, and the second excited singlet state (S(2)). This is demonstrated for acepentalene (I) and its perfluoro derivative (II) using the twin state concept for three states systems and based on the fact that the triplet T(1) is almost degenerate with the CI. The stable S(2) (C(3v) configuration) state exhibits unusual exaltation of Jahn-Teller active degenerate mode-ν(JT) = 2058 cm(-1) (∼500 cm(-1) higher than analogous e-mode of the symmetric (C(3v)) T(1) and the dianion I(-2) or any C-C vibration of the Jahn-Teller distorted (C(s)) ground state minimum). The acepentalene molecule, whose rigid structure and possibility to attain the relatively high symmetry C(3v) configuration, is a particularly suitable candidate for this purpose.

Calcium-containing diatomic dications in the gas phase

Sputtering (ion surface bombardment) of various calcium-containing powder samples with an energetic (17 keV), high-current (16)O(-) beam has produced the diatomic dications of CaSi(2+), CaP(2+), CaF(2+), CaH(2+), CaCl(2+), CaBr(2+) and CaI(2+). These molecular gas-phase species have been identified in positive ion mass spectra at half-integer m/z values; their ion flight times through a magnetic-sector mass spectrometer were roughly 10(-5) s. Most of them appear to be novel molecular ions; the stability of the latter four (CaH(2+), CaCl(2+), CaBr(2+) and CaI(2+)) had been demonstrated in previous theoretical studies, whereas only CaF(2+) and CaBr(2+) had been observed before. Here we combine the results of our experimental search with a detailed theoretical study of the remaining three systems CaSi(2+), CaP(2+) and CaF(2+). All electronic states correlating with the first dissociation channel are characterized using high level ab initio electronic structure calculations. In their ground states, we find CaSi(2+) to be a long-lived metastable molecule, whereas CaF(2+) and CaP(2+) are thermodynamically stable, with respective equilibrium internuclear distances of 6.253, 4.740, and 5.731 a(0). CaSi(2+) has a well depth of 7116 (0.88) cm(-1) (eV) and a dissociation asymptote 7956 (0.99) cm(-1) (eV) below the ground state minimum. The dissociation energy of CaF(2+) is estimated to be 3404 (0.42) cm(-1) (eV), whereas for CaP(2+) we found 2547 (0.32) cm(-1) (eV), and a barrier height of 8118 (1.01) cm(-1) (eV). Their adiabatic double ionisation energies are 22.87, 16.91, and 17.32 eV, respectively, for the F, Si, and P containing dications.

Molecular aniline clusters. II. The low-lying electronic excited states

The lowest electronically excited states of the aniline dimer and trimer related to the lowest π(∗)←π transition of the monomer are investigated by applying time-dependent coupled cluster theory, primarily at the level of the (spin-component-scaled) CC2 model. Minimum energy structures in the vicinity of the Franck-Condon points were determined on the individual potential energy surfaces. For the dimer we find an excimer and a head-to-tail configuration (with the monomers substantially displaced relative to the ground state minimum) for the lowest (dark) and second lowest (bright) states, respectively. The excitation is delocalized on both chromophores for both of these states. For the trimer three distinct minima with quite different hydrogen-bonding arrangements are found for the three lowest states. In strong contrast to the dimer the excitation here is clearly localized on the individual aniline chromophores for each of these three states. One of the three geometries is rather similar to the ground state minimum, while the two others are rather different and thus have presumably quite small Franck-Condon factors. It can be expected that only the electronic origin of the first conformer can eventually be detected in the absorption spectrum of the trimer, provided that it is separated by high-enough barriers from other, energetically lower configurations.

TRIAXIAL SHAPE IN Os–Pt REGION FROM GROUND STATES TO COLLECTIVE ROTATIONAL STATES

In order to study the deformation of the nuclei which belongs to the transitional region between strongly deformed and spherical nuclei, the calculations of the total Routhian surfaces by means of the pairing-deformation-frequency self-consistent cranked shell model were carried out in order to investigate shape evolution in even-mass Os and Pt isotopes (176–202 Os and 182–204 Pt ) starting from the static nuclear ground states. It is found that some nuclear ground states such as in 196 Os and 188–194 Pt are neither oblate nor prolate. Instead, the ground-state minima in these nuclei are axially asymmetric in shape. In addition, a complementary approach is used to extract equilibrium γ0 value (nonaxially symmetric equilibrium shape parameter), which support our predictions.

Direct restricted-step MCSCF calculations on the structure and spectrum of cyclobutadiene

The structure of cyclobutadiene in the ground 1A9 state and the lowest excited 3B1g, 1A1g, and 1B1g, states are investigated by means of a direct, step-restricted, second order MCSCF method. The barrier height for interconversion between the rectangular forms of the ground state minima is predicted to be 6.2kcal/mol and the singlet-triplet state separation at the ground state square confirmation is predicted to be 9.7 kcal/mol. All investigated excited states are fond to have square equilibrium geometry the and adiabatic transition energy to the 1A1g, state is predicted to be lower than that to the 1B1g, state.

A Computational Study on the Mechanism of Intramolecular Oxo−Hydroxy Phototautomerism Driven by Repulsive πσ* State

Potential energy (PE) surfaces of the lowest excited states of the 4(3H)-pyrimidinone/4-hydroxypyrimidine system were investigated with the aid of the CC2 and CASSCF methods of the electronic structure theory. These studies resulted in identification of a low-lying pi sigma* state, which is dissociative with respect to the stretching of the N-H or O-H bonds in the oxo and hydroxy structures of the compound, respectively. After initial excitation to the lowest local n pi* and/or pi pi* singlet states, the system can access the PE surface of the pi sigma* state by crossing a low barrier. It was computationally demonstrated that the system should evolve on the PE surface of the repulsive pi sigma* state toward a broad seam of intersection with the PE surface of the ground state. At the intersection, the nonadiabatic transition to the ground electronic state takes place and the system can either evolve to a minimum of the initially excited tautomer or to the ground-state minimum of the other tautomer. The steps listed above provide a mechanism of photoinduced dissociation-association (PIDA) phototautomerism, experimentally observed for a number of monomeric molecules, structurally similar to 4(3H)-pyrimidinone/4-hydroxypyrimidine. This mechanism describes a new class of intramolecular phototautomeric reactions driven by a repulsive pi sigma* state.

Two- and three-state conical intersections in the uracil cation

The excited states of the uracil radical cation have been studied theoretically using multireference configuration interaction and quasidegenerate perturbation theory methods. The ground state is an 2A″ state and three electronically excited states have been calculated with adiabatic energies 1.00 eV (D 1( 2A′)), 1.79 eV (D 2( 2A″)), and 2.31 eV (D 3( 2A′)). Two-state conical intersections between the ground state D 0 and the first excited state D 1 have been located 0.45 eV above the ground state minimum. Furthermore, three-state conical intersections between the states D 0, D 1 and D 2 have been located 0.95 eV above the D 0 minimum. These conical intersections may play an important role in photoionization processes in uracil.

A ground-state-directed optimization scheme for the Kohn–Sham energy

Kohn-Sham density-functional calculations are used in many branches of science to obtain information about the electronic structure of molecular systems and materials. Unfortunately, the traditional method for optimizing the Kohn-Sham energy suffers from fundamental problems that may lead to divergence or, even worse, convergence to an energy saddle point rather than to the ground-state minimum--in particular, for the larger and more complicated electronic systems that are often studied by Kohn-Sham theory nowadays. We here present a novel method for Kohn-Sham energy minimization that does not suffer from the flaws of the conventional approach, combining reliability and efficiency with linear complexity. In particular, the proposed method converges by design to a minimum, avoiding the sometimes spurious solutions of the traditional method and bypassing the need to examine the structure of the provided solution.

Two-Dimensional H2O−Cl2 and H2O−Br2 Potential Surfaces: An Ab Initio Study of Ground and Valence Excited Electronic States

All electron ab initio calculations for the interaction of H2O with Cl2 and Br2 are reported for the ground state and the lowest triplet and singlet Pi excited states as a function of both the X-X and O-X bond lengths (X = Cl or Br). For the ground state and lowest triplet state, the calculations are performed with the coupled cluster singles, doubles, and perturbative triple excitation level of correlation using an augmented triple-zeta basis set. For the 1Pi state the multireference average quadratic coupled cluster technique was employed. For several points on the potential, the calculations were repeated with the augmented quadruple-zeta basis set. The ground-state well depths were found to be 917 and 1,183 cm-1 for Cl2 and Br2, respectively, with the triple-zeta basis set, and they increased to 982 and 1,273 cm-1 for the quadruple-zeta basis set. At the geometry of the ground-state minimum, the lowest energy state corresponding to the unperturbed 1Pi states of the halogens increases in energy by 637 and 733 cm-1, respectively, relative to the ground-state dissociation limit of the H2O-X2 complex. Adding the attractive ground-state interaction energy to that of the repulsive excited state predicts a blue-shift, relative to that of the free halogen molecules, of approximately 1,600 cm-1 for H2O-Cl2 and approximately 2,000 cm-1 for H2O-Br2. These vertical blue-shifts for the dimers are greater than the shift of the band maximum upon solvation of either halogen in liquid water.

On the vibronic coupling approximation: A generally applicable approach for determining fully quadratic quasidiabatic coupled electronic state Hamiltonians

In this report we introduce an iterative procedure for constructing a quasidiabatic Hamiltonian representing N(state)-coupled electronic states in the vicinity of an arbitrary point in N(int)-dimensional nuclear coordinate space. The Hamiltonian, which is designed to compute vibronic spectra employing the multimode vibronic coupling approximation, includes all linear terms which are determined exactly using analytic gradient techniques. In addition, all [N(state)][N(int)] quadratic terms, where [n]=n(n+1)/2, are determined from energy gradient and derivative coupling information obtained from reliable multireference configuration interaction wave functions. The use of energy gradient and derivative coupling information enables the large number of second order parameters to be determined employing ab initio data computed at a limited number of points (N(int) being minimal) and assures a maximal degree of quasidiabaticity. Numerical examples are given in which quasidiabatic Hamiltonians centered around three points on the C(3)H(3)N(2) potential energy surface (the minimum energy point on the ground state surface and the minimum energy points on the two- and three-state seams of conical intersection) were computed and compared. A method to modify the conical intersection based Hamiltonians to better describe the region of the ground state minimum is introduced, yielding improved agreement with ab initio results, particularly in the case of the Hamiltonian defined at the two-state minimum energy crossing.

A complete active space self-consistent field study of the photochemistry of nitrosamine

Photodissociation mechanisms of nitrosamine (NH2NO) have been studied at the complete active space self-consistent field level of theory in conjunction with atomic-natural-orbital-type basis sets. In addition, the energies of all the critical points and the potential energy curves connecting them have been recomputed with the multiconfigurational second-order perturbation method. Ground state minimum of nitrosamine has a C1 nonplanar structure with the hydrogen atoms of the amino moiety out of the plane defined by the N-N-O bonds. Electronic transitions to the three lowest states are allowed by selection rules: (i) S0-->S3 (7.41 eV) has an oscillator strength of f=0.0006 and it is assigned as an (npO)0-->(piNO*)2 transition, (ii) S0-->S2 (5.86 eV) has an oscillator strength of f=0.14 and it is assigned as an npN-->piNO* transition, and (iii) S0-->S1 (2.98 eV) has an oscillator strength of f=0.002 and it is assigned as an npO-->piNO* transition. It is found that N-N bond cleavage is the most likely process in all the photochemical relevant states, namely, S1 (1 1A"), S2 (2 1A'), and T1 (1 3A"). While S1 and T1 yield exclusively homolytic dissociation: NH2NO-->NH2 (1 2B1)+NO(X 2Pi), on S2 the latter process constitutes the major path, but two additional minor channels are also available: adiabatic homolytic dissociation: NH2NO-->NH2 (1 2A1)+NO(X 2Pi), and adiabatic oxygen extrusion: NH2NO-->NH2N (1 3A1)+O(3P). The excited species NH2 (1 2A1) experiences a subsequent ultrafast decay to the ground state, the final products in all cases the fragments being in their lowest electronic state. We have not found a unimolecular mechanism connecting excited states with the ground state. In addition, homolytic dissociation in the ground state, tautomerizations to NHNOH and NHNHO, and intersystem crossings to T1 are considered. The most favorable process on this state is the isomerization to NHNOH.

Efficient Photochemical Merocyanine-to-Spiropyran Ring Closure Mechanism through an Extended Conical Intersection Seam. A Model CASSCF/CASPT2 Study

A mechanism of the thermal and photochemical bleaching of merocyanine to spiropyran is proposed on the basis of CASSCF/CASPT2 calculations on the 6-(2-propenyliden)cyclohexadienone model system. Our results suggest that this photochemical transformation takes place in two steps. First, the initially pumped 1(pi-pi) S2 undergoes radiationless decay to 1(n-pi) S1 via an extended S2/S1 conical intersection seam that runs approximately parallel to the trans-to-cis isomerization coordinate, a few kilocalories per mole higher in energy. Thus, S2 --> S1 internal conversion is possible at all values of the S2 trans-to-cis reaction coordinate. Second, on the S1 potential energy surface, there is a barrierless ring closure reaction path from the S1 cis minimum that leads to a peaked S1/S0 conical intersection where the deactivation to the ground state takes place. The inertia of the moving nuclei then drives the system toward the ground-state minimum of the 2H-chromene product. Thus, the extended seam topology of the S2/S1 conical intersection and the coordinate of the branching space of the S1/S0 conical intersection are essential to explain the efficiency and high speed of this reaction.

Photoinduced intramolecular charge transfer (ICT) reaction in trans-methyl p-(dimethylamino) cinnamate: A combined fluorescence measurement and quantum chemical calculations

The photophysical behaviour of trans-methyl p-(dimethylamino) cinnamate ( t-MDMAC) donor–acceptor system has been investigated by steady-state absorption and emission spectroscopy and quantum chemical calculations. The molecule t-MDMAC shows an emission from the locally excited state in non-polar solvents. In addition to weak local emission, a strong solvent dependent red shifted fluorescence in polar aprotic solvents is attributed to highly polar intramolecular charge transfer state. However, the formation of hydrogen-bonded clusters with polar protic solvents has been suggested from a linear correlation between the observed red shifted fluorescence band maxima with hydrogen bonding parameters ( α). Calculations by ab initio and density functional theory show that the lone pair electron at nitrogen center is out of plane of the benzene ring in the global minimum ground state structure. In the gas phase, a potential energy surface along the twist coordinate at the donor (–NMe 2) and acceptor (–CH = CHCOOMe) sites shows stabilization of S 1 state and destabilization S 2 and S 0 states. A similar potential energy calculation along the twist coordinate in acetonitrile solvent using non-equilibrium polarized continuum model also shows more stabilization of S 1 state relative to other states and supports solvent dependent red shifted emission properties. In all types of calculations it is found that the nitrogen lone pair is delocalized over the benzene ring in the global minimum ground state and is localized on the nitrogen centre at the 90° twisted configuration. The S 1 energy state stabilization along the twist coordinate at the donor site and localized nitrogen lone pair at the perpendicular configuration support well the observed dual fluorescence in terms of proposed twisted intramolecular charge transfer (TICT) model.

Application of equation-of-motion coupled-cluster methods to low-lying singlet and triplet electronic states of HBO and BOH

The equilibrium structures and physical properties of the X (1)sigma(+) linear electronic states, linear excited singlet and triplet electronic states of hydroboron monoxide (HBO) (A (1)sigma(-), B (1)delta, a (3)sigma(+), and b (3)delta) and boron hydroxide (BOH) (A (1)sigma(+), B (1)Pi, and b (3)Pi), and their bent counterparts (HBO a (3)A('), b (3)A("), A (1)A("), B (1)A(') and BOH X (1)A('), b (3)A('), c (3)A("), A (1)A('), B (1)A('), C (1)A(")) are investigated using excited electronic state ab initio equation-of-motion coupled-cluster (EOM-CC) methods. A new implementation of open-shell EOM-CC including iterative partial triple excitations (EOM-CC3) was tested. Coupled-cluster wave functions with single and double excitations (CCSD), single, double, and iterative partial triple excitations (CC3), and single, double, and full triple excitations (CCSDT) are employed with the correlation-consistent quadruple and quintuple zeta basis sets. The linear HBO X (1)sigma(+) state is predicted to lie 48.3 kcal mol(-1) (2.09 eV) lower in energy than the BOH X (1)sigma(+) linear stationary point at the CCSDT level of theory. The CCSDT BOH barrier to linearity is predicted to lie 3.7 kcal mol(-1) (0.16 eV). With a harmonic zero-point vibrational energy correction, the HBO X (1)sigma(+)-BOH X (1)A(') energy difference is 45.2 kcal mol(-1) (1.96 eV). The lowest triplet excited electronic state of HBO, a (3)A('), has a predicted excitation energy (T(e)) of 115 kcal mol(-1) (4.97 eV) from the HBO ground state minimum, while the lowest-bound BOH excited electronic state, b (3)A('), has a T(e) of 70.2 kcal mol(-1) (3.04 eV) with respect to BOH X (1)A('). The T(e) values predicted for the lowest singlet excited states are A (1)A(")<--X (1)sigma(+)=139 kcal mol(-1) (6.01 eV) for HBO and A (1)A(')<--X (1)A(')=102 kcal mol(-1) (4.42 eV) for BOH. Also for BOH, the triplet vertical transition energies are b (3)A(')<--X (1)A(')=71.4 kcal mol(-1) (3.10 eV) and c (3)A(")<--X (1)A(')=87.2 kcal mol(-1) (3.78 eV).

De novo prediction of the ground state structure of transition metal complexes using semiempirical and ab initio quantum mechanics. Coordination isomerism

The ground state coordination isomers for 30 different trigonal bipyramidal transition metal complexes have been predicted using different levels of quantum mechanics: semiempirical (PM3(tm)), ab initio (MP2//HF), pure (BPW91) and hybrid (B3PW91) density functional theory (DFT) methods. For species where these methods failed to reproduce crystallographic data, hybrid quantum mechanics/molecular mechanics (QM/MM) methods were used to study more exact experimental models. Literature deficiencies regarding ground state multiplicity of these species were supplemented by spin predictions using previously tested PM3(tm) methods. Geometry optimization calculations were performed for each possible coordination isomer. The predicted ground state minima provided by the different methods are compared to each other and with crystallographic data. Pure DFT functionals outperformed hybrid functionals and MP2//HF. The very rapid PM3(tm) parameterization method provided accurate predictions in comparison to other levels of theory. An integrated MM/PM3(tm)/DFT de novo scheme accurately reproduced crystallographic data for species where the individual methods failed.

Adiabatic time-dependent density functional methods for excited state properties

This work presents theory, implementation, and validation of excited state properties obtained from time-dependent density functional theory (TDDFT). Based on a fully variational expression for the excited state energy, a compact derivation of first order properties is given. We report an implementation of analytic excited state gradients and charge moments for local, gradient corrected, and hybrid functionals, as well as for the configuration interaction singles (CIS) and time-dependent Hartree–Fock (TDHF) methods. By exploiting analogies to ground state energy and gradient calculations, efficient techniques can be transferred to excited state methods. Benchmark results demonstrate that, for low-lying excited states, geometry optimizations are not substantially more expensive than for the ground state, independent of the molecular size. We assess the quality of calculated adiabatic excitation energies, structures, dipole moments, and vibrational frequencies by comparison with accurate experimental data for a variety of excited states and molecules. Similar trends are observed for adiabatic excitation energies as for vertical ones. TDDFT is more robust than CIS and TDHF, in particular, for geometries differing significantly from the ground state minimum. The TDDFT excited state structures, dipole moments, and vibrational frequencies are of a remarkably high quality, which is comparable to that obtained in ground state density functional calculations. Thus, yielding considerably more accurate results at similar computational cost, TDDFT rivals CIS as a standard method for calculating excited state properties in larger molecules.

Structural organization of a chain molecule with specific charge distribution: A molecular dynamics study

We report the results of molecular dynamics simulations of a charged bead-monomer chain molecule with charge distribution adopted from immunoglobulin-binding domain B1 of protein-g. The beads of the model are connected by invariable bonds and interact with each other via the Coulomb potential. To study the low-temperature conformational space of the designed model we use standard canonical, microcanonical and multicanonical molecular dynamics simulations. We find that at low temperature T = T c the chain undergoes a continuous freezing transition into one of many low-energy conformations. Below T c the molecule is a compact globule composed of an inner core, containing mostly charged monomers, and an outer corona, consisting of all the rest neutral units. All frozen conformations have almost equal potential energy but differ in structure. The potential energy surface of the model does not posses a pronounced ground-state minimum--an essential feature of protein-like heteropolymers.

C−H···Carboxylate Oxygen Hydrogen Bonding in Substrate Activation by Acyl-CoA Dehydrogenases: Synergy between the H-bonds

Several model systems for the α-carbonyl proton abstraction reaction were studied at the DFT level to examine the role of hydrogen bonding interactions in the enzyme catalytic activity of acyl-CoA dehydrogenase. Each system contains the thioester portion of acyl-CoA (CH3CH2(CO)SCH3) a weakly basic COO⊖ representing Glu376, H-donor molecules (methanol and methylformamide) with zero, one or two H-donors bound to the carbonyl oxygen of the thioester correspond to the Glu376 and FAD side-chains. Five transition structures (TS-2a, TS-2b, TS-2d, TS-2a*, and TS-2b*) for the α-proton abstraction step and their corresponding ground state minima have been fully optimized at the B3LYP/6-31+G(d,p) level of theory. Classical activation energy barriers of 5.8, 10.4, 7.7, 9.9, and 12.3 kcal/mol have been calculated for these model systems, respectively. ChelpG, NBO charge distribution analysis in conjunction with the structural properties of the transition structures suggest their enolate-like character and the “normal” character of hydrogen bonds to the carbonyl oxygen. These energetic consequences are discussed in terms of a synergistic interplay between charge-augmented (ionic) and neutral (“normal”) H-bonds. These interactions serve to disperse this negative charge not only in the ground state but also in the TS, thereby strengthening the internal H-bonds along the entire reaction coordinate. An additional unanticipated insight, with the identification of a strong ionic C−H···carboxylate-oxygen H-bond along the reaction coordinate provides an activating influence for the α-proton abstraction. The ionic C−H···⊖OOC H-bond identified here would not only orient the general base prior to proton abstraction, but also transfers electron density to the developing enolate anion along the reaction coordinate. This strengthens the H-bonds to the thioester carbonyl oxygen that prove so critical in the acidification of the α-proton and lowers the energy of the transition state for enolization. The correlation between the reaction endothermicity, barrier height and location of the transition state along the reaction coordinate is consistent with the classical Hammond postulate. Solvent modeling performed with the COSMO//B3LYP/6-31+G(d,p) method confirms the important role of desolvation of the carboxylate anion prior to the α-proton abstraction step.

Doping high Tc superconductors with oxygen and metallic atoms: A molecular dynamics study

Using classical molecular dynamics based on Lennard–Jones-like potentials, a mechanically stable YBa2Cu3O7 high Tc superconductor structure is generated. This process is controlled via interactive computer graphics. After doping atoms into or removing atoms from the sample using a recently implemented picking mechanism, the lattice oscillation energy is annihilated with a simulated annealing procedure. The remaining minimum ground state energy allows marking of the preferred doping location. Information on the doping mechanism is important because the magnetic and superconducting properties of these compounds depend very strongly on their oxygen content.

On the lowest electronic states of the C2F radical

The adiabatic energy surfaces of the lowest three electronic states [2(2A′ and 2A′)] and 2Σ+[2A′] of the C2F radical were investigated by the Hartree-Fock multiconfiguration self-consistent field (HF—MCSCF) ab initio method using a large set of atomic natural orbitals (ANO) and an extended configuration space, and the results were shown to be in agreement with the predictions of valence theory for this radical. The electronic ground state was found to have a bent equilibrium structure, hence contradicting the Walsh rule which predicts for the isoelectronic molecules a 2 linear state. The three states were found to be nearly degenerate and the potential energy surfaces of the two lowest electronic states exhibit an avoided crossing at an energy ∼2000 cm−1 above the ground-state minimum, lower than the highest vibrational fundamental. The strong adiabatic interaction which is responsible for the ordering of the electronic states and their equilibrium geometry involves not only the bending coordinate as normally found for Renner-Teller pairs of states, but also the C—C stretching coordinate, due to the near degeneracy of the 2Σ+ and the 2 lowest electronic states at linear geometries. © 1996 John Wiley & Sons, Inc.