Theoretical study of the adsorption and dissociation of oxygen on Pt(111) in the presence of homogeneous electric fields

The effect of homogeneous electric fields on the adsorption energies of atomic and molecular oxygen and the dissociation activation energy of molecular oxygen on Pt(111) were studied by density functional theory (DFT). Positive electric fields, corresponding to positively charged surfaces, reduce the adsorption energies of the oxygen species on Pt(111), whereas negative fields increase the adsorption energies. The magnitude of the energy change for a given field is primarily determined by the static surface dipole moment induced by adsorption. On 10-atom Pt(111) clusters, the adsorption energy of atomic oxygen decreased by ca. 0.25 eV in the presence of a 0.51 V/A (0.01 au) electric field. This energy change, however, is heavily dependent on the number of atoms in the Pt(111) cluster, as the static dipole moment decreases with cluster size. Similar calculations with periodic slab models revealed a change in energy smaller by roughly an order of magnitude relative to the 10-atom cluster results. Calculations with adsorbed molecular oxygen and its transition state for dissociation showed similar behavior. Additionally, substrate relaxation in periodic slab models lowers the static dipole moment and, therefore, the effect of electric field on binding energy. The results presented in this paper indicate that the electrostatic effect of electric fields at fuel cell cathodes may be sufficiently large to influence the oxygen reduction reaction kinetics by increasing the activation energy for dissociation.     

Modeling the influence of a laser pulse on the potential energy surface in optimal molecular control theory

Understanding and modeling the interaction between light and matter is essential to the theory of optical molecular control. While the effect of the electric field on a molecule's electronic structure is often not included in control theory, it can be modeled in an optimal control algorithm by a set or toolkit of potential energy surfaces indexed by discrete values of the electric field strength where the surfaces are generated by Born-Oppenheimer electronic structure calculations that directly include the electric field. Using a new optimal control algorithm with a trigonometric mapping to limit the maximum field strength explicitly, we apply the surface-toolkit method to control the hydrogen fluoride molecule. Potential energy surfaces in the presence and absence of the electric field are created with two-electron reduced-density-matrix techniques. The population dynamics show that adjusting for changes in the electronic structure of the molecule beyond the static dipole approximation can be significant for designing a field that drives a realistic quantum system to its target observable.     

Quantum control of molecular motion including electronic polarization effects with a two-stage toolkit

A method for incorporating strong electric field polarization effects into optimal control calculations is presented. A Born-Oppenheimer-type separation, referred to as the electric-nuclear Born-Oppenheimer (ENBO) approximation, is introduced in which variations of both the nuclear geometry and the external electric field are assumed to be slow compared with the speed at which the electronic degrees of freedom respond to these changes. This assumption permits the generation of a potential energy surface that depends not only on the relative geometry of the nuclei but also on the electric field strength and on the orientation of the molecule with respect to the electric field. The range of validity of the ENBO approximation is discussed in the paper. A two-stage toolkit implementation is presented to incorporate the polarization effects and reduce the cost of the optimal control dynamics calculations. As an illustration of the method, it is applied to optimal control of vibrational excitation in a hydrogen molecule aligned along the field direction. Ab initio configuration interaction calculations with a large orbital basis set are used to compute the H-H interaction potential in the presence of the electric field. The significant computational cost reduction afforded by the toolkit implementation is demonstrated.     

强外电场作用下二甲基硅酮的分子结构和激发态

The ground states of dimethyl siloxane under different intense electric fields ranging from - 0. 04 to 0. 04 a. u. are optimized using density functional theory DFT / B3P86 at 6-311 ++ G（d,p）level. The excitation energies and oscillator strengths under the same intense applied electric fields are calculated employing the revised hybrid CIS-DFT method. The result shows that the electronic state,molecular geometry,total energy,dipole moment and excitation energy are strongly dependent on the field strength and behave asymmetry to the direction of the applied electric field. As the electric field changes from - 0. 04 to 0. 04 a. u. ,the bond length of Si-O increases whereas the bond length of Si-C decreases because of the charge transfer induced by the applied electric field. The dipole moment of the ground state decreases linearly with the applied field strength. However,the dipole moment of molecule changes from positive to negative as the inverse electric field increase to - 0. 03 a. u. Further increase of the inverse electric field results in an increase of the total energy of the molecule. The dependence of the calculated excitation energies on the applied electric field strength is fitting well to the relationship proposed by Grozema. The excitation energies of the first five excited states of dimethyl siloxane decrease as the applied electric filed increases because the energy gap between the HOMO and LUMO become close with the field,which shows that the molecule is easy to be excited under electric field and hence can be easily dissociated.     

Quantum control of molecular vibrational and rotational excitations in a homonuclear diatomic molecule: a full three-dimensional treatment with polarization forces

The optimal control of the vibrational excitation of the hydrogen molecule [Balint-Kurti et al., J. Chem. Phys. 122, 084110 (2005)] utilizing polarization forces is extended to three dimensions. The polarizability of the molecule, to first and higher orders, is accounted for using explicit ab initio calculations of the molecular electronic energy in the presence of an electric field. Optimal control theory is then used to design infrared laser pulses that selectively excite the molecule to preselected vibrational-rotational states. The amplitude of the electric field of the optimized pulses is restricted so that there is no significant ionization during the process, and a new frequency sifting method is used to simplify the frequency spectrum of the pulse. The frequency spectra of the optimized laser pulses for processes involving rotational excitation are more complex than those relating to processes involving only vibrational excitation.     

First Principles Calculations of Electric Field Effect on the (6,0) Zigzag Single-Walled Silicon Carbide Nanotube for use in Nano-Electronic Circuits

Structural, electronic, and electrical responses of the H-capped (6,0) zigzag single-walled silicon carbide nanotube (SiCNT) was studied under the parallel and transverse electric fields with strengths 0–140 × 10^?4 a.u. by using density functional calculations. Analysis of the structural parameters indicates that resistance of the nanotube against the applied parallel electric field is more than resistance of the nanotube against the applied transverse electric field. The dipole moments, atomic charge variations, and total energy of the (6,0) zigzag SiCNT show increases with any increase in the applied external electric field strengths. The length, tip diameters, electronic spatial extent, and molecular volume of the nanotube do not change significantly with any increasing in the electric field strength. The energy gap of the nanotube increases with any increases in the electric field strength and its reactivity is decreased. Increase of the ionization potential, electron affinity, chemical potential, and HOMO and LOMO in the nanotube with increase of the applied external electric field strengths indicates that the properties of SiCNTs can be controlled by the proper external electric field for use in nano-electronic circuits.     

Solvent effects on the n-->pi* electronic transition in formaldehyde: a combined coupled cluster/molecular dynamics study

We present a study of the blueshift of the n-->pi* electronic transition in formaldehyde in aqueous solution using a combined coupled cluster/molecular mechanics model including mutual polarization effects in the Hamiltonian. In addition, we report ground and excited state dipole moments. Configurations are generated from molecular dynamics simulations with two different force fields, one with and one without an explicit polarization contribution. A statistical analysis using 1200 configurations is presented. Effects of explicit polarization contributions are found to be significant. It is found that the main difference in the effects on the excitation energies arises from the fact that the two force fields result in different liquid structures, and thus a different set of configurations is generated for the coupled cluster/molecular mechanics calculations.     

Vibrational polarization and opsin shift of retinal schiff bases: theoretical study

The changes in the electronic excitation energy arising from molecular structural displacement induced by external electric field (so-called vibrational polarization) are examined theoretically for the protonated and neutral 11-cis retinal Schiff bases. It is shown that the magnitude of the field-induced structural displacement is significantly large for the protonated species, so that the change in the electronic excitation energy arising from this structural displacement is of the same order of magnitude as that arising from the direct effect of electric field on the electronic wave function. These two effects contribute additively to the field-induced spectral shift. The intensity-carrying mode (ICM) theory is employed to extract a single vibrational mode (called primary infrared ICM) that is most important for the field-induced structural displacement. A simple one-dimensional model is constructed, and the extent to which we can interpret the field-induced spectral shift by such a model is examined. In the case of the neutral species, only a small change in the electronic excitation energy is induced by external electric field, mainly because the vibrational polarizability of this species is small. The meaning of these results in the spectral tuning of visual pigments is discussed.     

外电场作用下甲基乙烯基硅酮分子结构和电子光谱

采用密度泛函B3P86方法在6-311++G(d, p)基组水平上优化得到了沿分子轴方向不同外电场(0-0.04 a.u.)作用下, 甲基乙烯基硅酮分子的基态电子状态、几何结构、电偶极矩和分子总能量. 在优化构型下利用杂化CIS-DFT方法(CIS-B3P86)研究了同样外电场条件下对甲基乙烯基硅酮的激发能和振子强度的影响. 计算结果表明, 分子几何构型与电场大小呈现强烈的依赖关系, 分子偶极矩μ随电场的增加先减小后急剧增大. 电场为零时, 分子总能量为-483.5532137 a.u., 随着电场增加, 能量升高, 在F=0.02 a.u.时达到最大值-483.5393952 a.u., 此后, 继续增大电场系统总能量则开始降低. 激发能随电场增加急剧减小, 表明在电场作用下, 分子易于激发和离解.     

电场对(4, 0)Zigzag模型单壁碳纳米管的影响

The structural and electronic properties of a (4, 0) zigzag single-walled carbon nanotube (SWCNT) under parallel and transverse electric fields with strengths of 0-1.4×10~(-2) a.u. Were studied using the density functional theory (DFT) B3LYP/6-31G~* method. Results show that the properties of the SWCNT are dependent on the external electric field. The applied external electric field strongly affects the molecular dipole moments. The induced dipole moments increase linearly with increase in the electrical field intensities. This study shows that the application of parallel and transverse electric fields results in changes in the occupied and virtual molecular orbitals (Mos) but the energy gap between the highest occupied MO (HOMO) and the lowest unoccupied MO (LUMO) of this SWCNT is less sensitive to the electric field strength. The electronic spatial extent (ESE) and length of the SWCNT show small changes over the entire range of the applied electric field strengths. The natural bond orbital (NBO) electric charges on the atoms of the SWCNT show that increase in the external electric field strength increases the separation of the center of the positive and negative electric charges of the carbon nanotube.     

Computational Insights into the Role of External and Local Electric Fields in Macrocyclic Chemical and Biological Systems

The investigation of the role of the electric field in systems of widespread interest employing computational techniques is an emerging area of research. The outcome of applying an oriented external electric field (OEEF) on the geometric and electronic properties of the chemically unique π-conjugated cyclic carbon ring compounds has been explored with density functional theory (DFT). Distinct changes in the structural and electronic features of such ring compounds are observed upon the application of OEEFs. Importantly, the calculations indicate that a mixed aliphatic-aromatic conjugated ring converts from a singlet to a triplet after the application of an OEEF, suggesting potential applications in optoelectronics for such molecules, without the need for photochemically induced change in the spin state. Furthermore, the influence of built-in local electric fields (LEFs) present in naturally occurring macrocyclic systems such as valinomycin has also been explored. Static and ab initio molecular dynamics (AIMD) calculations indicate that LEFs are the primary driving factor in determining the energetically favoured position of counter anions such as chloride (Cl⁻) in the potassium (K⁺) and sodium (Na⁺) coordinated valinomycin macrocycle structures: they exist inside the cage in the case of K⁺ sequestration by valinomycin and outside for Na⁺. This divergence has been proposed to be the determining factor for the selectivity of the valinomycin macrocycle for binding a K⁺ cation over Na⁺.     

Electronic structures of organic molecule encapsulated BN nanotubes under transverse electric field

The electronic structures of boron nitride nanotubes (BNNTs) doped by different organic molecules under a transverse electric field were investigated via first-principles calculations. The external field reduces the energy gap of BNNT, thus makes the molecular bands closer to the BNNT band edges and enhances the charge transfers between BNNT and molecules. The effects of the electric field direction on the band structure are negligible. The electric field shielding effect of BNNT to the inside organic molecules is discussed. Organic molecule doping strongly modifies the optical property of BNNT, and the absorption edge is redshifted under static transverse electric field.     

On the performance of molecular polarization methods. II. Water and carbon tetrachloride close to a cation

Our initial study on the performance of molecular polarization methods close to a positive point charge [M. Masia, M. Probst, and R. Rey, J. Chem. Phys. 121, 7362 (2004)] is extended to the case in which a molecule interacts with a real cation. Two different methods (point dipoles and shell model) are applied to both the ion and the molecule. The results are tested against high-level ab initio calculations for a molecule (water or carbon tetrachloride) close to Li+, Na+, Mg2+, and Ca2+. The monitored observable is in all cases the dimer electric dipole as a function of the ion-molecule distance for selected molecular orientations. The moderate disagreement previously obtained for point charges at intermediate distances, and attributed to the linearity of current polarization methods (as opposed to the nonlinear effects evident in ab initio calculations), is confirmed for real cations as well. More importantly, it is found that at short separations the phenomenological polarization methods studied here substantially overestimate the dipole moment induced if the ion is described quantum chemically as well, in contrast to the dipole moment induced by a point-charge ion, for which they show a better degree of accord with ab initio results. Such behavior can be understood in terms of a decrease of atomic polarizabilities due to the repulsion between electronic charge distributions at contact separations. It is shown that a reparametrization of the Thole method for damping of the electric field, used in conjunction with any polarization scheme, allows to satisfactorily reproduce the dimer dipole at short distances. In contrast with the original approach (developed for intramolecular interactions), the present reparametrization is ion and method dependent, and corresponding parameters are given for each case.     

Towards a force field based on density fitting

Total intermolecular interaction energies are determined with a first version of the Gaussian electrostatic model (GEM-0), a force field based on a density fitting approach using s-type Gaussian functions. The total interaction energy is computed in the spirit of the sum of interacting fragment ab initio (SIBFA) force field by separately evaluating each one of its components: electrostatic (Coulomb), exchange repulsion, polarization, and charge transfer intermolecular interaction energies, in order to reproduce reference constrained space orbital variation (CSOV) energy decomposition calculations at the B3LYP/aug-cc-pVTZ level. The use of an auxiliary basis set restricted to spherical Gaussian functions facilitates the rotation of the fitted densities of rigid fragments and enables a fast and accurate density fitting evaluation of Coulomb and exchange-repulsion energy, the latter using the overlap model introduced by Wheatley and Price [Mol. Phys. 69, 50718 (1990)]. The SIBFA energy scheme for polarization and charge transfer has been implemented using the electric fields and electrostatic potentials generated by the fitted densities. GEM-0 has been tested on ten stationary points of the water dimer potential energy surface and on three water clusters (n = 16,20,64). The results show very good agreement with density functional theory calculations, reproducing the individual CSOV energy contributions for a given interaction as well as the B3LYP total interaction energies with errors below kBT at room temperature. Preliminary results for Coulomb and exchange-repulsion energies of metal cation complexes and coupled cluster singles doubles electron densities are discussed.     

Photoinduced dissociation of water and transport of hydrogen between silver clusters

Theoretical electronic structure calculations are reported for the dissociation of water adsorbed on a 31-atom silver cluster, Ag31, and subsequent transfer of a H to a second Ag31 cluster leaving OH on the first cluster. Both ground and excited electronic state processes are considered for two choices of Ag cluster separation, 6.35 and 7.94 A, on the basis of preliminary calculations for a range of separation distances. The excited electronic state of interest is formed by photoemission of an electron from one Ag cluster and transient attachment of the photoemitted electron to the adsorbed water molecule. A very large energy barrier is found for the ground-state process (3.53 eV at a cluster separation of 6.35 A), while the barrier in the excited state is small (0.38 eV at a cluster separation of 6.35 A). In the excited state, partial occupancy of an OH antibonding orbital facilitates OH stretch and concomitant movement of the negatively charged OH toward the electron-hole in the metal cluster. The excited-state pathway for dissociation of water and transfer of H begins with the formation of an excited electronic state at 3.59-3.82 eV. Stretch of the OH bond occurs with little change in energy (0.38-0.54 eV up to a stretch of 1.96 A). In this region of OH stretch the molecule must return to the ground-state potential energy surface to fully dissociate and to transfer H to the other Ag cluster. Geometry optimizations are carried out using a simplex algorithm and a semigrid method. These methods allow the total energy to be calculated directly using configuration interaction theory.     

Ab initio theoretical study of substituted dicarboxylic acids adsorbed on GaAs surfaces: correlation between microscopic properties and observed electrical behavior

Five substituted tartaric acid derivatives are studied using density functional theory, both isolated and adsorbed onto an oxidized GaAs cluster, to model molecular layers on semiconductor surfaces. The structures, energies, and electronic properties are computed to clarify the interactions responsible for the electric behavior of the modified surfaces, used in semiconductor/metal junction devices. The chemical structure of the molecule/GaAs adducts is optimized ab initio and discussed for the first time. A strong binding scheme is found, providing useful insights about the microscopic structure of the molecular layer. A widely used model based on molecular dipole layers is discussed and verified, by computing the dipole moment for the isolated systems and estimating the charge separation in the adducts; moreover the molecular orbitals energies are analyzed and correlated to the experimental measures of the modified surface electron affinity.     

Binding at molecule/gold transport interfaces. V. Comparison of different metals and molecular bridges

The geometric and electronic structural properties of symmetric and asymmetric metal cluster-molecule-cluster' complexes have been explored. The metals include Au, Ag, Pd, and Al, and both benzenedithiol and the three isometric forms of dicyanobenzene are included as bridging molecules. Calculated properties such as cluster-molecule interface geometry, electronic state, degree of metal --> molecule charge transfer, metal-molecule mixing in the highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) energy region, the HOMO-LUMO gap, cluster --> cluster' charge transfer as a function of external field strength and direction, and the form of the potential profile across such complexes have been examined. Attempts are made to correlate charge transport with the characteristics of the cluster-complex systems. Indications of rectification in complexes that are asymmetric in the molecule, clusters, and molecule-cluster interfaces are discussed. The results obtained here are only suggestive because of the limitations of the cluster-complex model as it relates to charge transport.