Visualization of Molecular Orbitals: Formaldehyde
ERIC Educational Resources Information Center
Olcott, Richard J.
1972-01-01
Describes a computer program that plots a solid" representation of molecular orbital charge density which can be used to analyze wave functions of molecules. Illustrated with diagrams for formaldehyde. (AL)
ISS EPS Orbital Replacement Unit Block Diagrams
NASA Technical Reports Server (NTRS)
Schmitz, Gregory V.
2001-01-01
The attached documents are being provided to Switching Power Magazine for information purposes. This magazine is writing a feature article on the International Space Station Electrical Power System, focusing on the switching power processors. These units include the DC-DC Converter Unit (DDCU), the Bi-directional Charge/Discharge Unit (BCDU), and the Sequential Shunt Unit (SSU). These diagrams are high-level schematics/block diagrams depicting the overall functionality of each unit.
Stability diagrams for a rigid gyrostat in a circular orbit
NASA Astrophysics Data System (ADS)
Hughes, P. C.; Golla, D. F.
1984-05-01
The dynamics of a rigid gyrostat travelling in a circular orbit are reviewed. The two torques to be balanced are the gyroscopic torque associated with the once per orbit precession of the in plane angular momentum vector, and the gravitational torque that results from the principal inertial axes not being aligned with the orbiting axes. Different methods of balancing these two torques lead to different classes of relative equilibria, viz., the cylindrical case, the conical case, the hyperbolic case, and the general case. The principal contribution is to present a large number of stability diagrams for the many different cases and subcases that arise. A large number of equilbria are, in fact, unstable. However, as shown here, they can usually be stabilized by adding a momentum wheel aligned with the pitch axis.
Floating orbital molecular dynamics simulations.
Perlt, Eva; Brüssel, Marc; Kirchner, Barbara
2014-04-21
We introduce an alternative ab initio molecular dynamics simulation as a unification of Hartree-Fock molecular dynamics and the floating orbital approach. The general scheme of the floating orbital molecular dynamics method is presented. Moreover, a simple but sophisticated guess for the orbital centers is provided to reduce the number of electronic structure optimization steps at each molecular dynamics step. The conservation of total energy and angular momentum is investigated in order to validate the floating orbital molecular dynamics approach with and without application of the initial guess. Finally, a water monomer and a water dimer are simulated, and the influence of the orbital floating on certain properties like the dipole moment is investigated. PMID:24600690
Imaging molecular orbitals using photoionization
NASA Astrophysics Data System (ADS)
Santra, Robin
2006-10-01
The interpretation of a recent experiment using high-order harmonic generation [Itatani et al., Nature 432 (2004) 867] as a measurement of the highest occupied molecular orbital of a molecule is conceptually problematic, even if the independent-particle picture is taken seriously. Guided by the relationship between the amplitude for one-photon-induced electron emission and the electron-ion recombination amplitude in the three-step model of high-order harmonic generation, it is argued that synchrotron-based photoionization might be a superior approach to imaging molecular orbitals. Within the Hartree-Fock independent-particle picture, the molecular-frame photoelectron angular distributions, measured as a function of photon energy, could be used to reconstruct all orbitals occupied in the Hartree-Fock ground state of the molecule investigated. It is suggested that laser alignment techniques could be employed to facilitate the measurement of the molecular-frame photoelectron angular distributions.
Ab initio molecular crystal structures, spectra, and phase diagrams.
Hirata, So; Gilliard, Kandis; He, Xiao; Li, Jinjin; Sode, Olaseni
2014-09-16
Conspectus Molecular crystals are chemists' solids in the sense that their structures and properties can be understood in terms of those of the constituent molecules merely perturbed by a crystalline environment. They form a large and important class of solids including ices of atmospheric species, drugs, explosives, and even some organic optoelectronic materials and supramolecular assemblies. Recently, surprisingly simple yet extremely efficient, versatile, easily implemented, and systematically accurate electronic structure methods for molecular crystals have been developed. The methods, collectively referred to as the embedded-fragment scheme, divide a crystal into monomers and overlapping dimers and apply modern molecular electronic structure methods and software to these fragments of the crystal that are embedded in a self-consistently determined crystalline electrostatic field. They enable facile applications of accurate but otherwise prohibitively expensive ab initio molecular orbital theories such as Møller-Plesset perturbation and coupled-cluster theories to a broad range of properties of solids such as internal energies, enthalpies, structures, equation of state, phonon dispersion curves and density of states, infrared and Raman spectra (including band intensities and sometimes anharmonic effects), inelastic neutron scattering spectra, heat capacities, Gibbs energies, and phase diagrams, while accounting for many-body electrostatic (namely, induction or polarization) effects as well as two-body exchange and dispersion interactions from first principles. They can fundamentally alter the role of computing in the studies of molecular crystals in the same way ab initio molecular orbital theories have transformed research practices in gas-phase physical chemistry and synthetic chemistry in the last half century. In this Account, after a brief summary of formalisms and algorithms, we discuss applications of these methods performed in our group as compelling
Classification of hyperbolic Dynkin diagrams, root lengths and Weyl group orbits
NASA Astrophysics Data System (ADS)
Carbone, Lisa; Chung, Sjuvon; Cobbs, Leigh; McRae, Robert; Nandi, Debajyoti; Naqvi, Yusra; Penta, Diego
2010-04-01
We give a criterion for a Dynkin diagram, equivalently a generalized Cartan matrix, to be symmetrizable. This criterion is easily checked on the Dynkin diagram. We obtain a simple proof that the maximal rank of a Dynkin diagram of compact hyperbolic type is 5, while the maximal rank of a symmetrizable Dynkin diagram of compact hyperbolic type is 4. Building on earlier classification results of Kac, Kobayashi-Morita, Li and Saçlio\\skew3\\tildeg lu, we present the 238 hyperbolic Dynkin diagrams in ranks 3-10, 142 of which are symmetrizable. For each symmetrizable hyperbolic generalized Cartan matrix, we give a symmetrization and hence the distinct lengths of real roots in the corresponding root system. For each such hyperbolic root system we determine the disjoint orbits of the action of the Weyl group on real roots. It follows that the maximal number of disjoint Weyl group orbits on real roots in a hyperbolic root system is 4.
How Different Variants of Orbit Diagrams Influence Student Explanations of the Seasons
ERIC Educational Resources Information Center
Lee, Victor R.
2010-01-01
The cause of the seasons is often associated with a very particular alternative conception: That the earth's orbit around the sun is highly elongated, and the differences in distance result in variations in temperature. It has been suggested that the standard diagrams used to depict the earth's orbit may be in some way responsible for the initial…
Orbital effect of magnetic field on the Majorana phase diagram
NASA Astrophysics Data System (ADS)
Nijholt, Bas; Akhmerov, Anton R.
2016-06-01
Studies of Majorana bound states in semiconducting nanowires frequently neglect the orbital effect of a magnetic field. Systematically studying its role leads us to several conclusions for designing Majoranas in this system. Specifically, we show that for experimentally relevant parameter values the orbital effect of a magnetic field has a stronger impact on the dispersion relation than the Zeeman effect. While Majoranas do not require the presence of only one dispersion subband, we observe that the size of the Majoranas becomes unpractically large, and the band gap unpractically small, when more than one subband is filled. Since the orbital effect of a magnetic field breaks several symmetries of the Hamiltonian, it leads to the appearance of large regions in parameter space with no band gap whenever the magnetic field is not aligned with the wire axis. The reflection symmetry of the Hamiltonian with respect to the plane perpendicular to the wire axis guarantees that the wire stays gapped in the topologically nontrivial region as long as the field is aligned with the wire.
Molecular Complexation and Phase Diagrams of Urea/PEG Mixtures
NASA Astrophysics Data System (ADS)
Fu, Guoepeng; Kyu, Thein
2014-03-01
Polyethylene glycol (PEG) and urea complexation has been known to form a stable crystal due to molecular complexation. The effect of molecular weight of PEG on the phase diagrams of its blends with urea has been explored. In the case of high molecular weight PEG8k/urea, the observed phase diagram is azeotrope, accompanied by eutectoid reactions in the submerged phases such as induced stable ``alpha'' phase crystals and metastable ``beta'' phase crystals. The metastable crystal can transform to stable crystal under a certain thermal annealing condition. However, the phase diagram of PEG1k/urea is of coexistence loop, whereas PEG400/urea exhibits eutectic character. Subsequently, the change of azeotrope to eutectic behavior with PEG molecular weight is analyzed in the context of the combined Flory-Huggins theory of liquid-liquid demixing and phase field theory of crystal solidification. Of particular interest is that only a very small urea amount (2 wt%) is needed to form a stable inclusion crystal via complexation with PEG. Potential application in lithium battery is discussed based on AC impedance spectroscopy and cyclic voltammetry. Supported by NSF-DMR 1161070.
Some Observations on Molecular Orbital Theory
ERIC Educational Resources Information Center
Journal of Chemical Education, 2005
2005-01-01
A few flawed predictions in the context of homonuclear diatomic molecules are presented in order to introduce students to molecular orbital (MO) theory. A common misrepresentation of the relationship between the energy of an atomic orbital and the energy of the MO associated with the atomic orbital is illustrated.
Argand diagram representation of orbiting resonance in proton-transfer collision
NASA Astrophysics Data System (ADS)
Lutrus, C. K.; Suck Salk, S. H.
1991-05-01
Dynamic resonance in atom-diatomic molecule collisions has been relatively well studied compared to orbiting resonance. We discuss orbiting resonance on reactive scattering involving proton (charge) transfer. Resonance structure is predicted to exist at forward-scattering angles in both the state-to-state angular distribution and the Argand diagram for the proton-transfer collision system of He+H+2-->HeH++H. The present study demonstrates the possibility of orbiting resonance particularly in proton (charge) -transfer reaction involving atom-diatomic molecule systems.
ERIC Educational Resources Information Center
Bo¨rgel, Jonas; Campbell, Michael G.; Ritter, Tobias
2016-01-01
The presentation of d-orbital splitting diagrams for square planar transition metal complexes in textbooks and educational materials is often inconsistent and therefore confusing for students. Here we provide a concise summary of the key features of orbital splitting diagrams for square planar complexes, which we propose may be used as an updated…
A Simple Huckel Molecular Orbital Plotter
ERIC Educational Resources Information Center
Ramakrishnan, Raghunathan
2013-01-01
A program is described and presented to readily plot the molecular orbitals from a Huckel calculation. The main features of the program and the scope of its applicability are discussed through some example organic molecules. (Contains 2 figures.)
ERIC Educational Resources Information Center
Orenha, Renato P.; Galembeck, Sérgio E.
2014-01-01
This computational experiment presents qualitative molecular orbital (QMO) and computational quantum chemistry exercises of NO, NO[superscript+], and NO[superscript-]. Initially students explore several properties of the target molecules by Lewis diagrams and the QMO theory. Then, they compare qualitative conclusions with EHT and DFT calculations…
Exact phase diagram of multi-orbital spin-fermion model for hole doped cuprates
NASA Astrophysics Data System (ADS)
Meetei, Oinam Nganba; Fischer, Mark; Lawler, Michael; Paramekanti, Arun; Kim, Eun-Ah
2015-03-01
Recent experiments revealing the ubiquitous presence of spin and charge ordered states in hole-doped cuprates have placed the study of broken symmetry states at the center of high Tc superconductivity research. Here we aim to understand the phase diagram of broken symmetry states using a simple model that captures the essence of hole doped cuprates. The model consists of itinerant quantum holes on oxygen p-orbitals coupled to classical Cu spins. It is amenable to sign problem free Monte-Carlo simulation allowing us to study finite temperature properties as well as unbiased determination of ground state spin and charge configuration. As a function of system parameters, we obtain a rich phase diagram. Our analysis provides a transparent and unifying picture for various charge and spin ordered states as arising from frustration of antiferromagnetic order due to hole doping, through exact finite temperature phase diagram of the model. Supported by the U.S. Department of Energy Grant DE-SC0010313.
Spin–orbit interaction mediated molecular dissociation
Kokkonen, E. Jänkälä, K.; Kettunen, J. A.; Heinäsmäki, S.; Karpenko, A.; Huttula, M.; Löytynoja, T.
2014-05-14
The effect of the spin–orbit interaction to photofragmentation is investigated in the mercury(II) bromide (HgBr{sub 2}) molecule. Changes in the fragmentation between the two spin–orbit components of Hg 5d photoionization, as well as within the molecular-field-splitted levels of these components are observed. Dissociation subsequent to photoionization is studied with synchrotron radiation and photoelectron-photoion coincidence spectroscopy. The experimental results are accompanied by relativistic ab initio analysis of the photoelectron spectrum.
Simple orbital theory for the molecular electrician.
Ernzerhof, Matthias
2011-07-01
Theories of molecular electronic devices (MEDs) are quite involved in general. However, various prominent features of MEDs can be understood drawing only on elementary quantum theory. To support this point of view, we provide a two component orbital theory that enables one to reproduce various important features of MEDs. In this theory, the device orbitals are divided into two components, each of which is obtained from simple rules. To illustrate our two-component model, we apply it to explain, among other things, the conductance suppression in cross-conjugated systems and the dependence of the conductance on the contact position in aromatic systems. PMID:21744885
Phase diagram and collective excitations in an excitonic insulator from an orbital physics viewpoint
NASA Astrophysics Data System (ADS)
Nasu, Joji; Watanabe, Tsutomu; Naka, Makoto; Ishihara, Sumio
2016-05-01
An excitonic-insulating system is studied from a viewpoint of the orbital physics in strongly correlated electron systems. An effective model Hamiltonian for low-energy electronic states is derived from the two-orbital Hubbard model with a finite-energy difference corresponding to the crystalline-field splitting. The effective model is represented by the spin operators and the pseudospin operators for the spin-state degrees of freedom. The ground-state phase diagram is analyzed by the mean-field approximation. In addition to the low-spin state and high-spin state phases, two kinds of the excitonic-insulating phases emerge as a consequence of the competition between the crystalline-field effect and the Hund coupling. Transitions to the excitonic phases are classified to an Ising-type transition resulted from a spontaneous breaking of the Z2 symmetry. Magnetic structures in the two excitonic-insulating phases are different from each other: an antiferromagnetic order and a spin nematic order. Collective excitations in each phase are examined using the generalized spin-wave approximation. Characteristics in the Goldstone modes in the excitonic-insulating phases are studied through the calculations of the dynamical correlation functions for the spins and pseudospin operators. Both the transverse and longitudinal spin excitation modes are active in the two excitonic-insulating phases in contrast to the low-spin state and high-spin state phases. Relationships of the present results to the perovskite cobalt oxides are discussed.
A Comparison of Molecular Vibrational Theory to Huckel Molecular Orbital Theory.
ERIC Educational Resources Information Center
Keeports, David
1986-01-01
Compares the similar mathematical problems of molecular vibrational calculations (at any intermediate level of sophistication) and molecular orbital calculations (at the Huckel level). Discusses how the generalizations of Huckel treatment of molecular orbitals apply to vibrational theory. (TW)
Rotation and Anisotropic Molecular Orbital Effect in a Single H2TPP Molecule Transistor
NASA Astrophysics Data System (ADS)
Sakata, Shuichi; Yoshida, Kenji; Kitagawa, Yuichi; Ishii, Kazuyuki; Hirakawa, Kazuhiko
2013-12-01
Electron transport through a single molecule is determined not only by the intrinsic properties of the molecule but also by the configuration of the molecule with respect to the lead electrodes. Here, we show how electron transport through a single H2TPP molecule is modulated by changes in the configuration. The Coulomb stability diagram of a single H2TPP molecule transistor exhibited a few different patterns in different measurement scans. Furthermore, the sample exhibited negative differential resistance, the magnitude of which changed with the pattern in the Coulomb stability diagram. Such behavior can be explained by the rotation of the molecule with anisotropic molecular orbitals in the gap electrodes induced by electrical stress. Moreover, we find that the energy separations between molecular orbitals are also affected by the rotation, confirming that the metal-molecule interface configuration renormalizes the electronic levels in the molecule.
Meyer, Benjamin; Guillot, Benoît; Ruiz-Lopez, Manuel F; Genoni, Alessandro
2016-03-01
Despite more and more remarkable computational ab initio results are nowadays continuously obtained for large macromolecular systems, the development of new linear-scaling techniques is still an open and stimulating field of research in theoretical chemistry. In this family of methods, an important role is occupied by those strategies based on the observation that molecules are generally constituted by recurrent functional units with well-defined intrinsic features. In this context, we propose to exploit the notion of extremely localized molecular orbitals (ELMOs) that, due to their strict localization on small molecular fragments (e.g., atoms, bonds, or functional groups), are in principle transferable from one molecule to another. Accordingly, the construction of orbital libraries to almost instantaneously build up approximate wave functions and electron densities of very large systems becomes conceivable. In this work, the ELMOs transferability is further investigated in detail and, furthermore, suitable rules to construct model molecules for the computation of ELMOs to be stored in future databanks are also defined. The obtained results confirm the reliable transferability of the ELMOs and show that electron densities obtained from the transfer of extremely localized molecular orbitals are very close to the corresponding Hartree-Fock ones. These observations prompt us to construct new ELMOs databases that could represent an alternative/complement to the already popular pseudoatoms databanks both for determining electron densities and for refining crystallographic structures of very large molecules. PMID:26799516
Molecular orbitals for properties and spectroscopies
Robert, Vincent; Domingo, Alex; Braunstein, Pierre; Danopoulos, Andreas; Monakhov, Kirill
2015-12-31
The description and clarification of spectroscopies and properties goes through ab initio calculations. Wave function based calculations (CASSCF/CASPT2) are particularly appealing since they offer spectroscopic accuracy and means of interpretation. we performed such calculations to elucidate the origin of unusual structural changes and intramolecular electron transfer phenomenon. Based on optimized molecular orbitals and a reading of the multireference wave function, it is suggested that intimate interactions are likely to considerably modify the standard pictures. A so-called PIMA (polarization-induced metalâĹŠarene) interaction similar to the more familiar anion-π interaction is responsible for a significant deviation from sp{sup 3} geometry and an energetic stabilization of 50 kJ/mol in Cr(II) benzyl organometallic complexes. In a similar fashion, it is proposed that the energetic profile of the IVCT (inter valence charge transfer) exhibits strong similarities to the Marcus’ theory, suggesting a response behaviour of the ensemble of electrons as electron transfer occurs in Fe{sup 2+}/Fe{sup 3+} bimetallic compound. The electronic reorganization induced by the IVCT process accounts for 11.8 eV, a very large effect that reduces the transfer energy down to 0.89 eV, in very good agreement with experiments.
Molecular-orbital model for metal-sapphire interfacial strength
NASA Technical Reports Server (NTRS)
Johnson, K. H.; Pepper, S. V.
1982-01-01
Self-consistent-field X-Alpha scattered-wave cluster molecular-orbital models have been constructed for transition and noble metals (Fe, Ni, Cu, and Ag) in contact with a sapphire (Al2O3) surface. It is found that a chemical bond is established between the metal d-orbital electrons and the nonbonding 2p-orbital electrons of the oxygen anions on the Al2O3 surface. An increasing number of occupied metal-sapphire antibonding molecular orbitals explains qualitatively the observed decrease of contact shear strength through the series Fe, Ni, Cu, and Ag.
Orbital Energy Levels in Molecular Hydrogen. A Simple Approach.
ERIC Educational Resources Information Center
Willis, Christopher J.
1988-01-01
Described are the energetics involved in the formation of molecular hydrogen using concepts that should be familiar to students beginning the study of molecular orbital theory. Emphasized are experimental data on ionization energies. Included are two-electron atomic and molecular systems. (CW)
BetaVoid: molecular voids via beta-complexes and Voronoi diagrams.
Kim, Jae-Kwan; Cho, Youngsong; Laskowski, Roman A; Ryu, Seong Eon; Sugihara, Kokichi; Kim, Deok-Soo
2014-09-01
Molecular external structure is important for molecular function, with voids on the surface and interior being one of the most important features. Hence, recognition of molecular voids and accurate computation of their geometrical properties, such as volume, area and topology, are crucial, yet most popular algorithms are based on the crude use of sampling points and thus are approximations even with a significant amount of computation. In this article, we propose an analytic approach to the problem using the Voronoi diagram of atoms and the beta-complex. The correctness and efficiency of the proposed algorithm is mathematically proved and experimentally verified. The benchmark test clearly shows the superiority of BetaVoid to two popular programs: VOIDOO and CASTp. The proposed algorithm is implemented in the BetaVoid program which is freely available at the Voronoi Diagram Research Center (http://voronoi.hanyang.ac.kr). PMID:24677176
Periodic orbits of the hydrogen molecular ion and their quantization
Duan, Y.; Yuan, J.; Bao, C.
1995-11-01
In a classical study of the hydrogen molecular ion beyond the Born-Oppenheimer approximation (BOA), we have found that segments of trajectories resemble that of the Born-Oppenheimer approximation periodic orbits. The importance of this fact to the classical understanding of chemical bonding leads us to a systematic study of the periodic orbits of the planar hydrogen molecular ion within the BOA. Besides introducing a classification scheme for periodic orbits, we discuss the convergence properties of families of periodic orbits and their bifurcation patterns according to their types. Semiclassical calculations of the density of states based on these periodic orbits yield results in agreement with the exact quantum eigenvalues of the hydrogen molecular ion system.
NASA Astrophysics Data System (ADS)
Dong, Dong; Zou, Xu-Bo; Guo, Guang-Can
2015-07-01
We studied the system of pure Rashba spin-orbit coupled Bose gas with an in-plane magnetic field. Based on the mean field theory, we obtained the zero temperature phase diagram of the system which exhibits three phases, plane wave (PW) phase, striped wave (SW) phase, and zero momentum (ZM) phase. It was shown that with a growing in-plane field, both SW and ZM phases will eventually turn into the PW phase. Furthermore, we adopted the Bogoliubov theory to study the excitation spectrum as well as the sound speed. Project supported by the National Natural Science Foundation of China (Grant No. 10774088).
NASA Technical Reports Server (NTRS)
Adler, David S.; Roberts, William W., Jr.
1992-01-01
Techniques which use longitude-velocity diagrams to identify molecular cloud complexes in the disk of the Galaxy are investigated by means of model Galactic disks generated from N-body cloud-particle simulations. A procedure similar to the method used to reduce the low-level emission in Galactic l-v diagrams is employed to isolate complexes of emission in the model l-v diagram (LVCs) from the 'background'clouds. The LVCs produced in this manner yield a size-line-width relationship with a slope of 0.58 and a mass spectrum with a slope of 1.55, consistent with Galactic observations. It is demonstrated that associations identified as LVCs are often chance superpositions of clouds spread out along the line of sight in the disk of the model system. This indicates that the l-v diagram cannot be used to unambiguously determine the location of molecular cloud complexes in the model Galactic disk. The modeling results also indicate that the existence of a size-line-width relationship is not a reliable indicator of the physical nature of cloud complexes, in particular, whether the complexes are gravitationally bound objects.
Local Molecular Orbitals from a Projection onto Localized Centers.
Heßelmann, Andreas
2016-06-14
A localization method for molecular orbitals is presented which exploits the locality of the eigenfunctions associated with the largest eigenvalues of the matrix representation of spatially localized functions. Local molecular orbitals are obtained by a projection of the canonical orbitals onto the set of the eigenvectors which correspond to the largest eigenvalues of these matrices. Two different types of spatially localized functions were chosen in this work, a two-parameter smooth-step-type function and the weight functions determined by a Hirshfeld partitioning of the molecular volume. It is shown that the method can provide fairly local occupied molecular orbitals if the positions of the set of local functions are set to the molecular bond centers. The method can also yield reasonably well-localized virtual molecular orbitals, but here, a sensible choice of the positions of the functions are the atomic sites and the locality then depends more strongly on the shape of the set of local functions. The method is tested for a range of polypeptide molecules in two different conformations, namely, a helical and a β-sheet conformation. Futhermore, it is shown that an adequate locality of the occupied and virtual orbitals can also be obtained for highly delocalized systems. PMID:27164445
Molecular shield - An orbiting low-density materials laboratory
NASA Technical Reports Server (NTRS)
Melfi, L. T., Jr.; Outlaw, R. A.; Hueser, J. E.; Brock, F. J.
1976-01-01
Analysis of a molecular shield orbited at 200 km utilizes the kinetic theory of a drifting Maxwellian gas, applied to a hemispherical shell geometry containing internal sources. The molecular shield provides very low gas density conditions for materials experiments at low gravity, while the hemispherical geometry minimizes the internal surface/volume ratio. Deployment of the shield in orbit is described. Contributions to density by shield outgassing, by experiment outgassing, and by interaction with the orbiter are discussed separately. A jettisonable closure plate sealing the hemisphere minimizes any risk of experiment contamination during deployment.
Molecular orbital analysis of dicarbido-transition-metal cluster compounds
Halet, J.; Mingos, D.M.P.
1988-01-01
Molecular orbital calculations on dicarbido-transition-metal carbonyl cluster compounds have shown that the bonding between C/sub 2/ and the metal cage results primarily from electron donation from the C/sub 2/ sigma/sub rho/- and ..pi..-bonding molecular orbitals and back donation from filled metallic molecular orbitals to the C/sub 2/ ..pi..* orbitals. The bonding therefore follows closely the Chatt-Dewar-Ducanson model that has been established previously for ethyne and ethene complexes but not for interstitial moieties. The C-C separation in the dicarbido clusters depends critically on the geometric constraints imposed by the metal cage and the extent of forward and back donation. In these clusters where the carbon atoms are in adjacent trigonal-prismatic sites the calculated formal bond order is between 1.0 and 1.5, which agrees well with the observed C-C bond lengths.
Ground-state phase diagram of a spin-orbit-coupled bosonic superfluid in an optical lattice
NASA Astrophysics Data System (ADS)
Chen, Zhu; Liang, Zhaoxin
2016-01-01
In recent experiments, spin-orbit-coupled (SOC) bosonic gases in an optical lattice have been successfully prepared into any Bloch band [Hamner et al., Phys. Rev. Lett. 114, 070401 (2015), 10.1103/PhysRevLett.114.070401], which promises a viable contender in the competitive field of simulating gauge-related phenomena. However, the ground-state phase diagram of such systems in the superfluid regime is still lacking. Here we present a detailed study of the phase diagram in an optically trapped Bose gas with equal-weight Rashba and Dresselhaus SO coupling. We identify four different quantum phases, which include three normal phases and a mixed phase, by considering the wave vector k1, the longitudinal <σz> , and the transverse <σx> spin polarizations as three order parameters. The ground state of normal phases is a Bloch wave with a single wave vector k1, which can position in arbitrary regions in the Brillouin zone. By contrast, the ground state of the mixed phase is a superposition of two Bloch waves with opposite k1, which, remarkably, may lack periodicity even though the system's Hamiltonian is periodic. This mixed phase in the lattice setting can be seen as the counterpart of the stripe phase associated with the uniform SOC gas. Furthermore, due to the lattice-renormalized SOC, the phase diagram of the model system becomes significantly different from the uniform case when the lattice strength grows. Finally, a scheme for experimentally probing the mixed phase using Bragg spectroscopy is proposed.
Exploring chemistry with the fragment molecular orbital method.
Fedorov, Dmitri G; Nagata, Takeshi; Kitaura, Kazuo
2012-06-01
The fragment molecular orbital (FMO) method makes possible nearly linear scaling calculations of large molecular systems, such as water clusters, proteins and DNA. In particular, FMO has been widely used in biochemical applications involving protein-ligand binding and drug design. The method has been efficiently parallelized suitable for petascale computing. Many commonly used wave functions and solvent models have been interfaced with FMO. We review the historical background of FMO, and summarize its method development and applications. PMID:22410762
Coulomb-corrected molecular orbital tomography of nitrogen
Zhai, Chunyang; He, Lixin; Lan, Pengfei; Zhu, Xiaosong; Li, Yang; Wang, Feng; Shi, Wenjing; Zhang, Qingbin; Lu, Peixiang
2016-01-01
High-order harmonic generation (HHG) from aligned molecules has provided a promising way to probe the molecular orbital with an Ångström resolution. This method, usually called molecular orbital tomography (MOT) replies on a simple assumption of the plane-wave approximation (PW), which has long been questioned due to that PW approximation is known to be valid in the keV energy region. However, the photon energy is usually no more than 100 eV in HHG. In this work, we experimentally reconstruct the highest occupied molecular orbital (HOMO) of nitrogen (N2) by using a Coulomb-corrected MOT (CCMOT) method. In our scheme, the molecular continuum states are described by a Coulomb wave function instead of the PW approximation. With CCMOT, the reconstructed orbital is demonstrated to agree well with the theoretical prediction and retain the main features of the HOMO of N2. Compared to the PW approximation method, the CCMOT shows a significant improvement in eliminating the artificial structures caused by PW approximation. PMID:27000666
Coulomb-corrected molecular orbital tomography of nitrogen
NASA Astrophysics Data System (ADS)
Zhai, Chunyang; He, Lixin; Lan, Pengfei; Zhu, Xiaosong; Li, Yang; Wang, Feng; Shi, Wenjing; Zhang, Qingbin; Lu, Peixiang
2016-03-01
High-order harmonic generation (HHG) from aligned molecules has provided a promising way to probe the molecular orbital with an Ångström resolution. This method, usually called molecular orbital tomography (MOT) replies on a simple assumption of the plane-wave approximation (PW), which has long been questioned due to that PW approximation is known to be valid in the keV energy region. However, the photon energy is usually no more than 100 eV in HHG. In this work, we experimentally reconstruct the highest occupied molecular orbital (HOMO) of nitrogen (N2) by using a Coulomb-corrected MOT (CCMOT) method. In our scheme, the molecular continuum states are described by a Coulomb wave function instead of the PW approximation. With CCMOT, the reconstructed orbital is demonstrated to agree well with the theoretical prediction and retain the main features of the HOMO of N2. Compared to the PW approximation method, the CCMOT shows a significant improvement in eliminating the artificial structures caused by PW approximation.
Coulomb-corrected molecular orbital tomography of nitrogen.
Zhai, Chunyang; He, Lixin; Lan, Pengfei; Zhu, Xiaosong; Li, Yang; Wang, Feng; Shi, Wenjing; Zhang, Qingbin; Lu, Peixiang
2016-01-01
High-order harmonic generation (HHG) from aligned molecules has provided a promising way to probe the molecular orbital with an Ångström resolution. This method, usually called molecular orbital tomography (MOT) replies on a simple assumption of the plane-wave approximation (PW), which has long been questioned due to that PW approximation is known to be valid in the keV energy region. However, the photon energy is usually no more than 100 eV in HHG. In this work, we experimentally reconstruct the highest occupied molecular orbital (HOMO) of nitrogen (N2) by using a Coulomb-corrected MOT (CCMOT) method. In our scheme, the molecular continuum states are described by a Coulomb wave function instead of the PW approximation. With CCMOT, the reconstructed orbital is demonstrated to agree well with the theoretical prediction and retain the main features of the HOMO of N2. Compared to the PW approximation method, the CCMOT shows a significant improvement in eliminating the artificial structures caused by PW approximation. PMID:27000666
THE FUELING DIAGRAM: LINKING GALAXY MOLECULAR-TO-ATOMIC GAS RATIOS TO INTERACTIONS AND ACCRETION
Stark, David V.; Kannappan, Sheila J.; Eckert, Kathleen D.; Wei, Lisa H.; Baker, Andrew J.; Leroy, Adam K.; Vogel, Stuart N.
2013-05-20
To assess how external factors such as local interactions and fresh gas accretion influence the global interstellar medium of galaxies, we analyze the relationship between recent enhancements of central star formation and total molecular-to-atomic (H{sub 2}/H I) gas ratios, using a broad sample of field galaxies spanning early-to-late type morphologies, stellar masses of 10{sup 7.2}-10{sup 11.2} M{sub Sun }, and diverse stages of evolution. We find that galaxies occupy several loci in a ''fueling diagram'' that plots H{sub 2}/H I ratio versus mass-corrected blue-centeredness, a metric tracing the degree to which galaxies have bluer centers than the average galaxy at their stellar mass. Spiral galaxies of all stellar masses show a positive correlation between H{sub 2}/H I ratio and mass-corrected blue-centeredness. When combined with previous results linking mass-corrected blue-centeredness to external perturbations, this correlation suggests a systematic link between local galaxy interactions and molecular gas inflow/replenishment. Intriguingly, E/S0 galaxies show a more complex picture: some follow the same correlation, some are quenched, and a distinct population of blue-sequence E/S0 galaxies (with masses below key scales associated with transitions in gas richness) defines a separate loop in the fueling diagram. This population appears to be composed of low-mass merger remnants currently in late- or post-starburst states, in which the burst first consumes the H{sub 2} while the galaxy center keeps getting bluer, then exhausts the H{sub 2}, at which point the burst population reddens as it ages. Multiple lines of evidence suggest connected evolutionary sequences in the fueling diagram. In particular, tracking total gas-to-stellar mass ratios within the fueling diagram provides evidence of fresh gas accretion onto low-mass E/S0s emerging from their central starburst episodes. Drawing on a comprehensive literature search, we suggest that virtually all galaxies
Conformation effects on the molecular orbitals of serine
NASA Astrophysics Data System (ADS)
Wang, Ke-Dong; Ma, Peng-Fei; Shan, Xu
2011-03-01
This paper calculates the five most stable conformers of serine with Hartree—Fock theory, density functional theory (B3LYP), Møller—Plesset perturbation theory (MP4(SDQ)) and electron propagation theory with the 6-311++G(2d,2p) basis set. The calculated vertical ionization energies for the valence molecular orbitals of each conformer are in agreement with the experimental data, indicating that a range of molecular conformations would coexist in an equilibrium sample. Information of the five outer valence molecular orbitals for each conformer is explored in coordinate and momentum spaces using dual space analysis to investigate the conformational processes, which are generated from the global minimum conformer Ser1 by rotation of C2-C3 (Ser4), C1-C2 (Ser5) and C1-O2 (Ser2 and Ser3). Orbitals 28a, 27a and 26a are identified as the fingerprint orbitals for all the conformational processes. Project supported by the Doctoral Research Fund of Henan Normal University, China (Grant No. 525449).
NASA Astrophysics Data System (ADS)
Abramo, M. C.; Caccamo, C.; Costa, D.; Munaò, G.
2014-09-01
We report an atomistic molecular dynamics determination of the phase diagram of a rigid-cage model of C36. We first show that free energies obtained via thermodynamic integrations along isotherms displaying "van der Waals loops," are fully reproduced by those obtained via isothermal-isochoric integration encompassing only stable states. We find that a similar result also holds for isochoric paths crossing van der Waals regions of the isotherms, and for integrations extending to rather high densities where liquid-solid coexistence can be expected to occur. On such a basis we are able to map the whole phase diagram of C36, with resulting triple point and critical temperatures about 1770 K and 2370 K, respectively. We thus predict a 600 K window of existence of a stable liquid phase. Also, at the triple point density, we find that the structural functions and the diffusion coefficient maintain a liquid-like character down to 1400-1300 K, this indicating a wide region of possible supercooling. We discuss why all these features might render possible the observation of the melting of C36 fullerite and of its liquid state, at variance with what previously experienced for C60.
Optimization of selected molecular orbitals in group basis sets.
Ferenczy, György G; Adams, William H
2009-04-01
We derive a local basis equation which may be used to determine the orbitals of a group of electrons in a system when the orbitals of that group are represented by a group basis set, i.e., not the basis set one would normally use but a subset suited to a specific electronic group. The group orbitals determined by the local basis equation minimize the energy of a system when a group basis set is used and the orbitals of other groups are frozen. In contrast, under the constraint of a group basis set, the group orbitals satisfying the Huzinaga equation do not minimize the energy. In a test of the local basis equation on HCl, the group basis set included only 12 of the 21 functions in a basis set one might ordinarily use, but the calculated active orbital energies were within 0.001 hartree of the values obtained by solving the Hartree-Fock-Roothaan (HFR) equation using all 21 basis functions. The total energy found was just 0.003 hartree higher than the HFR value. The errors with the group basis set approximation to the Huzinaga equation were larger by over two orders of magnitude. Similar results were obtained for PCl(3) with the group basis approximation. Retaining more basis functions allows an even higher accuracy as shown by the perfect reproduction of the HFR energy of HCl with 16 out of 21 basis functions in the valence basis set. When the core basis set was also truncated then no additional error was introduced in the calculations performed for HCl with various basis sets. The same calculations with fixed core orbitals taken from isolated heavy atoms added a small error of about 10(-4) hartree. This offers a practical way to calculate wave functions with predetermined fixed core and reduced base valence orbitals at reduced computational costs. The local basis equation can also be used to combine the above approximations with the assignment of local basis sets to groups of localized valence molecular orbitals and to derive a priori localized orbitals. An
Liu, Shubin E-mail: schauer@unc.edu; Schauer, Cynthia K. E-mail: schauer@unc.edu
2015-02-07
To have a quantitative understanding about the origin of conformation stability for molecular systems is still an unaccomplished task. Frontier orbital interactions from molecular orbital theory and energy partition schemes from density functional reactivity theory are the two approaches available in the literature that can be used for this purpose. In this work, we compare the performance of these approaches for a total of 48 simple molecules. We also conduct studies to flexibly bend bond angles for water, carbon dioxide, borane, and ammonia molecules to obtain energy profiles for these systems over a wide range of conformations. We find that results from molecular orbital interactions using frontier occupied orbitals such as the highest occupied molecular orbital and its neighbors are only qualitatively, at most semi-qualitatively, trustworthy. To obtain quantitative insights into relative stability of different conformations, the energy partition approach from density functional reactivity theory is much more reliable. We also find that the electrostatic interaction is the dominant descriptor for conformational stability, and steric and quantum effects are smaller in contribution but their contributions are indispensable. Stable molecular conformations prefer to have a strong electrostatic interaction, small molecular size, and large exchange-correlation effect. This work should shed new light towards establishing a general theoretical framework for molecular stability.
Intracellular molecular distributions in spacecraft experiments in orbit around Earth
NASA Astrophysics Data System (ADS)
Haranas, Ioannis; Gkigkitzis, Ioannis; Zouganelis, George D.
2012-04-01
It is possible that the nucleolous inside the cell plays the role of a "gravity receptor". Furthermore, cells up to 10 μm in diameter can demonstrate some effect due to the redistribution of mitochondria or nucleolous. Effects of gravity should be present in various cell systems where larger objects such as the ribosomes move from cell to cell. In this paper we study the effects of gravity on cells. In particular, we examine the resulting intracellular molecular distribution due to Brownian motion and the ordered distribution of molecules under the action of gravity, where n0 is the number per unit volume at certain level, and n is the number per unit volume above that level. This is an experiment that takes place at a certain orbital altitude in a spacecraft in orbit around Earth, where the acceleration due to the central field is corrected for the oblateness and also the rotation of the Earth. We found that equatorial circular and elliptical orbits have the highest n/n0 ratios. This experiment takes place in circular and elliptical orbits, with eccentricities e = 0, 0.1 and involves a bacterial cell at an orbital altitude of 300 km. We found that n/n0 = 1.00299 and 1.0037 respectively, which is still a 0.6-0.7 % higher than n/n0 = 0.0996685 calculated on the surface of the Earth. Examining mitochondria in similar orbital experiments we found that equatorial orbits result to higher n/n0 ratios. In particular, we found that n/n0 = 8.38119, where an elliptical orbit of eccentricity e = 0.1 results to n/n0 = 13.8525. Both are high above 100%, signifying the importance of Brownian motion over gravity. Our results are of interest to biomedical applications. Molecular concentrations are important for various processes such as the embryogenesis, positional homeostasis and its relation to cell energy expenditure, cell torque, cell deformation, and more. These results indicate that statistical molecular distributions play an important role for the recognition of a
A Simple Demonstration of Atomic and Molecular Orbitals Using Circular Magnets
ERIC Educational Resources Information Center
Chakraborty, Maharudra; Mukhopadhyay, Subrata; Das, Ranendu Sekhar
2014-01-01
A quite simple and inexpensive technique is described here to represent the approximate shapes of atomic orbitals and the molecular orbitals formed by them following the principles of the linear combination of atomic orbitals (LCAO) method. Molecular orbitals of a few simple molecules can also be pictorially represented. Instructors can employ the…
Super-atom molecular orbital excited states of fullerenes.
Johansson, J Olof; Bohl, Elvira; Campbell, Eleanor E B
2016-09-13
Super-atom molecular orbitals are orbitals that form diffuse hydrogenic excited electronic states of fullerenes with their electron density centred at the centre of the hollow carbon cage and a significant electron density inside the cage. This is a consequence of the high symmetry and hollow structure of the molecules and distinguishes them from typical low-lying molecular Rydberg states. This review summarizes the current experimental and theoretical studies related to these exotic excited electronic states with emphasis on femtosecond photoelectron spectroscopy experiments on gas-phase fullerenes.This article is part of the themed issue 'Fullerenes: past, present and future, celebrating the 30th anniversary of Buckminster Fullerene'. PMID:27501970
Analyzing molecular static linear response properties with perturbed localized orbitals
NASA Astrophysics Data System (ADS)
Autschbach, Jochen; King, Harry F.
2010-07-01
Perturbed localized molecular orbitals (LMOs), correct to first order in an applied static perturbation and consistent with a chosen localization functional, are calculated using analytic derivative techniques. The formalism is outlined for a general static perturbation and variational localization functionals. Iterative and (formally) single-step approaches are compared. The implementation employs an iterative sequence of 2×2 orbital rotations. The procedure is verified by calculations of molecular electric-field perturbations. Boys LMO contributions to the electronic static polarizability and the electric-field perturbation of the ⟨r2⟩ expectation value are calculated and analyzed for ethene, ethyne, and fluoroethene (H2CCHF). For ethene, a comparison is made with results from a Pipek-Mezey localization. The calculations show that a chemically intuitive decomposition of the calculated properties is possible with the help of the LMO contributions and that the polarizability contributions in similar molecules are approximately transferable.
Emulating Molecular Orbitals and Electronic Dynamics with Ultracold Atoms
NASA Astrophysics Data System (ADS)
Lühmann, Dirk-Sören; Weitenberg, Christof; Sengstock, Klaus
2015-07-01
In recent years, ultracold atoms in optical lattices have proven their great value as quantum simulators for studying strongly correlated phases and complex phenomena in solid-state systems. Here, we reveal their potential as quantum simulators for molecular physics and propose a technique to image the three-dimensional molecular orbitals with high resolution. The outstanding tunability of ultracold atoms in terms of potential and interaction offer fully adjustable model systems for gaining deep insight into the electronic structure of molecules. We study the orbitals of an artificial benzene molecule and discuss the effect of tunable interactions in its conjugated π electron system with special regard to localization and spin order. The dynamical time scales of ultracold atom simulators are on the order of milliseconds, which allows for the time-resolved monitoring of a broad range of dynamical processes. As an example, we compute the hole dynamics in the conjugated π system of the artificial benzene molecule.
GAUSSIAN 76: An ab initio Molecular Orbital Program
DOE R&D Accomplishments Database
Binkley, J. S.; Whiteside, R.; Hariharan, P. C.; Seeger, R.; Hehre, W. J.; Lathan, W. A.; Newton, M. D.; Ditchfield, R.; Pople, J. A.
1978-01-01
Gaussian 76 is a general-purpose computer program for ab initio Hartree-Fock molecular orbital calculations. It can handle basis sets involving s, p and d-type Gaussian functions. Certain standard sets (STO-3G, 4-31G, 6-31G*, etc.) are stored internally for easy use. Closed shell (RHF) or unrestricted open shell (UHF) wave functions can be obtained. Facilities are provided for geometry optimization to potential minima and for limited potential surface scans.
A low Earth orbit molecular beam space simulation facility
NASA Technical Reports Server (NTRS)
Cross, J. B.
1984-01-01
A brief synopsis of the low Earth orbit (LEO) satellite environment is presented including neutral and ionic species. Two ground based atomic and molecular beam instruments are described which are capable of simulating the interaction of spacecraft surfaces with the LEO environment and detecting the results of these interactions. The first detects mass spectrometrically low level fluxes of reactively and nonreactively surface scattered species as a function of scattering angle and velocity while the second ultrahigh velocity (UHV) molecular beam, laser induced fluorescence apparatus is capable of measuring chemiluminescence produced by either gas phase or gas-surface interactions. A number of proposed experiments are described.
Tomographic Imaging of Molecular Orbitals in Length and Velocity Form
Zwan, Elmar V. van der; Lein, Manfred
2007-11-29
Recently Itatani et al. [Nature 432, 876 (2004)] introduced the new concept of molecular orbital tomography, where high harmonic generation (HHG) is used to image electronic wave functions. We describe an alternative reconstruction form, using momentum instead of dipole matrix elements for the electron recombination step in HHG. We show that using this velocity-form reconstruction, one obtains better results than using the original length-form reconstruction. We provide numerical evidence for our claim that one has to resort to extremely short pulses to perform the reconstruction for an orbital with arbitrary symmetry. The numerical evidence is based on the exact solution of the time-dependent Schroedinger equation for 2D model systems to simulate the experiment. Furthermore we show that in the case of cylindrically symmetric orbitals, such as the N{sub 2} orbital that was reconstructed in the original work, one can obtain the full 3D wave function and not only a 2D projection of it.
Localization of molecular orbitals: from fragments to molecule.
Li, Zhendong; Li, Hongyang; Suo, Bingbing; Liu, Wenjian
2014-09-16
Conspectus Localized molecular orbitals (LMO) not only serve as an important bridge between chemical intuition and molecular wave functions but also can be employed to reduce the computational cost of many-body methods for electron correlation and excitation. Therefore, how to localize the usually completely delocalized canonical molecular orbitals (CMO) into confined physical spaces has long been an important topic: It has a long history but still remains active to date. While the known LMOs can be classified into (exact) orthonormal and nonorthogonal, as well as (approximate) absolutely localized MOs, the ways for achieving these can be classified into two categories, a posteriori top-down and a priori bottom-up, depending on whether they invoke the global CMOs (or equivalently the molecular density matrix). While the top-down approaches have to face heavy tasks of minimizing or maximizing a given localization functional typically of many adjacent local extrema, the bottom-up ones have to invoke some tedious procedures for first generating a local basis composed of well-defined occupied and unoccupied subsets and then maintaining or resuming the locality when solving the Hartree-Fock/Kohn-Sham (HF/KS) optimization condition. It is shown here that the good of these kinds of approaches can be combined together to form a very efficient hybrid approach that can generate the desired LMOs for any kind of gapped molecules. Specifically, a top-down localization functional, applied to individual small subsystems only, is minimized to generate an orthonormal local basis composed of functions centered on the preset chemical fragments. The familiar notion for atomic cores, lone pairs, and chemical bonds emerges here automatically. Such a local basis is then employed in the global HF/KS calculation, after which a least action is taken toward the final orthonormal localized molecular orbitals (LMO), both occupied and virtual. This last step is very cheap, implying that, after
Molecular orbital analysis of the hydrogen bonded water dimer
NASA Astrophysics Data System (ADS)
Wang, Bo; Jiang, Wanrun; Dai, Xin; Gao, Yang; Wang, Zhigang; Zhang, Rui-Qin
2016-02-01
As an essential interaction in nature, hydrogen bonding plays a crucial role in many material formations and biological processes, requiring deeper understanding. Here, using density functional theory and post-Hartree-Fock methods, we reveal two hydrogen bonding molecular orbitals crossing the hydrogen-bond’s O and H atoms in the water dimer. Energy decomposition analysis also shows a non-negligible contribution of the induction term. Our finding sheds light on the essential understanding of hydrogen bonding in ice, liquid water, functional materials and biological systems.
Molecular orbital analysis of the hydrogen bonded water dimer
Wang, Bo; Jiang, Wanrun; Dai, Xin; Gao, Yang; Wang, Zhigang; Zhang, Rui-Qin
2016-01-01
As an essential interaction in nature, hydrogen bonding plays a crucial role in many material formations and biological processes, requiring deeper understanding. Here, using density functional theory and post-Hartree-Fock methods, we reveal two hydrogen bonding molecular orbitals crossing the hydrogen-bond’s O and H atoms in the water dimer. Energy decomposition analysis also shows a non-negligible contribution of the induction term. Our finding sheds light on the essential understanding of hydrogen bonding in ice, liquid water, functional materials and biological systems. PMID:26905305
Kafafi, S.A. ); LaFemina, J.P. ); Nauss, J.L. )
1990-11-21
Full geometry optimizations using molecular mechanics and the quantum chemical AM1 method have been carried out to determine the minimum energy conformation of pyromellitic dianhydride-oxydianiline polyimide (PMDA-ODA PI). The phenyl-imide twist angle for this compound was determined to be {approximately}30. These computations also provided a quantitative determination of the energy gap (7 eV), electron affinity ({minus}2 eV), and ionization potential (8.97 eV). Computations on the PMDA-ODA PI radical anion provided an estimate of the hopping barrier for an electron to hop from one chain to another (3.2 eV), the mechanism believed responsible for photoconduction. Moreover, the use of qualitative molecular orbital theory (QMOT) arguments provided an interpretation of these results in a simple molecular orbital framework.
da Silva, Rodrigo R; Ramalho, Teodorico C; Santos, Joana M; Figueroa-Villar, J Daniel
2006-01-26
We carried out Hartree-Fock (HF) and density functional theory calculations for 61 compounds, the conjugated bases of carboxylic acids, phenols, and alcohols, and analyzed their acid-base behavior using molecular orbital (MO) energies and their dependence on solvent effects. Despite the well-known correlation between highest-occupied MO (HOMO) energies and pKa, we observed that HOMO energies are inadequate to describe the acid-base behavior of these compounds. Therefore, we established a criterion to identify the best frontier MO for describing pKa values and also to understand why the HOMO approach fails. The MO that fits our criterion provided very good correlations with pKa values, much better than those obtained by HOMO energies. Since they are the frontier molecular orbitals that drive the acid-base reactions in each compound, they were called frontier effective-for-reaction MOs, or FERMOs. By use of the FERMO concept, the reactions that are HOMO driven, and those that are not, can be better explained, independently from the calculation method used, as both HF and Kohn-Sham methodologies lead to the same FERMO. PMID:16420004
The activation strain model and molecular orbital theory
Wolters, Lando P; Bickelhaupt, F Matthias
2015-01-01
The activation strain model is a powerful tool for understanding reactivity, or inertness, of molecular species. This is done by relating the relative energy of a molecular complex along the reaction energy profile to the structural rigidity of the reactants and the strength of their mutual interactions: ΔE(ζ) = ΔEstrain(ζ) + ΔEint(ζ). We provide a detailed discussion of the model, and elaborate on its strong connection with molecular orbital theory. Using these approaches, a causal relationship is revealed between the properties of the reactants and their reactivity, e.g., reaction barriers and plausible reaction mechanisms. This methodology may reveal intriguing parallels between completely different types of chemical transformations. Thus, the activation strain model constitutes a unifying framework that furthers the development of cross-disciplinary concepts throughout various fields of chemistry. We illustrate the activation strain model in action with selected examples from literature. These examples demonstrate how the methodology is applied to different research questions, how results are interpreted, and how insights into one chemical phenomenon can lead to an improved understanding of another, seemingly completely different chemical process. WIREs Comput Mol Sci 2015, 5:324–343. doi: 10.1002/wcms.1221 PMID:26753009
Charge transfer processes: the role of optimized molecular orbitals.
Meyer, Benjamin; Domingo, Alex; Krah, Tim; Robert, Vincent
2014-08-01
The influence of the molecular orbitals on charge transfer (CT) reactions is analyzed through wave function-based calculations. Characteristic CT processes in the organic radical 2,5-di-tert-butyl-6-oxophenalenoxyl linked with tetrathiafulvalene and the inorganic crystalline material LaMnO3 show that changes in the inner shells must be explicitly taken into account. Such electronic reorganization can lead to a reduction of the CT vertical transition energy up to 66%. A state-specific approach accessible through an adapted CASSCF (complete active space self-consistent field) methodology is capable of reaching good agreement with the experimental spectroscopy of CT processes. A partitioning of the relaxation energy in terms of valence- and inner-shells is offered and sheds light on their relative importance. This work paves the way to the intimate description of redox reactions using quantum chemistry methods. PMID:24781811
P and N compensation in diamond molecular orbital theory
NASA Astrophysics Data System (ADS)
Anderson, Alfred B.; Kostadinov, Lubomir N.
1997-01-01
Cluster models and the atom superposition and electron delocalization molecular orbital theory calculations lead to an explanation for the ability of nitrogen to cause phosphorous incorporation in low pressure grown diamond films as observed recently by Cao and coworkers. The theory shows that substitutional N compensates substitutional P, creating stable P+-N- disubstitutional pairs. These ionized systems are calculated to be deep donors, which explains the absence of measurable electrical conductivity or phosphorous induced luminescence. The possibility of creating donor P defects by the annealing reaction P-N+N→P+N-N is discussed. The issues of atom size and electronegativity and their influence on donor capability are addressed. It is shown that the difference between substitutional P, a shallow donor, and substitutional N, a deep donor, is predominantly due to the larger size of P; its lower electronegativity makes a relatively small contribution.
Molecular electric moments calculated by using natural orbital functional theory.
Mitxelena, Ion; Piris, Mario
2016-05-28
The molecular electric dipole, quadrupole, and octupole moments of a selected set of 21 spin-compensated molecules are determined employing the extended version of the Piris natural orbital functional 6 (PNOF6), using the triple-ζ Gaussian basis set with polarization functions developed by Sadlej, at the experimental geometries. The performance of the PNOF6 is established by carrying out a statistical analysis of the mean absolute errors with respect to the experiment. The calculated PNOF6 electric moments agree satisfactorily with the corresponding experimental data and are in good agreement with the values obtained by accurate ab initio methods, namely, the coupled-cluster single and doubles and multi-reference single and double excitation configuration interaction methods. PMID:27250280
Theoretical analysis of the density within an orbiting molecular shield
NASA Technical Reports Server (NTRS)
Hueser, J. E.; Brock, F. J.
1976-01-01
An analytical model based on the kinetic theory of a drifting Maxwellian gas is used to determine the nonequilibrium molecular density distribution within a hemispherical shell open aft with its axis parallel to its velocity. Separate numerical results are presented for the primary and secondary density distribution components due to the drifting Maxwellian gas for speed ratios between 2.5 and 10. An analysis is also made of the density component due to gas desorbed from the wall of the hemisphere, and numerical results are presented for the density distribution. It is shown that the adsorption process may be completely ignored. The results are applicable to orbital trajectories in any planet-atmosphere system and interplanetary transfer trajectories. Application to the earth's atmosphere is mentioned briefly.
Sulfur at nickel-alumina interfaces - Molecular orbital theory
NASA Technical Reports Server (NTRS)
Hong, S. Y.; Anderson, Alfred B.; Smialek, James L.
1990-01-01
Previous studies on Al-Ni alloys containing sulfur as an impurity suggest that, when S is in the interface between a metal and an oxide scale, it weakens the chemical bonding between them. This paper investigates factors responsible for this effect, using a molecular orbital theory to predict sulfur structures and electronic properties on the Ni-Al2O3 interface. It is shown that, in absence of S, the basal plane of Al2O3 will bind strongly through the Al(3+) cation surface to Ni (111). When segregated S impurity is present on the Ni surface, there are too few interfacial AlS bonds to effect good adhesion, leading to an inhibition of the oxide scale adhesion in NiCrAl alloys.
Molecular electric moments calculated by using natural orbital functional theory
NASA Astrophysics Data System (ADS)
Mitxelena, Ion; Piris, Mario
2016-05-01
The molecular electric dipole, quadrupole, and octupole moments of a selected set of 21 spin-compensated molecules are determined employing the extended version of the Piris natural orbital functional 6 (PNOF6), using the triple-ζ Gaussian basis set with polarization functions developed by Sadlej, at the experimental geometries. The performance of the PNOF6 is established by carrying out a statistical analysis of the mean absolute errors with respect to the experiment. The calculated PNOF6 electric moments agree satisfactorily with the corresponding experimental data and are in good agreement with the values obtained by accurate ab initio methods, namely, the coupled-cluster single and doubles and multi-reference single and double excitation configuration interaction methods.
Band Formation in a Molecular Quantum Well via 2D Superatom Orbital Interactions
Dougherty, D. B.; Feng, Min; Petek, Hrvoje; Yates, John T.; Zhao, Jin
2012-12-28
By scanning tunneling microscopy and spectroscopy, we study nearly free electron band formation of the σ*lowest unoccupied molecular orbital of C₆F₆ on a Cu(111) surface. In fractal islands, the lowest unoccupied molecular orbital energy systematically stabilizes with the number of interacting near-neighbor C₆F₆ molecules. Density functional theory calculations reveal the origin of effective intermo- lecular orbital overlap in the previously unrecognized superatom character of the σ*orbital of ₆F₆ molecules. The discovery of superatom orbitals in planar molecules offers a new universal principle for effective band formation, which can be exploited in designing organic semiconductors with nearly free electron properties
Analyzing and Interpreting NMR Spin-Spin Coupling Constants Using Molecular Orbital Calculations
ERIC Educational Resources Information Center
Autschbach, Jochen; Le Guennic, Boris
2007-01-01
Molecular orbital plots are used to analyze and interpret NMR spin-spin coupling constants, also known as J coupling constants. Students have accepted the concept of contributions to molecular properties from individual orbitals without the requirement to provide explicit equations.
Offenbacher, Hannes; Lüftner, Daniel; Ules, Thomas; Reinisch, Eva Maria; Koller, Georg; Puschnig, Peter; Ramsey, Michael G.
2015-01-01
The frontier orbitals of molecules are the prime determinants of their chemical, optical and electronic properties. Arguably, the most direct method of addressing the (filled) frontier orbitals is ultra-violet photoemission spectroscopy (UPS). Although UPS is a mature technique from the early 1970s on, the angular distribution of the photoemitted electrons was thought to be too complex to be analysed quantitatively. Recently angle resolved UPS (ARUPS) work on conjugated molecules both, in ordered thick films and chemisorbed monolayers, has shown that the angular (momentum) distribution of the photocurrent from orbital emissions can be simply understood. The approach, based on the assumption of a plane wave final state is becoming known as orbital tomography. Here we will demonstrate, with selected examples of pentacene (5A) and sexiphenyl (6P), the potential of orbital tomography. First it will be shown how the full angular distribution of the photocurrent (momentum map) from a specific orbital is related to the real space orbital by a Fourier transform. Examples of the reconstruction of 5A orbitals will be given and the procedure for recovering the lost phase information will be outlined. We then move to examples of sexiphenyl where we interrogate the original band maps of thick sexiphenyl in the light of our understanding of orbital tomography that has developed since then. With comparison to theoretical simulations of the molecular band maps, the molecular conformation and orientation will be concluded. New results for the sexiphenyl monolayer on Al(1 1 0) will then be presented. From the band maps it will be concluded that the molecule is planarised and adopts a tilted geometry. Finally the momentum maps down to HOMO-11 will be analysed and real space orbitals reconstructed. PMID:26752804
Nakata, Hiroya; Schmidt, Michael W; Fedorov, Dmitri G; Kitaura, Kazuo; Nakamura, Shinichiro; Gordon, Mark S
2014-10-16
The fully analytic energy gradient has been developed and implemented for the restricted open-shell Hartree–Fock (ROHF) method based on the fragment molecular orbital (FMO) theory for systems that have multiple open-shell molecules. The accuracy of the analytic ROHF energy gradient is compared with the corresponding numerical gradient, illustrating the accuracy of the analytic gradient. The ROHF analytic gradient is used to perform molecular dynamics simulations of an unusual open-shell system, liquid oxygen, and mixtures of oxygen and nitrogen. These molecular dynamics simulations provide some insight about how triplet oxygen molecules interact with each other. Timings reveal that the method can calculate the energy gradient for a system containing 4000 atoms in only 6 h. Therefore, it is concluded that the FMO-ROHF method will be useful for investigating systems with multiple open shells.
Amide N-oxides: an ab initio molecular orbital study
NASA Astrophysics Data System (ADS)
Greenberg, Arthur; DuBois, Thomas D.
2001-06-01
There are no known examples of amide N-oxides. The present study employs ab initio molecular orbital calculations at the 6-3G ∗ level to explore potential target molecules in this class. Bridgehead bicyclic lactams appear to be attractive targets for oxidation to form the corresponding N-oxides because they have reduced (or zero) amide resonance energy. The amide N-oxide linkage is predicted to have a ca. 9-10 kcal/mol rotational barrier due to eclipsing of nonbonded oxygen atoms in the transition state. The linkage has a nearly flat conformational ( ΦON-CO) profile in the range 120-240° and this suggests that a very sterically hindered acyclic amide N-oxide may be a practical synthetic target. The oxidation of strained amides is calculated to be highly exothermic if dimethyldioxirane is employed. This reagent is predicted to react appreciably exothermically with normal, stable amides such as N, N-dimethylacetamide, thus offering the potential for generating and studying such relatively unstable amide N-oxides at low temperatures.
Towards simple orbital-dependent density functionals for molecular dissociation
NASA Astrophysics Data System (ADS)
Zhang, Igor Ying; Richter, Patrick; Scheffler, Matthias
2015-03-01
Density functional theory (DFT) is one of the leading first-principles electronic-structure theories. However, molecular dissociation remains a challenge, because it requires a well-balanced description of the drastically different electronic structure at different bond lengths. One typical and well-documented case is the dissociation of both H2+ and H2, for which all popular DFT functionals fail. We start from the Bethe-Goldstone equation to propose a simple orbital-dependent correlation functional which generalizes the linear adiabatic connection approach. The resulting scheme is based on second-order perturbation theory (PT2), but includes the self-consistent coupling of electron-hole pairs, which ensures the correct H2 dissociation limit and gives a finite correlation energy for systems with a (near)-degenerate energy gap. This coupling PT2-like (CPT2) approximation delivers a significant improvement over all existing functionals for both H2 and H2+ dissociation. We will demonstrate the reason for this improvement analytically for H2 in a minimal basis.
[Applications of the Fragment Molecular Orbital Method in Drug Discovery].
Ishikawa, Takeshi
2016-01-01
Recently, ab initio quantum mechanical calculations have been applied to large molecules, including biomolecular systems. The fragment molecular orbital (FMO) method is one of the most efficient approaches for the quantum mechanical investigation of such molecules. In the FMO method, dividing a target molecule into small fragments reduces computational effort. The clear definition of inter-fragment interaction energy (IFIE) as an expression of total energy is another valuable feature of the FMO method because it provides the ability to analyze interactions in biomolecules. Thus, the FMO method is expected to be useful for drug discovery. This study demonstrates applications of the FMO method related to drug discovery. First, IFIE, according to FMO calculations, was used in the optimization of drug candidates for the development of anti-prion compounds. The second example involved interaction analysis of the human immunodeficiency virus type 1 (HIV-1) protease and a drug compound that used a novel analytical method for dispersion interaction, i.e., fragment interaction analysis based on LMP2 (FILM). PMID:26725679
ERIC Educational Resources Information Center
Johnson, Adam R.
2013-01-01
A molecular orbital (MO) diagram, especially its frontier orbitals, explains the bonding and reactivity for a chemical compound. It is therefore important for students to learn how to construct one. The traditional methods used to derive these diagrams rely on linear algebra techniques to combine ligand orbitals into symmetry-adapted linear…
Gilbert, Kathleen M; Skawinski, William J; Misra, Milind; Paris, Kristina A; Naik, Neelam H; Buono, Ronald A; Deutsch, Howard M; Venanzi, Carol A
2004-11-01
Methylphenidate (MP) binds to the cocaine binding site on the dopamine transporter and inhibits reuptake of dopamine, but does not appear to have the same abuse potential as cocaine. This study, part of a comprehensive effort to identify a drug treatment for cocaine abuse, investigates the effect of choice of calculation technique and of solvent model on the conformational potential energy surface (PES) of MP and a rigid methylphenidate (RMP) analogue which exhibits the same dopamine transporter binding affinity as MP. Conformational analysis was carried out by the AM1 and AM1/SM5.4 semiempirical molecular orbital methods, a molecular mechanics method (Tripos force field with the dielectric set equal to that of vacuum or water) and the HF/6-31G* molecular orbital method in vacuum phase. Although all three methods differ somewhat in the local details of the PES, the general trends are the same for neutral and protonated MP. In vacuum phase, protonation has a distinctive effect in decreasing the regions of space available to the local conformational minima. Solvent has little effect on the PES of the neutral molecule and tends to stabilize the protonated species. The random search (RS) conformational analysis technique using the Tripos force field was found to be capable of locating the minima found by the molecular orbital methods using systematic grid search. This suggests that the RS/Tripos force field/vacuum phase protocol is a reasonable choice for locating the local minima of MP. However, the Tripos force field gave significantly larger phenyl ring rotational barriers than the molecular orbital methods for MP and RMP. For both the neutral and protonated cases, all three methods found the phenyl ring rotational barriers for the RMP conformers/invertamers (denoted as cte, tte, and cta) to be: cte, tte > MP > cta. Solvation has negligible effect on the phenyl ring rotational barrier of RMP. The B3LYP/6-31G* density functional method was used to calculate the
Quantum Monte Carlo study of the phase diagram of solid molecular hydrogen at extreme pressures.
Drummond, N D; Monserrat, Bartomeu; Lloyd-Williams, Jonathan H; López Ríos, P; Pickard, Chris J; Needs, R J
2015-01-01
Establishing the phase diagram of hydrogen is a major challenge for experimental and theoretical physics. Experiment alone cannot establish the atomic structure of solid hydrogen at high pressure, because hydrogen scatters X-rays only weakly. Instead, our understanding of the atomic structure is largely based on density functional theory (DFT). By comparing Raman spectra for low-energy structures found in DFT searches with experimental spectra, candidate atomic structures have been identified for each experimentally observed phase. Unfortunately, DFT predicts a metallic structure to be energetically favoured at a broad range of pressures up to 400 GPa, where it is known experimentally that hydrogen is non-metallic. Here we show that more advanced theoretical methods (diffusion quantum Monte Carlo calculations) find the metallic structure to be uncompetitive, and predict a phase diagram in reasonable agreement with experiment. This greatly strengthens the claim that the candidate atomic structures accurately model the experimentally observed phases. PMID:26215251
Dependence of the conductance change on the molecular orbitals in Ag and Au electrodes
NASA Astrophysics Data System (ADS)
Kim, Hyun Sung; Kim, Taekyeong
2016-01-01
The conductance change in single-molecule junctions due to the molecular orbitals and the metal's Fermi energy was investigated by using a scanning tunneling microscopy break-junction technique with Ag and Au electrodes. 4,4'-diaminobiphenyl and 4,4'-dicyanobiphenyl as the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) junctions were used. The amine-linked HOMO-conducting junction had a lower conductance, and cyano-linked LUMO-conducting junction had a higher conductance for Ag electrode compared to Au electrode. These results are attributed to the difference between the metal electrode Fermi energy and the molecular orbital level in the metal-molecule junction. Furthermore, 2,7-diaminofluorene exhibited a higher conductance but the identical molecular plateau length for the Ag electrodes compared to that of 4,4'-diaminobiphenyl indicating that the twist angle of the molecular backbone affects the conductance.
Molecular structures of carotenoids as predicted by MNDO-AM1 molecular orbital calculations
NASA Astrophysics Data System (ADS)
Hashimoto, Hideki; Yoda, Takeshi; Kobayashi, Takayoshi; Young, Andrew J.
2002-02-01
Semi-empirical molecular orbital calculations using AM1 Hamiltonian (MNDO-AM1 method) were performed for a number of biologically important carotenoid molecules, namely all- trans-β-carotene, all- trans-zeaxanthin, and all- trans-violaxanthin (found in higher plants and algae) together with all- trans-canthaxanthin, all- trans-astaxanthin, and all- trans-tunaxanthin in order to predict their stable structures. The molecular structures of all- trans-β-carotene, all- trans-canthaxanthin, and all- trans-astaxanthin predicted based on molecular orbital calculations were compared with those determined by X-ray crystallography. Predicted bond lengths, bond angles, and dihedral angles showed an excellent agreement with those determined experimentally, a fact that validated the present theoretical calculations. Comparison of the bond lengths, bond angles and dihedral angles of the most stable conformer among all the carotenoid molecules showed that the displacements are localized around the substituent groups and hence around the cyclohexene rings. The most stable conformers of all- trans-zeaxanthin and all- trans-violaxanthin gave rise to a torsion angle around the C6-C7 bond to be ±48.7 and -84.8°, respectively. This difference is a key factor in relation to the biological function of these two carotenoids in plants and algae (the xanthophyll cycle). Further analyses by calculating the atomic charges and using enpartment calculations (division of bond energies between component atoms) were performed to ascribe the cause of the different observed torsion angles.
Hybridization of organic molecular orbitals with substrate states at interfaces: PTCDA on silver.
Ziroff, J; Forster, F; Schöll, A; Puschnig, P; Reinert, F
2010-06-11
We demonstrate the application of orbital k-space tomography for the analysis of the bonding occurring at metal-organic interfaces. Using angle-resolved photoelectron spectroscopy, we probe the spatial structure of the highest occupied molecular orbital and the former lowest unoccupied molecular orbital (LUMO) of one monolayer 3, 4, 9, 10-perylene-tetracarboxylic-dianhydride (PTCDA) on Ag(110) and (111) surfaces and, in particular, the influence of the hybridization between the orbitals and the electronic states of the substrate. We are able to quantify and localize the substrate contribution to the LUMO and thus prove the metal-molecule hybrid character of this complex state. PMID:20867234
Moving Beyond the Single Center--Ways to Reinforce Molecular Orbital Theory in an Inorganic Course
ERIC Educational Resources Information Center
Cass, Marion E.; Hollingsworth, William E.
2004-01-01
It is suggested that molecular theory should be taught earlier in the inorganic chemistry curriculum even in the introductory chemistry course in order to integrate molecular orbital arguments more effectively throughout the curriculum. The method of teaching relies on having access to molecular modeling software as having access to such software…
Abramo, M. C.; Caccamo, C. Costa, D.; Munaò, G.
2014-09-07
We report an atomistic molecular dynamics determination of the phase diagram of a rigid-cage model of C{sub 36}. We first show that free energies obtained via thermodynamic integrations along isotherms displaying “van der Waals loops,” are fully reproduced by those obtained via isothermal-isochoric integration encompassing only stable states. We find that a similar result also holds for isochoric paths crossing van der Waals regions of the isotherms, and for integrations extending to rather high densities where liquid-solid coexistence can be expected to occur. On such a basis we are able to map the whole phase diagram of C{sub 36}, with resulting triple point and critical temperatures about 1770 K and 2370 K, respectively. We thus predict a 600 K window of existence of a stable liquid phase. Also, at the triple point density, we find that the structural functions and the diffusion coefficient maintain a liquid-like character down to 1400–1300 K, this indicating a wide region of possible supercooling. We discuss why all these features might render possible the observation of the melting of C{sub 36} fullerite and of its liquid state, at variance with what previously experienced for C{sub 60}.
A growth diagram for plasma-assisted molecular beam epitaxy of In-face InN
Gallinat, C. S.; Koblmueller, G.; Brown, J. S.; Speck, J. S.
2007-09-15
We investigated the role of temperature and In/N flux ratios to determine suitable growth windows for the plasma-assisted molecular beam epitaxy of In-face (0001) InN. Under vacuum, InN starts decomposing at 435 deg. C as defined by the release of N{sub 2} from the InN crystal and a buildup of an In adlayer and liquid In droplets on the sample surface. At temperatures greater than 470 deg. C, InN decomposition was characterized by a release of both In vapor and N{sub 2} in the absence of a significant accumulation of an In adlayer. No growth was observed at substrate temperatures above 500 deg. C or at temperatures in which the decomposition rates were higher than the growth rates. A growth diagram was then constructed consisting of two growth regimes: the 'In-droplet regime' characterized by step-flow growth and relatively flat surfaces and the ''N-rich regime'' characterized by rough, three-dimensional surfaces. The growth diagram can then be used to predict the surface structure of films grown at varying substrate temperatures and In fluxes. A 2.5 monolayer In adlayer was observed during In-droplet growth, suggesting that an In wetting layer was necessary for step-flow growth.
NASA Astrophysics Data System (ADS)
Bâldea, Ioan
2014-11-01
Estimating the relative alignment between the frontier molecular orbitals (MOs) that dominates the charge transport through single-molecule junctions represents a challenge for theory. This requires approaches beyond the widely employed framework provided by the density functional theory, wherein the Kohn-Sham ‘orbitals’ are treated as if they were real MOs, which is not the case. In this paper, we report results obtained by means of quantum chemical calculations, including the equation-of-motion coupled-cluster singles and doubles, which is the state-of-the-art of quantum chemistry for medium-size molecules like those considered here. These theoretical results are validated against data on the MO energy offset relative to the electrodes’ Fermi energy extracted from experiments for junctions based on 4,4’-bipyridine and 1,4-dicyanobenzene.
MOLE: a Voronoi diagram-based explorer of molecular channels, pores, and tunnels.
Petrek, Martin; Kosinová, Pavlína; Koca, Jaroslav; Otyepka, Michal
2007-11-01
We have developed an algorithm, "MOLE," for the rapid, fully automated location and characterization of molecular channels, tunnels, and pores. This algorithm has been made freely available on the Internet (http://mole.chemi.muni.cz/) and overcomes many of the shortcomings and limitations of the recently developed CAVER software. The core of our MOLE algorithm is a Dijkstra's path search algorithm, which is applied to a Voronoi mesh. Tests on a wide variety of biomolecular systems including gramicidine, acetylcholinesterase, cytochromes P450, potassium channels, DNA quadruplexes, ribozymes, and the large ribosomal subunit have demonstrated that the MOLE algorithm performs well. MOLE is thus a powerful tool for exploring large molecular channels, complex networks of channels, and molecular dynamics trajectories in which analysis of a large number of snapshots is required. PMID:17997961
Accurate description of phase diagram of clathrate hydrates at the molecular level
NASA Astrophysics Data System (ADS)
Belosludov, Rodion V.; Subbotin, Oleg S.; Mizuseki, Hiroshi; Kawazoe, Yoshiyuki; Belosludov, Vladimir R.
2009-12-01
In order to accurately estimate the thermodynamic properties of hydrogen clathrate hydrates, we developed a method based on the solid solution theory of van der Waals and Platteeuw. This model allows one to take into account the influence of guest molecules on the host lattice and guest-guest interactions—especially when more than one guest molecule occupies a cage. The free energies, equations of state, and chemical potentials of hydrogen and mixed propane-hydrogen clathrate hydrates of cubic structure II with different cage fillings have been estimated using this approach. Moreover, the proposed theory has been used for construction p -T phase diagrams of hydrogen hydrate and mixed hydrogen-propane hydrates in a wide range of pressures and temperatures. For the systems with well defined interactions the calculated curves of "guest gas-hydrate-ice Ih" equilibrium agree with the available experimental data. We also believe that the present model allows one not only to calculate the hydrogen storage ability of known hydrogen hydrate but also predict this value for structures that have not yet been realized by experiment.
Polarization and molecular-orbital dependence of strong-field enhanced ionization
NASA Astrophysics Data System (ADS)
Lai, Wei; Guo, Chunlei
2016-04-01
In this work we perform a polarization dependence study of enhanced ionization (EI) in diatomic molecules. We find that EI exists when the field polarization is parallel to the molecular axis but disappears when polarization is perpendicular. We further study EI with circular polarization and find that EI exists with circular polarization indicating that rescattering does not play a significant role for EI. Furthermore, we study molecular orbital effect on EI. We find that EI exists in σ type but not π type outmost molecular orbitals.
Orbital symmetry and interference effects in molecular high-order harmonic generation
Lagmago Kamta, G.; Bandrauk, A. D.
2009-10-15
We investigate harmonic generation from H{sub 2}{sup +} molecules driven by intense few-cycle laser pulses whose linearly polarization axis makes an arbitrary angle {chi} with respect to the molecular axis. The H{sub 2}{sup +} molecule is considered initially in various orbitals with nodal planes. It is found that a strong enhancement of high-order harmonics (HOHs) occurs when the laser polarization axis overlaps with major axes of electron distribution in the active orbital, while broad suppression of HOHs occurs when the laser polarization axis is parallel to a nodal plane of the active molecular orbital. We show that this harmonic suppression is enhanced by destructive interferences when the nodal and the laser polarization axes are parallel to the internuclear axis, which leads to a shortening of the harmonic cutoff. It follows that the orientation dependence of HOHs intensities mimics the electronic density in active orbitals through the angular dependence of ionization and recombination processes.
Molecular orbital model of optical centers in bismuth-doped glasses.
Kustov, E F; Bulatov, L I; Dvoyrin, V V; Mashinsky, V M
2009-05-15
Spectroscopic properties of optical fibers with a bismuth-doped silicate glass core are explained on the basis of molecular orbital theory and a solution of the Schrödinger equation, which takes into account the exchange, the spin-orbital, and the glass field potential interactions of s, p, and d electron shells of bismuth with s(sigma), p(sigma), and p(pi) orbits of oxygen atoms. The approach can explain the IR luminescence properties of other optical centers formed by other atoms with the same structure of electron shells as the bismuth atom. The model of transitions based on intramolecular charge transfer between molecular orbital and metallic states is proposed. PMID:19448817
Norris, Scott A; Samela, Juha; Bukonte, Laura; Backman, Marie; Djurabekova, Flyura; Nordlund, Kai; Madi, Charbel S; Brenner, Michael P; Aziz, Michael J
2011-01-01
Energetic particle irradiation can cause surface ultra-smoothening, self-organized nanoscale pattern formation or degradation of the structural integrity of nuclear reactor components. A fundamental understanding of the mechanisms governing the selection among these outcomes has been elusive. Here we predict the mechanism governing the transition from pattern formation to flatness using only parameter-free molecular dynamics simulations of single-ion impacts as input into a multiscale analysis, obtaining good agreement with experiment. Our results overturn the paradigm attributing these phenomena to the removal of target atoms via sputter erosion: the mechanism dominating both stability and instability is the impact-induced redistribution of target atoms that are not sputtered away, with erosive effects being essentially irrelevant. We discuss the potential implications for the formation of a mysterious nanoscale topography, leading to surface degradation, of tungsten plasma-facing fusion reactor walls. Consideration of impact-induced redistribution processes may lead to a new design criterion for stability under irradiation. PMID:21505432
Liu, Xiaojing; Sun, Xin; Xu, Xijin; Sun, Ping
2015-11-01
By using atomic force microscopy, birth and spread has proved to be the primary growth mechanism for L-prolinium phosphate (LPP). The phenomenon of newly formed islands expanding to the edge of the preceding terrace was observed. The optimized molecular structure and the molecular properties were calculated by density functional theory method. Natural bond orbital analysis was carried out to demonstrate the various inter and intramolecular interactions that are responsible for the stabilization of LPP leading to high NLO activity. Molecular electrostatic potential, frontier molecular orbital analysis and thermodynamic properties were investigated to get a better insight of the molecular properties. Global and local reactivity descriptors were computed to predict the reactivity and reactive sites on the molecules. Non-linear optical (NLO) properties such as the total dipole moment (μ) and first order hyperopolarizability (β) were also calculated to predict NLO behavior. PMID:26067937
Solution of multi-center molecular integrals of Slater-type orbitals
NASA Technical Reports Server (NTRS)
Tai, H.
1989-01-01
The troublesome multi-center molecular integrals of Slater-type orbitals (STO) in molecular physics calculations can be evaluated by using the Fourier transform and proper coupling of the two center exchange integrals. A numerical integration procedure is then readily rendered to the final expression in which the integrand consists of well known special functions of arguments containing the geometrical arrangement of the nuclear centers and the exponents of the atomic orbitals. A practical procedure was devised for the calculation of a general multi-center molecular integrals coupling arbitrary Slater-type orbitals. Symmetry relations and asymptotic conditions are discussed. Explicit expressions of three-center one-electron nuclear-attraction integrals and four-center two-electron repulsion integrals for STO of principal quantum number n=2 are listed. A few numerical results are given for the purpose of comparison.
Intramolecular charge ordering in the multi molecular orbital system (TTM-TTP)I3
NASA Astrophysics Data System (ADS)
Bonnet, Marie-Laure; Robert, Vincent; Tsuchiizu, Masahisa; Omori, Yukiko; Suzumura, Yoshikazu
2010-06-01
Starting from the structure of the (TTM-TTP)I3 molecular-based material, we examine the characteristics of frontier molecular orbitals using ab initio (CASSCF/CASPT2) configurations interaction calculations. It is shown that the singly occupied and second-highest-occupied molecular orbitals are close to each other, i.e., this compound should be regarded as a two-orbital system. By dividing virtually the [TTM-TTP] molecule into three fragments, an effective model is constructed to rationalize the origin of this picture. In order to investigate the low-temperature, symmetry breaking experimentally observed in the crystal, the electronic distribution in a pair of [TTM-TTP] molecules is analyzed from CASPT2 calculations. Our inspection supports and explains the speculated intramolecular charge ordering which is likely to give rise to low-energy magnetic properties.
Fragmented Molecular Orbital with Diffusion Monte Carlo for large molecular systems
NASA Astrophysics Data System (ADS)
Benali, Anouar; Pruitt, Spencer R.; Fedorov, Dmitri G.
Performing accurate quantum mechanics (QM) calculations on larger and larger systems, while maintaining a high level of accuracy is an ongoing effort in many ab initio fields. Many different attempts have been made to develop highly scalable and accurate methods. The fragment molecular orbital (FMO) method is an ab initio method capable of taking advantage of modern supercomputers, such as the Blue Gene Q system Mira at the Argonne National Laboratory Leadership Computing Facility (ALCF). FMO is based on dividing molecules into fragments and performing ab initio calculations on fragments, their dimers and, optionally, trimers. This decomposition makes it possible to perform QM calculations of real size biological molecules. In contrast to many other fragment-based methods, the effect of the environment is rigorously accounted for by computing the electrostatic potential (ESP) due to remaining fragments that are not explicitly included in a given monomer, dimer, or trimer calculation. The use of highly accurate levels of theory, such as Diffusion Monte Carlo (DMC-QMC), in conjunction with FMO allows for the goal of highly scalable and accurate all electron calculations demonstrated in this study, on a variety of relevant systems (H2O)[3-6] and protein using GAMESS and QMCPACK.
Terazulene Isomers: Polarity Change of OFETs through Molecular Orbital Distribution Contrast.
Yamaguchi, Yuji; Takubo, Maki; Ogawa, Keisuke; Nakayama, Ken-Ichi; Koganezawa, Tomoyuki; Katagiri, Hiroshi
2016-09-01
Intermolecular orbital coupling is fundamentally important to organic semiconductor performance. Recently, we reported that 2,6':2',6″-terazulene (TAz1) exhibited excellent performance as an n-type organic field-effect transistor (OFET) via molecular orbital distribution control. To validate and develop this concept, here we present three other terazulene regioisomers, which have three azulene molecules connected at the 2- or 6-position along the long axis of the azulene, thus constructing a linear expanded π-conjugation system: 2,2':6',2″-terazulene (TAz2), 2,2':6',6″-terazulene (TAz3), and 6,2':6',6″-terazulene (TAz4). TAz2 and TAz3 exhibit ambipolar characteristics; TAz4 exhibits clear n-type transistor behavior as an OFET. The lowest unoccupied molecular orbitals (LUMOs) of all terazulenes are fully delocalized over the entire molecule. In contrast, the highest occupied molecular orbitals (HOMOs) of TAz2 and TAz3 are delocalized over the 2,2'-biazulene units; the HOMOs of TAz4 are localized at one end of the azulene unit. These findings confirm that terazulene isomers which are simple hydrocarbon compounds are versatile materials with a tunable-polarity FET characteristic that depends on the direction of the azulene unit and the related contrast of the molecular orbital distribution in the terazulene backbone. PMID:27511286
Molecular orbital ordering in titania and the associated semiconducting behavior
Park, Joseph; Ok, Kyung-Chul; Park, Jin-Seong; Du Ahn, Byung; Hun Lee, Je; Park, Jae-Woo; Chung, Kwun-Bum
2011-10-03
RF-sputtered TiO{sub x} layers were thermally treated and the associated thin-film transistor properties were studied. X-ray diffraction and x-ray absorption spectroscopy analyses indicate that as-grown amorphous TiO{sub x} films crystallize to anatase at temperatures above 450 deg. C in air. Thin-film transistors incorporating anatase active layers exhibit n-type behavior, with field effect mobility values near 0.11 cm{sup 2}/Vs when annealed at 550 deg. C. Such a phenomenon is suggested to originate from the ordering of Ti 3d orbitals upon crystallization, and the mobility enhancement at higher annealing temperatures may be attributed to the reduced grain boundary scattering of carriers by virtue of enlarged average grain size.
Pruitt, Spencer R; Nakata, Hiroya; Nagata, Takeshi; Mayes, Maricris; Alexeev, Yuri; Fletcher, Graham; Fedorov, Dmitri G; Kitaura, Kazuo; Gordon, Mark S
2016-04-12
The analytic first derivative with respect to nuclear coordinates is formulated and implemented in the framework of the three-body fragment molecular orbital (FMO) method. The gradient has been derived and implemented for restricted second-order Møller-Plesset perturbation theory, as well as for both restricted and unrestricted Hartree-Fock and density functional theory. The importance of the three-body fully analytic gradient is illustrated through the failure of the two-body FMO method during molecular dynamics simulations of a small water cluster. The parallel implementation of the fragment molecular orbital method, its parallel efficiency, and its scalability on the Blue Gene/Q architecture up to 262 144 CPU cores are also discussed. PMID:26913837
ERIC Educational Resources Information Center
Seddon, G. M.; Shubber, K. E.
1984-01-01
Investigated use of overhead transparencies containing a sequence of diagrams to represent a three-dimensional structure at different stages during a rotation. Significant learning occurred among 120 Bahraini boys (aged 15-16) when the transparencies contained multi-colored diagrams which were exposed simultaneously or individually in a cumulative…
NASA Technical Reports Server (NTRS)
Luther, George W., III
1987-01-01
In this paper, molecular orbital theory is used to explain a heterogeneous reaction mechanism for both pyrite oxidation and reduction. The mechanism demonstrates that the oxidation of FeS2 by Fe(3+) may occur as a result of three important criteria: (1) the presence of a suitable oxidant having a vacant orbital (in case of liquid phase) or site (solid phase) to bind to the FeS2 via sulfur; (2) the initial formation of a persulfido (disulfide) bridge between FeS2 and the oxidant, and (3) an electron transfer from a pi(asterisk) orbital in S2(2-) to a pi or pi(asterisk) orbital of the oxidant.
Multi-electron coincidence spectroscopy: double photoionization from molecular inner-shell orbitals
NASA Astrophysics Data System (ADS)
Hikosaka, Y.; Lablanquie, P.; Penent, F.; Nakano, M.; Ito, K.
2014-04-01
We have studied double photoionization from molecular inner-shell orbitals and investigated the properties of the resultant double core-hole states in molecules, by multi-electron coincidence spectroscopy with a magnetic bottle electron spectrometer. A brief summary of our previous studies is presented.
Martay, Hugo E. L.; England, Duncan G.; McCabe, David J.; Walmsley, Ian A.
2011-04-15
The complexity of ultrafast molecular photoionization presents an obstacle to the modeling of pump-probe experiments. Here, a simple optimized model of atomic rubidium is combined with a molecular dynamics model to predict quantitatively the results of a pump-probe experiment in which long-range rubidium dimers are first excited, then ionized after a variable delay. The method is illustrated by the outline of two proposed feasible experiments and the calculation of their outcomes. Both of these proposals use Feshbach {sup 87}Rb{sub 2} molecules. We show that long-range molecular pump-probe experiments should observe spin-orbit precession given a suitable pump pulse, and that the associated high-frequency beat signal in the ionization probability decays after a few tens of picoseconds. If the molecule was to be excited to only a single fine-structure state, then a low-frequency oscillation in the internuclear separation would be detectable through the time-dependent ionization cross section, giving a mechanism that would enable observation of coherent vibrational motion in this molecule.
Atomlike, Hollow-Core–Bound Molecular Orbitals of C₆₀
Feng, Min; Zhao, Jin; Petek, Hrvoje
2008-04-17
The research described in this product was performed in part in the Environmental Molecular Sciences Laboratory, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory. The atomic electron orbitals that underlie molecular bonding originate from the central Coulomb potential of the atomic core. We used scanning tunneling microscopy and density functional theory to explore the relation between the nearly spherical shape and unoccupied electronic structure of buckminsterfullerene (C60) molecules adsorbed on copper surfaces. Besides the known p* antibonding molecular orbitals of the carbon-atom framework, above 3.5 electron volts we found atomlike orbitals bound to the core of the hollow C60 cage. These “superatom” states hybridize like the s and p orbitals of hydrogen and alkali atoms into diatomic molecule-like dimers and free-electron bands of one-dimensional wires and two-dimensional quantum wells in C60 aggregates. We attribute the superatom states to the central potential binding an electron to its screening charge, a property expected for hollow-shell molecules derived from layered materials.
A Frontier Molecular Orbital determination of the active sites on dispersed metal catalysts
Augustine, R.L.; Lahanas, K.M.
1992-01-01
An angular overlap calculation has been used to determine the s, p and d orbital energy levels of the different types of surface sites present on a dispersed metal catalysts. The basis for these calculations is the reported finding that a large number of catalyzed reactions take place on single atom active sites on the metal surface. Thus, these sites can be considered as surface complexes made up of the central active atom surrounded by near-neighbor metal atom ligands'' with localized surface orbitals perturbed only by these ligands''. These complexes'' are based on a twelve coordinate species with the ligands'' attached to the t{sub 2g} orbitals and the coordinate axes coincident with the direction of the e{sub g} orbitals on the central atom. These data can permit a Frontier Molecular Orbital treatment of specific site activities as long as the surface orbital availability for overlap with adsorbed substrates is considered along with its energy value and symmetry.
A Frontier Molecular Orbital determination of the active sites on dispersed metal catalysts
Augustine, R.L.; Lahanas, K.M.
1992-11-01
An angular overlap calculation has been used to determine the s, p and d orbital energy levels of the different types of surface sites present on a dispersed metal catalysts. The basis for these calculations is the reported finding that a large number of catalyzed reactions take place on single atom active sites on the metal surface. Thus, these sites can be considered as surface complexes made up of the central active atom surrounded by near-neighbor metal atom ``ligands`` with localized surface orbitals perturbed only by these ``ligands``. These ``complexes`` are based on a twelve coordinate species with the ``ligands`` attached to the t{sub 2g} orbitals and the coordinate axes coincident with the direction of the e{sub g} orbitals on the central atom. These data can permit a Frontier Molecular Orbital treatment of specific site activities as long as the surface orbital availability for overlap with adsorbed substrates is considered along with its energy value and symmetry.
NASA Astrophysics Data System (ADS)
Kruijssen, J. M. Diederik; Dale, James E.; Longmore, Steven N.
2015-02-01
We recently proposed that the star-forming potential of dense molecular clouds in the Central Molecular Zone (CMZ, i.e. the central few 100 pc) of the Milky Way is intimately linked to their orbital dynamics, potentially giving rise to an absolute-time sequence of star-forming clouds. In this paper, we present an orbital model for the gas stream(s) observed in the CMZ. The model is obtained by integrating orbits in the empirically constrained gravitational potential and represents a good fit (χ _red^2=2.0) to the observed position-velocity distribution of dense (n > several 103 cm-3) gas, reproducing all of its key properties. The orbit is also consistent with observational constraints not included in the fitting process, such as the 3D space velocities of Sgr B2 and the Arches and Quintuplet clusters. It differs from previous, parametric models in several respects: (1) the orbit is open rather than closed due to the extended mass distribution in the CMZ, (2) its orbital velocity (100-200 km s-1) is twice as high as in previous models, and (3) Sgr A* coincides with the focus of the (eccentric) orbit rather than being offset. Our orbital solution supports the recently proposed scenario in which the dust ridge between G0.253+0.016 (`the Brick') and Sgr B2 represents an absolute-time sequence of star-forming clouds, of which the condensation was triggered by the tidal compression during their most recent pericentre passage. We position the clouds on a common timeline and find that their pericentre passages occurred 0.30-0.74 Myr ago. Given their short free-fall times (tff ˜ 0.34 Myr), the quiescent cloud G0.253+0.016 and the vigorously star-forming complex Sgr B2 are separated by a single free-fall time of evolution, implying that star formation proceeds rapidly once collapse has been initiated. We provide the complete orbital solution, as well as several quantitative predictions of our model (e.g. proper motions and the positions of star formation `hotspots'). The
Proppe, Jonny; Herrmann, Carmen
2015-02-01
Common trends in communication through molecular bridges are ubiquitous in chemistry, such as the frequently observed exponential decay of conductance/electron transport and of exchange spin coupling with increasing bridge length, or the increased communication through a bridge upon closing a diarylethene photoswitch. For antiferromagnetically coupled diradicals in which two equivalent spin centers are connected by a closed-shell bridge, the molecular orbitals (MOs) whose energy splitting dominates the coupling strength are similar in shape to the MOs of the dithiolated bridges, which in turn can be used to rationalize conductance. Therefore, it appears reasonable to expect the observed common property trends to result from common orbital trends. We illustrate based on a set of model compounds that this assumption is not true, and that common property trends result from either different pairs of orbitals being involved, or from orbital energies not being the dominant contribution to property trends. For substituent effects, an effective modification of the π system can make a comparison difficult. PMID:25382464
Attia, Ali Kamal; Souaya, Eglal R.; Soliman, Ethar A.
2015-01-01
Purpose: Thermal analysis techniques have been used to study the thermal behavior of dapoxetine and vardenafil hydrochlorides and confirmed using semi-empirical molecular orbital calculations. Methods: Thermogravimetric analysis, derivative thermogravimetry, differential thermal analysis and differential scanning calorimetry were used to determine the thermal behavior and purity of the drugs under investigation. Thermodynamic parameters such as activation energy, enthalpy, entropy and Gibbs free energy were calculated. Results: Thermal behavior of DAP and VAR were confirmed using by semi-empirical molecular orbital calculations. The purity values were found to be 99.97% and 99.95% for dapoxetine and vardenafil hydrochlorides, respectively. The purity of dapoxetine and vardenafil hydrochlorides is similar to that found by reported methods according to DSC data. Conclusion: Thermal analysis justifies its application in quality control of pharmaceutical compounds due to its simplicity, sensitivity and low operational costs. PMID:26819925
Mulholland, J.A.; Mukherjee, J.; Wornat, M.J.; Sarofim, A.F.; Rutledge, G.C. . Dept. of Chemical Engineering)
1993-08-01
The pyrolysis of pure anthracene at temperatures between 1,200 and 1,500 K produced all six bianthryl isomers whose relative yields appear to be related to steric factors. To evaluate the hypothesis that thermodynamic factors govern the product distribution of bianthryls in this system, the relative enthalpies and entropies of biaryl isomers were estimated by molecular orbital modeling, using the semiempirical AM1 (Austin Model 1). Computational analysis of several isomer sets demonstrates that the relative stabilities of a large number of biaryl isomers are determined largely by steric interactions caused by structural features defined as bays, coves, and fjords. These steric factors affect both the degree of biaryl twist in the preferred conformation and the freedom of internal rotation. Molecular orbital modeling supports the hypothesis that a thermodynamic distribution of bianthryl isomers is produced by anthracene pyrolysis.
NASA Astrophysics Data System (ADS)
Romere, P. O.
1982-03-01
A proposed configuration for a Space Operations Center is presented in its eight stages of buildup. The on orbit aerodynamic force and moment characteristics were calculated for each stage based upon free molecular flow theory. Calculation of the aerodynamic characteristics was accomplished through the use of an orbital aerodynamic computer program, and the computation method is described with respect to the free molecular theory used. The aerodynamic characteristics are presented in tabulated form for each buildup stage at angles of attack from 0 to 360 degrees and roll angles from -60 to +60 degrees. The reference altitude is 490 kilometers, however, the data should be applicable for altitudes below 490 kilometers down to approximately 185 kilometers.
Effect of vacuum processing on outgassing within an orbiting molecular shield
NASA Technical Reports Server (NTRS)
Outlaw, R. A.
1982-01-01
The limiting hydrogen number density in an orbiting molecular shield is highly dependent on the outgassing rates from the materials of construction for the shield, experimental apparatus, and other hardware contained within the shield. Ordinary degassing temperatures used for ultrahigh vacuum studies (less than 450 C) are not sufficient to process metals so that the contribution to the number density within the shield due to outgassing is less than the theoretically attainable level (approximately 200 per cu. cm). Pure aluminum and type 347 stainless steel were studied as candidate shield materials. Measurements of their hydrogen concentration and diffusion coefficients were made, and the effects of high temperature vacuum processing (greater than 600 C) on their resulting outgassing rates was determined. The densities in a molecular shield due to the outgassing from either metal were substantially less ( 0.003) than the density due to the ambient atomic hydrogen flux at an orbital altitude of 500 km.
NASA Technical Reports Server (NTRS)
Romere, P. O.
1982-01-01
A proposed configuration for a Space Operations Center is presented in its eight stages of buildup. The on orbit aerodynamic force and moment characteristics were calculated for each stage based upon free molecular flow theory. Calculation of the aerodynamic characteristics was accomplished through the use of an orbital aerodynamic computer program, and the computation method is described with respect to the free molecular theory used. The aerodynamic characteristics are presented in tabulated form for each buildup stage at angles of attack from 0 to 360 degrees and roll angles from -60 to +60 degrees. The reference altitude is 490 kilometers, however, the data should be applicable for altitudes below 490 kilometers down to approximately 185 kilometers.
Imoto, Hideo; Saito, Taro; Adachi, Hirohiko
1995-04-26
Discrete variational-{Chi}{alpha} molecular orbital methods were applied to octahedral cluster complexes, [Mo{sub 6}X{sub 8}-(PH{sub 3}){sub 6}](X = S and Se). This structure is of interest due to its role in superconductivity of Chevrel plates. Level energies are discussed and factors contributing to their separations are categorized. Agreement with empirical XPS data is excellent.
Spectral studies and molecular orbital PPP-calculations of some azo-dyes
NASA Astrophysics Data System (ADS)
Shalabi, A. S.; Dessouki, H. A.; Issa, Y. M.; Ahmed, I. S.
2002-10-01
The UV, IR and 1H-NMR spectra of some 4-( R-phenyl azo) 1-hydroxy 2-naphthoic acid derivatives are studied. The effects of substituent groups and the solvent polarity on electronic spectral, IR bands and 1H-NMR proton chemical shifts are considered, the molecular orbital calculations obtained are rationalized quantitatively with that obtained practically using the PPP-model with configuration interaction (CI).
Mass transport properties of Pu/DT mixtures from orbital free molecular dynamics simulations
Kress, Joel David; Ticknor, Christopher; Collins, Lee A.
2015-09-16
Mass transport properties (shear viscosity and diffusion coefficients) for Pu/DT mixtures were calculated with Orbital Free Molecular Dynamics (OFMD). The results were fitted to simple functions of mass density (for ρ=10.4 to 62.4 g/cm^{3}) and temperature (for T=100 up to 3,000 eV) for Pu/DT mixtures consisting of 100/0, 25/75, 50/50, and 75/25 by number.
Sherman, David M.
1986-01-01
A molecular orbital description, based on spin-unrestricted X??-scattered wave calculations, is given for the electronic structures of mixed valence iron oxides and silicates. The cluster calculations show that electron hopping and optical intervalence charge-transger result from weak FeFe bonding across shared edges of FeO6 coordination polyhedra. In agreement with Zener's double exchange model, FeFe bonding is found to stabilize ferromagnetic coupling between Fe2+ and Fe3+ cations. ?? 1986.
The interaction of hydrogen with an Fe vacancy: a molecular orbital simulation
NASA Astrophysics Data System (ADS)
Irigoyen, B.; Ferullo, R.; Castellani, N.; Juan, A.
1996-05-01
The semi-empirical atom superposition and electron delocalization molecular orbital (ASED-MO) theory was used to study the H - Fe interaction present during the adsorption of monotomic hydrogen on Fe(110) and the subsequent absorption of the hydrogen atom in a bulk Fe vacancy. For that purpose calculations were carried out using clusters with 53 and 50 atoms, respectively. Results indicate that, in general, the hydrogen - iron interaction is stronger near the surface and near the vacancy.
Ab initio derivation of multi-orbital extended Hubbard model for molecular crystals
NASA Astrophysics Data System (ADS)
Tsuchiizu, Masahisa; Omori, Yukiko; Suzumura, Yoshikazu; Bonnet, Marie-Laure; Robert, Vincent
2012-01-01
From configuration interaction (CI) ab initio calculations, we derive an effective two-orbital extended Hubbard model based on the gerade (g) and ungerade (u) molecular orbitals (MOs) of the charge-transfer molecular conductor (TTM-TTP)I3 and the single-component molecular conductor [Au(tmdt)2]. First, by focusing on the isolated molecule, we determine the parameters for the model Hamiltonian so as to reproduce the CI Hamiltonian matrix. Next, we extend the analysis to two neighboring molecule pairs in the crystal and we perform similar calculations to evaluate the inter-molecular interactions. From the resulting tight-binding parameters, we analyze the band structure to confirm that two bands overlap and mix in together, supporting the multi-band feature. Furthermore, using a fragment decomposition, we derive the effective model based on the fragment MOs and show that the staking TTM-TTP molecules can be described by the zig-zag two-leg ladder with the inter-molecular transfer integral being larger than the intra-fragment transfer integral within the molecule. The inter-site interactions between the fragments follow a Coulomb law, supporting the fragment decomposition strategy.
Pseudo-symmetry analysis of the d-block molecular orbitals in four-coordinate complexes.
Falceto, Andrés; Casanova, David; Alemany, Pere; Alvarez, Santiago
2013-06-01
A rigorous definition of the concept of pseudo-symmetry, which is as important to chemistry as the concepts of symmetry implemented through group theory, should allow us to apply those group theoretical tools to molecules that are significantly distorted from those ideal symmetries best known and understood by the chemical community. In this paper, we consider four-coordinate transition-metal complexes with geometries along the interconversion path between the square and the tetrahedron and show how their molecular orbitals can be expressed in terms of either the tetrahedral or tetragonal symmetry groups. Furthermore, we analyze how the intensity of a d-d absorption band can be related to the degree of symmetry loss of the d-block molecular orbitals by means of their decomposition in terms of contributions from different pseudo-symmetry representations. As a final example, we also show how the substitution of a single ligand in a square planar complex affects the symmetry of the molecular orbitals and the absorption intensity associated to an electronic transition. PMID:23668721
Molecular orbital study of some eight-coordinate sulfur chelate complexes of molybdenum
Perkins, P.G.; Schultz, F.A.
1983-03-30
A number of molybdenum complexes involving the formal oxidation states Mo(IV) and Mo(V) have been studied by a self-consistent-field molecular orbital technique. All the complexes were of dodecahedral geometry and had eight sulfurs chelated to the central metal atom. In all, a series of five tetrakis complexes was studied, including the ligands dithiocarbamate (dtc), thioxanthate (txn), 1,1-dicyano-2,2-ethylenedithiolate (i-mnt), 1-cyano-1-carbethoxy-2,2-ethylenedithiolate (ced), and 1,1-dicarbethoxy-2,2-ethylenedithiolate (ded). The 4d orbitals were included on molybdenum, and the empty 3d levels on all sulfur atoms. The results show that the highest occupied molecular orbital in each case has over 90% metal d/sub xy/ character. Further, the energy of this orbital is linearly related to the reversible half-wave potentials for Mo(IV) ..-->.. Mo(V) and Mo(V) ..-->.. Mo(VI) oxidations of the complexes. A further irreversible oxidation observed experimentally also is closely related to the calculated energy levels. Relationships between the calculated results and Mo 3d/sub 5///sub 2/ X-ray photoelectron binding energies, EPR parameters, and charge-transfer absorption energies are discussed. Electrochemical and spectroscopic properties of these MoS/sub 8/ complexes can be understood in terms of a manifold of orbital energies that retain approximately constant spacings between one another and that move up or down in absolute energy in response to the charge donated or withdrawn by the ligands.
Meyer, Benjamin; Guillot, Benoît; Ruiz-Lopez, Manuel F; Jelsch, Christian; Genoni, Alessandro
2016-03-01
Due to both technical and methodological difficulties, determining and analyzing charge densities of very large molecular systems represents a serious challenge that, in the crystallographers community, has been mainly tackled by observing that the so-called pseudoatoms of the electron density multipole expansions are reliably transferable from molecule to molecule. This has led to the construction of pseudoatoms databanks that have allowed successful refinements of crystallographic structures of macromolecules, while taking into account their corresponding reconstructed electron distributions. A recent alternative/complement to the previous approach is represented by techniques based on extremely localized molecular orbitals (ELMOs) that, due to their strict localization on small molecular fragments (e.g., atoms, bonds, and functional groups), are also in principle exportable from system to system. The ELMOs transferability has been already tested in detail, and, in this work, it has been compared to the one of the pseudoatoms. To accomplish this task, electron distributions obtained both through the transfer of pseudoatoms and through the transfer of extremely localized molecular orbitals have been analyzed, especially taking into account topological properties and similarity indexes. The obtained results indicate that all the considered reconstruction methods give completely reasonable and similar charge densities, and, consequently, the new ELMOs libraries will probably represent new useful tools not only for refining crystal structures but also for computing approximate electronic properties of very large molecules. PMID:26799595
Néel, Nicolas; Lattelais, Marie; Bocquet, Marie-Laure; Kröger, Jörg
2016-02-23
Single-molecule chemistry with a scanning tunneling microscope has preponderantly been performed on metal surfaces. The molecule-metal hybridization, however, is often detrimental to genuine molecular properties and obscures their changes upon chemical reactions. We used graphene on Ir(111) to reduce the coupling between Ir(111) and adsorbed phthalocyanine molecules. By local electron injection from the tip of a scanning tunneling microscope the two pyrrolic H atoms were removed from single phthalocyanines. The detachment of the H atom pair induced a strong modification of the molecular electronic structure, albeit with no change in the adsorption geometry. Spectra and maps of the differential conductance combined with density functional calculations unveiled the entire depopulation of the highest occupied molecular orbital upon H abstraction. Occupied π states of intact molecules are proposed to be emptied via intramolecular electron transfer to dangling σ states of H-free N atoms. PMID:26812093
An analytical molecular orbital approach in tetrathiafulvalene tetracyanoquinodimethane (TTF-TCNQ)
NASA Astrophysics Data System (ADS)
Tada, Tomofumi; Aoki, Yuriko; Imamura, Akira
In low-dimensional molecular crystals exhibiting the Peierls instability, intercolumnar interactions play an important role in electrical conduction. In this work, we propose an efficient method based on molecular orbital theory for studying the Peierls instability in molecular crystals composed of mixed-valence complexes. A perturbational approach to general N-merized systems is incorporated and applied to tetrathiafulvalene tetracyanoquinodimethane (TTF-TCNQ) crystals. From the results of the calculations, a mixed state of several commensurate distortions is proposed as another aspect of the incommensurate distortion in TTF-TCNQ under atmospheric pressure, and TTF+0.5-TCNQ-0.5 crystal is found to be expected as an electronic conductor under high pressure.
Transport properties of lithium hydride at extreme conditions from orbital-free molecular dynamics
NASA Astrophysics Data System (ADS)
Burakovsky, L.; Ticknor, C.; Kress, J. D.; Collins, L. A.; Lambert, F.
2013-02-01
We have performed a systematic study of lithium hydride (LiH), using orbital-free molecular dynamics, with a focus on mass transport properties such as diffusion and viscosity by extending our previous studies at the lower end of the warm, dense matter regime to cover a span of densities from ambient to 10-fold compressed and temperatures from 10 eV to 10 keV. We determine analytic formulas for self- and mutual-diffusion coefficients, and viscosity, which are in excellent agreement with our molecular dynamics results, and interpolate smoothly between liquid and dense plasma regimes. In addition, we find the orbital-free calculations begin to agree with the Brinzinskii-Landau formula above about 250 eV at which point the medium becomes fully ionized. A binary-ion model based on a bare Coulomb interaction within a neutralizing background with the effective charges determined from a regularization prescription shows good agreement above about 100 eV with the orbital-free results. Finally, we demonstrate the validity of a pressure-based mixing rule in determining the transport properties from the pure-species quantities.
NASA Technical Reports Server (NTRS)
Richmond, R. G.; Kelso, R. M.
1980-01-01
A concern has arisen regarding the emissive distribution of water molecules from the shuttle orbiter flash evaporator system (FES). The role of the orbiter fuselage and elevon in affecting molecular scattering distributions was nuclear. The effect of these components were evaluated. Molecular distributions of the water vapor effluents from the FE were measured. These data were compared with analytically predicted values and the resulting implications were calculated.
NASA Astrophysics Data System (ADS)
Chaston, Scot
1999-02-01
Thermodynamic data such as equilibrium constants, standard cell potentials, molar enthalpies of formation, and standard entropies of substances can be a very useful basis for an organized presentation of knowledge in diverse areas of applied chemistry. Thermodynamic data can become particularly useful when incorporated into thermodynamic diagrams that are designed to be easy to recall, to serve as a basis for reconstructing previous knowledge, and to determine whether reactions can occur exergonically or only with the help of an external energy source. Few students in our chemistry-based courses would want to acquire the depth of knowledge or rigor of professional thermodynamicists. But they should nevertheless learn how to make good use of thermodynamic data in their professional occupations that span the chemical, biological, environmental, and medical laboratory fields. This article discusses examples of three thermodynamic diagrams that have been developed for this purpose. They are the thermodynamic energy account (TEA), the total entropy scale, and the thermodynamic scale diagrams. These diagrams help in the teaching and learning of thermodynamics by bringing the imagination into the process of developing a better understanding of abstract thermodynamic functions, and by allowing the reader to keep track of specialist thermodynamic discourses in the literature.
A molecular orbital study of H interaction with an edge dislocation in Fe bcc
NASA Astrophysics Data System (ADS)
Irigoyen, B.; Ferullo, R.; Castellani, N.; Juan, A.
1997-07-01
We have calculated the total energy curve for an H atom near an edge dislocation interacting with the surface. We have used a cluster to simulate the Fe(100) surface and an edge (100) dislocation. The semi-empirical method based on atom superposition and electron delocalization molecular orbital (ASED-MO) theory was employed. In the ground state, the distorted cluster expands while the H occupies an eccentric position. The H energy in that position is 0.15 eV lower than at the surface. The electronic structure was analysed by local density of states (LDOS). The width of the metal band is increased near the dislocation void because of a higher hybridization and the interaction with the H 1s orbital.
NASA Astrophysics Data System (ADS)
Fukushima, Kimichika
1993-07-01
Cluster calculations using the molecular orbital method were performed to search for impurities which would act as pinning centers for magnetic flux in copper oxide superconductors. Electronic structures were first investigated for interstitial helium atoms and for He atoms (including He clusters) substituted for oxygen and copper atoms. Calculations were then done for F atoms substituted for O atoms, and C, Mg, Al, Zn, Ga, and Cd atoms substituted for Cu. With these impurities, the energies of orbitals mainly attributable to the impurities differ from the Fermi energy. The Mulliken populations show that most charge carriers do not exist on the impurity sites (though some charge carriers enter the C sites), so it can be expected that the impurities considered here will act as pinning centers.
Multi-Orbital Molecular Compound (TTM-TTP)I3: Effective Model and Fragment Decomposition
NASA Astrophysics Data System (ADS)
Tsuchiizu, Masahisa; Omori, Yukiko; Suzumura, Yoshikazu; Bonnet, Marie-Laure; Robert, Vincent; Ishibashi, Shoji; Seo, Hitoshi
2011-01-01
The electronic structure of the molecular compound (TTM-TTP)I3, which exhibits a peculiar intra-molecular charge ordering, has been studied using multi-configuration ab initio calculations. First we derive an effective Hubbard-type model based on the molecular orbitals (MOs) of TTM-TTP; we set up a two-orbital Hamiltonian for the two MOs near the Fermi energy and determine its full parameters: the transfer integrals, the Coulomb and exchange interactions. The tight-binding band structure obtained from these transfer integrals is consistent with the result of the direct band calculation based on density functional theory. Then, by decomposing the frontier MOs into two parts, i.e., fragments, we find that the stacked TTM-TTP molecules can be described by a two-leg ladder model, while the inter-fragment Coulomb energies are scaled to the inverse of their distances. This result indicates that the fragment picture that we proposed earlier [M.-L. Bonnet et al.: J. Chem. Phys. 132 (2010) 214705] successfully describes the low-energy properties of this compound.
Detection of a Molecular Disk Orbiting the Nearby, "old," Classical T Tauri Star MP Muscae
NASA Astrophysics Data System (ADS)
Kastner, Joel H.; Hily-Blant, Pierry; Sacco, G. G.; Forveille, Thierry; Zuckerman, B.
2010-11-01
We have used the Atacama Pathfinder Experiment 12 m telescope to detect circumstellar CO emission from MP Muscae (MP Mus; K1 IVe), a nearby (D ~ 100 pc), actively accreting, ~7 Myr old pre-main-sequence (pre-MS) star. The CO emission line profile measured for MP Mus is indicative of an orbiting disk with radius ~120 AU, assuming that the central star mass is 1.2 M sun and the disk inclination is i ~ 30°. The inferred disk molecular gas mass is ~3 M ⊕. MP Mus thereby joins TW Hya and V4046 Sgr as the only late-type (low-mass), pre-MS star systems within ~100 pc of Earth that are known to retain orbiting, molecular disks. We also report the nondetection (with the Institut de Radio Astronomie Millimetrique 30 m telescope) of CO emission from another 10 nearby (D <~ 100 pc), dusty, young (age ~10-100 Myr) field stars of spectral type A-G. We discuss the implications of these results for the timescales for stellar and Jovian planet accretion from, and dissipation of, molecular disks around young stars. This research is based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere, Chile, proposal number 385.C-0143, with the Atacama Pathfinder Experiment (APEX).
DETECTION OF A MOLECULAR DISK ORBITING THE NEARBY, 'OLD', CLASSICAL T TAURI STAR MP MUSCAE
Kastner, Joel H.; Sacco, G. G.; Hily-Blant, Pierry; Forveille, Thierry; Zuckerman, B.
2010-11-10
We have used the Atacama Pathfinder Experiment 12 m telescope to detect circumstellar CO emission from MP Muscae (MP Mus; K1 IVe), a nearby (D {approx} 100 pc), actively accreting, {approx}7 Myr old pre-main-sequence (pre-MS) star. The CO emission line profile measured for MP Mus is indicative of an orbiting disk with radius {approx}120 AU, assuming that the central star mass is 1.2 M {sub sun} and the disk inclination is i {approx} 30{sup 0}. The inferred disk molecular gas mass is {approx}3 M {sub +}. MP Mus thereby joins TW Hya and V4046 Sgr as the only late-type (low-mass), pre-MS star systems within {approx}100 pc of Earth that are known to retain orbiting, molecular disks. We also report the nondetection (with the Institut de Radio Astronomie Millimetrique 30 m telescope) of CO emission from another 10 nearby (D {approx_lt} 100 pc), dusty, young (age {approx}10-100 Myr) field stars of spectral type A-G. We discuss the implications of these results for the timescales for stellar and Jovian planet accretion from, and dissipation of, molecular disks around young stars.
Radical damage in lipids investigated with the fragment molecular orbital method
NASA Astrophysics Data System (ADS)
Green, Mandy C.; Nakata, Hiroya; Fedorov, Dmitri G.; Slipchenko, Lyudmila V.
2016-05-01
To quantify the thermodynamics for hydrogen abstraction lipids, the fragment molecular orbital method (FMO) is used to calculate structures and energies of the reactants and products. The analytic second derivative is developed for the open-shell Hartree-Fock formulation of FMO and used to calculate zero point energy corrections. The accuracy of FMO is evaluated for a lipid model and the errors in reaction energies are found not to exceed 0.5 kcal/mol. The reaction energies determined for multiple sites in two lipids are used to discuss likely sites and pathways of radical initiation in membranes.
Li, Shaowei; Yuan, Dingwang; Yu, Arthur; Czap, Gregory; Wu, Ruqian; Ho, W
2015-05-22
A hydrogen molecule can diffuse freely on the surface and be trapped above an adsorbed molecule within the junction of a scanning tunneling microscope. The trapped dihydrogen exhibits the properties of a free rotor. Here we show that the intermolecular interaction between dihydrogen and Mg-porphyrin (MgP) can be visualized by imaging j=0 to 2 rotational excitation of dihydrogen. The interaction leads to a weakened H-H bond and modest electron donation from the dihydrogen to the lowest unoccupied molecular orbital of MgP, a process similarly observed for the interaction between dihydrogen and an adsorbed Au atom. PMID:26047242
NASA Astrophysics Data System (ADS)
Watanabe, Hirofumi; Okiyama, Yoshio; Nakano, Tatsuya; Tanaka, Shigenori
2010-11-01
We developed FMO-PB method, which incorporates solvation effects into the Fragment Molecular Orbital calculation with the Poisson-Boltzmann equation. This method retains good accuracy in energy calculations with reduced computational time. We calculated the solvation free energies for polyalanines, Alpha-1 peptide, tryptophan cage, and complex of estrogen receptor and 17 β-estradiol to show the applicability of this method for practical systems. From the calculated results, it has been confirmed that the FMO-PB method is useful for large biomolecules in solution. We also discussed the electric charges which are used in solving the Poisson-Boltzmann equation.
Counterpoise-corrected interaction energy analysis based on the fragment molecular orbital scheme
NASA Astrophysics Data System (ADS)
Okiyama, Yoshio; Fukuzawa, Kaori; Yamada, Haruka; Mochizuki, Yuji; Nakano, Tatsuya; Tanaka, Shigenori
2011-06-01
Basis set superposition error (BSSE) correction with counterpoise (CP) procedure under the environmental electrostatic potential is newly introduced to interfragment interaction energy (IFIE), which is important for interaction analysis in the fragment molecular orbital method. The CP correction for IFIE is applied to a stacked dimer of base pair and a protein-ligand complex of estrogen receptor and 17β-estradiol with scaled third-order Møller-Plesset perturbation theory. The BSSEs amount to about quarter of IFIE for hydrogen-bonding and electrostatic interactions and half or even more for dispersion interactions. Estimation of IFIE with the CP correction is therefore preferred for the quantitative discussion.
Multicenter molecular integrals for Slater orbitals of higher principal quantum numbers
NASA Technical Reports Server (NTRS)
Tai, H.
1989-01-01
As was shown earlier by Tai (1979), by using the Fourier-transform technique and properly coupling a pair of two-center exchange integrals, the multicenter molecular integrals can be cast into a simple expression upon which numerical procedures can be directly applied. In this paper, the procedure of Tai is extended to integrals involving orbitals with arbitrarily higher principal quantum number. The derivation is outlined, and the explicit expressions are presented for a three-center nuclear attraction integral and a four-center two-electron Coulomb repulsion integral of arbitrary higher states.
NASA Astrophysics Data System (ADS)
Li, Shaowei; Yuan, Dingwang; Yu, Arthur; Czap, Gregory; Wu, Ruqian; Ho, W.
2015-05-01
A hydrogen molecule can diffuse freely on the surface and be trapped above an adsorbed molecule within the junction of a scanning tunneling microscope. The trapped dihydrogen exhibits the properties of a free rotor. Here we show that the intermolecular interaction between dihydrogen and Mg-porphyrin (MgP) can be visualized by imaging j =0 to 2 rotational excitation of dihydrogen. The interaction leads to a weakened H-H bond and modest electron donation from the dihydrogen to the lowest unoccupied molecular orbital of MgP, a process similarly observed for the interaction between dihydrogen and an adsorbed Au atom.
Dyadyusha, G.G.; Ushomirskii, M.M.
1986-09-01
A method previously proposed for determining the energy structure of a polymethyne dye with any terminal groups is used in considering formulas for the molecularorbital coefficients and the differences in the distribution on the atoms in the polymethyne chain for localized and delocalized energy levels, as well as the accuracy in calculating the molecular-orbital coefficients by means of a finite number of effective additive parameters. It is found that the localized states are important to the electron-density distribution on the chain atoms characteristic of the polymethyne dyes.
Fedorov, Dmitri G; Kitaura, Kazuo
2009-11-01
We have examined the role of the exchange in describing the electrostatic potential in the fragment molecular orbital method and showed that it should be included in the total Fock matrix to obtain an accurate one-electron spectrum; however, adding it to the Fock matrices of individual fragments and pairs leads to very large errors. For the error analysis we have used the virial theorem; numerical tests have been performed for solvated phenol at the Hartree-Fock level with the 6-31G( *) and 6-311G( * *) basis sets. PMID:19894991
Chiang, Ming-Hsi; Antonio, Mark R; Soderholm, L
2004-11-01
The ground state electronic properties of metal-exchanged Preyssler heteropolyoxoanions [M(n+)P(5)W(30)O(110)](n-15), in which the encapsulated M(n+) ions are the spherical, diamagnetic ions Na(+), Ca(2+), Sr(2+), Y(3+), La(3+) and Th(4+), are studied using a combination of electrochemical, optical, and NMR experiments. We have designed experiments that focus on the influence of the charge (n) of the encapsulated cations, which themselves have no redox response, and its effect upon the W-O framework MOs. As n increases, the cluster anions accept electrons into their LUMOs with increasing ease, and their lowest-energy LMCT bands reveal a corresponding blue shift, which is indicative of an increase of the LUMO-HOMO energy splitting with increasing n. (183)W NMR spectra are used to identify the atomic origin of the LUMO states, which are shown to be composed primarily of orbitals from the ring of 5 W atoms near M(n+). The cation charge correlates directly and linearly with the half-wave potentials of the first redox couples, the LMCT band energies, and the W chemical shifts. We have combined this suite of experimental results to construct an energy level diagram of the frontier MOs for the Preyssler cluster anions. In so doing, we provide a fundamental perspective that is not otherwise available on the cation's role with specific regard to the electronic behavior of the W-O orbitals. These results are expected to provide benchmarking information as theorists begin to study these large POM systems. PMID:15510277
NASA Astrophysics Data System (ADS)
Bauer, C.; Bauduin, P.; Dufrêche, J. F.; Zemb, T.; Diat, O.
2012-11-01
We consider the class of surfactants called "extractants" since they specifically interact with some cations and are used in liquid-liquid separation processes. We review here features of water-poor reverse micelles in water/oil/ extractant systems as determined by combined structural studies including small angle scattering techniques on absolute scale. Origins of instabilities, liquid-liquid separation as well as emulsification failure are detected. Phase diagrams contain the same multi-phase domains as classical microemulsions, but special unusual features appear due to the high spontaneous curvature directed towards the polar cores of aggregates as well as rigidity of the film made by extracting molecules.
Modeling of hydroxyapatite-peptide interaction based on fragment molecular orbital method
NASA Astrophysics Data System (ADS)
Kato, Koichiro; Fukuzawa, Kaori; Mochizuki, Yuji
2015-06-01
We have applied the four-body corrected fragment molecular orbital (FMO4) calculations to analyze the interaction between a designed peptide motif (Glu1-Ser2-Gln3-Glu4-Ser5) and the hydroxyapatite (HA) solid mimicked by a cluster model consisting of 1408 atoms. To incorporate statistical fluctuations, a total of 30 configurations were generated through classical molecular dynamics simulation with water molecules and were subjected to FMO4 calculations at the MP2 level. It was found that Ser5 plays a leading role in interacting with the phosphate moieties of HA via charge transfer and also that negatively charged Glu1 and Glu4 provide electrostatic stabilizations with the calcium ions.
Liquid Be, Ca and Ba. An orbital-free ab-initio molecular dynamics study
Rio, B. G. del; González, L. E.
2015-08-17
Several static and dynamic properties of liquid beryllium (l-Be), liquid calcium (l-Ca) and liquid barium (l-Ba) near their triple point have been evaluated by the orbital-free ab initio molecular dynamics method (OF-AIMD), where the interaction between valence electrons and ions is described by means of local pseudopotentials. These local pseudopotentials used were constructed through a force-matching process with those obtained from a Kohn-Sham ab initio molecular dynamics study (KS-AIMD) of a reduced system with non-local pseudopotentials. The calculated static structures show good agreement with the available experimental data, including an asymmetric second peak in the structure factor which has been linked to the existence of a marked icosahedral short-range order in the liquid. As for the dynamic properties, we obtain collective density excitations whose associated dispersion relations exhibit a positive dispersion.
Chen, Ji; Ren, Xinguo; Li, Xin-Zheng; Alfè, Dario; Wang, Enge
2014-07-14
The finite-temperature phase diagram of hydrogen in the region of phase IV and its neighborhood was studied using the ab initio molecular dynamics (MD) and the ab initio path-integral molecular dynamics (PIMD). The electronic structures were analyzed using the density-functional theory (DFT), the random-phase approximation, and the diffusion Monte Carlo (DMC) methods. Taking the state-of-the-art DMC results as benchmark, comparisons of the energy differences between structures generated from the MD and PIMD simulations, with molecular and dissociated hydrogens, respectively, in the weak molecular layers of phase IV, indicate that standard functionals in DFT tend to underestimate the dissociation barrier of the weak molecular layers in this mixed phase. Because of this underestimation, inclusion of the quantum nuclear effects (QNEs) in PIMD using electronic structures generated with these functionals leads to artificially dissociated hydrogen layers in phase IV and an error compensation between the neglect of QNEs and the deficiencies of these functionals in standard ab initio MD simulations exists. This analysis partly rationalizes why earlier ab initio MD simulations complement so well the experimental observations. The temperature and pressure dependencies for the stability of phase IV were also studied in the end and compared with earlier results. PMID:25028021
Chen, Ji; Ren, Xinguo; Li, Xin-Zheng; Alfè, Dario; Wang, Enge
2014-07-14
The finite-temperature phase diagram of hydrogen in the region of phase IV and its neighborhood was studied using the ab initio molecular dynamics (MD) and the ab initio path-integral molecular dynamics (PIMD). The electronic structures were analyzed using the density-functional theory (DFT), the random-phase approximation, and the diffusion Monte Carlo (DMC) methods. Taking the state-of-the-art DMC results as benchmark, comparisons of the energy differences between structures generated from the MD and PIMD simulations, with molecular and dissociated hydrogens, respectively, in the weak molecular layers of phase IV, indicate that standard functionals in DFT tend to underestimate the dissociation barrier of the weak molecular layers in this mixed phase. Because of this underestimation, inclusion of the quantum nuclear effects (QNEs) in PIMD using electronic structures generated with these functionals leads to artificially dissociated hydrogen layers in phase IV and an error compensation between the neglect of QNEs and the deficiencies of these functionals in standard ab initio MD simulations exists. This analysis partly rationalizes why earlier ab initio MD simulations complement so well the experimental observations. The temperature and pressure dependencies for the stability of phase IV were also studied in the end and compared with earlier results.
ERIC Educational Resources Information Center
Ruddick, Kristie R.; Parrill, Abby L.; Petersen, Richard L.
2012-01-01
In this study, a computational molecular orbital theory experiment was implemented in a first-semester honors general chemistry course. Students used the GAMESS (General Atomic and Molecular Electronic Structure System) quantum mechanical software (as implemented in ChemBio3D) to optimize the geometry for various small molecules. Extended Huckel…
Shinagawa, Y
1992-02-01
Quantum chemical studies were applied to analyze the activities of catechol O-methyltransferase (COMT) inhibitors. Molecular orbital calculations of inhibitor molecules were made by semi-empirical molecular orbital calculations, CNDO/2 (complete neglect of differential overlap) methods. Regression analysis among theoretical reaction indices based on the frontier electron theory and COMT inhibitory activities were carried out. The COMT inhibitory actions of two series of inhibitors, a series of 1,5-substituted 3,4-dihydroxy benzenes and a series of substituted 3-hydroxy-4-methoxy benzenes, were investigated. The resulting regression equations contain two common reaction indices as regression variables: the electron density on the oxygen atom of the hydroxyl group and the super-delocalizability on the 5th carbon atom of the benzene ring. These two atomic positions are considered to play an important role in the interaction of these inhibitors with COMT. The hydroxyl of atomic position 3 is probably indispensable to the COMT inhibitory action by these inhibitors. PMID:1507526
Why many semiempirical molecular orbital theories fail for liquid water and how to fix them.
Welborn, Matthew; Chen, Jiahao; Wang, Lee-Ping; Van Voorhis, Troy
2015-05-01
Water is an extremely important liquid for chemistry and the search for more accurate force fields for liquid water continues unabated. Neglect of diatomic differential overlap (NDDO) molecular orbital methods provide and intriguing generalization of classical force fields in this regard because they can account both for bond breaking and electronic polarization of molecules. However, we show that most standard NDDO methods fail for water because they give an incorrect description of hydrogen bonding, water's key structural feature. Using force matching, we design a reparameterized NDDO model and find that it qualitatively reproduces the experimental radial distribution function of water, as well as various monomer, dimer, and bulk properties that PM6 does not. This suggests that the apparent limitations of NDDO models are primarily due to poor parameterization and not to the NDDO approximations themselves. Finally, we identify the physical parameters that most influence the condensed phase properties. These results help to elucidate the chemistry that a semiempirical molecular orbital picture of water must capture. We conclude that properly parameterized NDDO models could be useful for simulations that require electronically detailed explicit solvent, including the calculation of redox potentials and simulation of charge transfer and photochemistry. PMID:25766721
Orbital alignment at the internal interface of arylthiol functionalized CdSe molecular hybrids
Li, Zhi; Schlaf, Rudy; Mazzio, Katherine A.; Okamoto, Ken; Luscombe, Christine K.
2015-04-21
Organic-inorganic nanoparticle molecular hybrid materials are interesting candidates for improving exciton separation in organic solar cells. The orbital alignment at the internal interface of cadmium selenide (ArS-CdSe) hybrid materials functionalized with covalently attached arylthiolate moieties was investigated through X-ray photoemission spectroscopy (XPS) and ultraviolet photoemission spectroscopy (UPS). A physisorbed interface between arylthiol (ArSH) ligands and CdSe nanoparticles was also investigated for comparison. This interface was created via a multi-step thin film deposition procedure in-vacuo, where the surface was characterized after each experimental step. This enabled the direct comparison of ArSH/CdSe interfaces produced via physisorption and ArS-CdSe covalently attached hybrid materials, which rely on a chemical reaction for their synthesis. All material depositions were performed using an electrospray deposition, which enabled the direct injection of solution-originating molecular species into the vacuum system. This method allows XPS and UPS measurements to be performed immediately after deposition without exposure to the atmosphere. Transmission electron microscopy was used to determine the morphology and particle size of the deposited materials. Ultraviolet-visible spectroscopy was used to estimate the optical band gap of the CdSe nanoparticles and the HOMO-LUMO gap of the ArSH ligands. These experiments showed that hybridization via covalent bonds results in an orbital realignment at the ArSH/CdSe interface in comparison to the physisorbed interface. The orbital alignment within the hybrid caused a favorable electron injection barrier, which likely facilitates exciton-dissociation while preventing charge-recombination.
A regression based model was developed to determine whether highest occupied molecular orbital (HOMO) energies, calculated using Kohn-Sham orbital density functional theory (DFT) could be used to estimate the OH rate constants of hydrofluorocarbons (HFCS) and hydrofluoroethers (H...
Danel, J.-F.; Blottiau, P.; Kazandjian, L.; Piron, R.; Torrent, M.
2014-10-15
The applicability of quantum molecular dynamics to the calculation of the equation of state of a dense plasma is limited at high temperature by computational cost. Orbital-free molecular dynamics, based on a semiclassical approximation and possibly on a gradient correction, is a simulation method available at high temperature. For a high-Z element such as lutetium, we examine how orbital-free molecular dynamics applied to the equation of state of a dense plasma can be regarded as the limit of quantum molecular dynamics at high temperature. For the normal mass density and twice the normal mass density, we show that the pressures calculated with the quantum approach converge monotonically towards those calculated with the orbital-free approach; we observe a faster convergence when the orbital-free approach includes the gradient correction. We propose a method to obtain an equation of state reproducing quantum molecular dynamics results up to high temperatures where this approach cannot be directly implemented. With the results already obtained for low-Z plasmas, the present study opens the way for reproducing the quantum molecular dynamics pressure for all elements up to high temperatures.
NASA Astrophysics Data System (ADS)
McBride, Carl; Vega, Carlos
2002-12-01
A study of a rigid fully flexible fused hard sphere model [C. McBride, C. Vega, and L. G. MacDowell, Phys. Rev. E 64, 011703 (2001)] is extended to the smectic and solid branches of the phase diagram. Computer simulations have been performed for a completely rigid model composed of 15 fused hard spheres (15+0), a model of 15 fused hard spheres of which 2 monomers at one end of the model form a flexible tail (13+2), and a model consisting of 15 fused hard spheres with 5 monomers forming a flexible tail (10+5). For the 15+0 model the phase sequence isotropic-nematic-smectic A-columnar is found on compression, and the sequence solid-smectic A-nematic-isotropic on expansion. For the 13+2 model the phase sequence isotropic-nematic-smectic C is found on compression, and the sequence solid-smectic A-nematic-isotropic on expansion. For the 10+5 model the phase sequence isotropic-glass is found on compression. The expansion runs displayed the phase sequence solid-smectic A-isotropic. The introduction of flexibility was seen to stabilize the smectic A phase at the expense of the nematic phase.
Ghasemi, A S; Mashhadban, F; Hoseini-Alfatemi, S M; Sharifi-Rad, J
2015-01-01
Today the use of nanotubes (CNTs) is widely spread a versatile vector for drug delivery that can officiate as a platform for transporting a variety of bioactive molecules, such as drugs. In the present study, the interaction between the nanotube and anticancer drugs is investigated. Density functional theory (DFT) calculations were using the Gauss view and the complexes were optimized by B3LYP method using B3LYP/6-31G (d, p) and B3LYP/6-311++G (d, p) basis set in the gas phase and water solution at 298.15K. The calculated hikes' occupied molecular orbital (HOMO) and the lowest unoccupied (LUMO) energies Show that charge transfer occurs within the molecule. Furthermore, the effects of interactions on the natural bond orbital analysis (NBO) have been used to a deeper investigation into the studied compounds. These factors compete against each other to determine the adsorption behavior of the tube computer simulation is seen to be capable to optimize anticancer drug design. This review article mainly concentrates on the different protocols of loading anticancer drugs onto CNTs as well as how to control the anticancer drug release and cancer treatment. PMID:26718433
Molecular Orbital Simulations of Metal 1s2p Resonant Inelastic X-ray Scattering.
Guo, Meiyuan; Källman, Erik; Sørensen, Lasse Kragh; Delcey, Mickaël G; Pinjari, Rahul V; Lundberg, Marcus
2016-07-28
For first-row transition metals, high-resolution 3d electronic structure information can be obtained using resonant inelastic X-ray scattering (RIXS). In the hard X-ray region, a K pre-edge (1s→3d) excitation can be followed by monitoring the dipole-allowed Kα (2p→1s) or Kβ (3p→1s) emission, processes labeled 1s2p or 1s3p RIXS. Here the restricted active space (RAS) approach, which is a molecular orbital method, is used for the first time to study hard X-ray RIXS processes. This is achieved by including the two sets of core orbitals in different partitions of the active space. Transition intensities are calculated using both first- and second-order expansions of the wave vector, including, but not limited to, electric dipoles and quadrupoles. The accuracy of the approach is tested for 1s2p RIXS of iron hexacyanides [Fe(CN)6](n-) in ferrous and ferric oxidation states. RAS simulations accurately describe the multiplet structures and the role of 2p and 3d spin-orbit coupling on energies and selection rules. Compared to experiment, relative energies of the two [Fe(CN)6](3-) resonances deviate by 0.2 eV in both incident energy and energy transfer directions, and multiplet splittings in [Fe(CN)6](4-) are reproduced within 0.1 eV. These values are similar to what can be expected for valence excitations. The development opens the modeling of hard X-ray scattering processes for both solution catalysts and enzymatic systems. PMID:27398775
Zhou, Xin
1998-11-30
Photoemission intensities from the molecular orbitals of c(2x2)CO/Pt(111) over a wide photon energy range were measured and analyzed by the same methods developed for structural studies using core levels. The 4{sigma} orbital center of gravity is found to be concentrated between the C and O atoms, while that of the 5{sigma} orbital lies between the C atom and the Pt surface. The C 1s photoelectron diffraction was used to determine the adsorption geometry. The earlier ambiguity that multiple scattering is needed to correctly model a {chi} curve while single scattering is sufficient for understanding major peaks in the ARPEFS-FTS is clarified by studying the clean Ni(111) surface. In the normal emission case, several different combinations of scattering events have similar path length differences (PLDs), and can either cancel each other or enhance the corresponding FT peak. In the off-normal case the degeneracy is greatly reduced due to the lower degree of symmetry. In normal emission AR PEFS, up to third order multiple scattering is needed to describe fully both the {chi} curve and its FT spectrum. To improve the spectral resolution in the ARPEFS-FT analysis, several new spectral analysis methods are introduced. With both autocorrelation autoregression (ACAR) and autocorrelation eigenvector (ACE), we can produce a reliable power spectrum by following the order-closing procedure. The best spectra are usually obtained when the autocorrelation sequence is computed with lags up to half the data range. A simple way of determining surface adsorption sites is proposed as follows: First use a single scattering cluster for possible adsorption sites to construct the geometrical PLDs from the strong backscattering events; then compare these PLDs with those obtained from the ARPEFS-FT analysis of the experimental data. After the preferred adsorption site is determined, fine tune the interlayer distances according to the positional R-factor.
Saju, K K; Jayadas, N H; Vidyanand, S; James, J
2011-03-01
It has been established that the adhesion of cells on to the surfaces of orthopaedic implants depends on the ability of the surfaces to accommodate protein molecules. Hydroxyapatite coating and anodizing are the most common methods to make TiAl6V4 implants (Ti) more biocompatible. In this paper Spartan 02, a molecular dynamics software, is used to analyze and predict the bonding characteristics of Extra cellular matrix protein sequence arginine-glycine-aspartic acid (RGD) on a Hyrdoxyapatite (HA) coated Ti and an anodized Ti surface based on the property of its constituent atoms, their polarity (net electrostatic charge, Qr), the energies of the molecular orbital E_HOMO (energy of the highest occupied molecular orbital), and E_LUMO (energy of the lowest unoccupied molecular orbital). The results show favourable criterion for formation of bonding between the HOMO orbital of the HA coated and anodized surfaces and LUMO orbital of the glycine strand from the RGD unit. The mechanism of bonding of individual atoms to form primary calcium oxide compounds is likely only in the case of HA coated surfaces . The surface texture of the anodized Ti with inherent porosities appear more responsible for the adsorption of proteins on to them by mechanical interlocking than the formation of any intermediate calcium oxide compounds. PMID:21485326
Rakhmilevitch, David; Sarkar, Soumyajit; Bitton, Ora; Kronik, Leeor; Tal, Oren
2016-03-01
Molecular junctions based on ferromagnetic electrodes allow the study of electronic spin transport near the limit of spintronics miniaturization. However, these junctions reveal moderate magnetoresistance that is sensitive to the orbital structure at their ferromagnet-molecule interfaces. The key structural parameters that should be controlled in order to gain high magnetoresistance have not been established, despite their importance for efficient manipulation of spin transport at the nanoscale. Here, we show that single-molecule junctions based on nickel electrodes and benzene molecules can yield a significant anisotropic magnetoresistance of up to ∼200% near the conductance quantum G0. The measured magnetoresistance is mechanically tuned by changing the distance between the electrodes, revealing a nonmonotonic response to junction elongation. These findings are ascribed with the aid of first-principles calculations to variations in the metal-molecule orientation that can be adjusted to obtain highly spin-selective orbital hybridization. Our results demonstrate the important role of geometrical considerations in determining the spin transport properties of metal-molecule interfaces. PMID:26926769
Terao, Jun; Wadahama, Akihisa; Matono, Akitoshi; Tada, Tomofumi; Watanabe, Satoshi; Seki, Shu; Fujihara, Tetsuaki; Tsuji, Yasushi
2013-01-01
The feasibility of using π-conjugated polymers as next-generation electronic materials is extensively studied; however, their charge mobilities are lower than those of inorganic materials. Here we demonstrate a new design principle for increasing the intramolecular charge mobility of π-conjugated polymers by covering the π-conjugated chain with macrocycles and regularly localizing π-molecular orbitals to realize an ideal orbital alignment for charge hopping. Based on theoretical predictions, insulated wires containing meta-junctioned poly(phenylene–ethynylene) as the backbone units were designed and synthesized. The zigzag wires exhibited higher intramolecular charge mobility than the corresponding linear wires. When the length of the linear region of the zigzag wires was increased to 10 phenylene–ethynylene units, the intramolecular charge mobility increased to 8.5 cm2 V−1 s−1. Theoretical analysis confirmed that this design principle is suitable for obtaining ideal charge mobilities in π-conjugated polymer chains and that it provides the most effective pathways for inter-site hopping processes. PMID:23575695
Predicting the solid solubility limit in high-entropy alloys using the molecular orbital approach
NASA Astrophysics Data System (ADS)
Sheikh, Saad; Klement, Uta; Guo, Sheng
2015-11-01
High-entropy alloys (HEAs) are currently at the research frontier of metallic materials. Understanding the solid solubility limit in HEAs, such a highly concentrated multicomponent alloy system, is scientifically intriguing. It is also technically important to achieve desirable mechanical properties by controlling the formation of topologically or geometrically closed packed phases. Previous approaches to describe the solid solubilities in HEAs could not accurately locate the solubility limit and have to utilize at least two parameters. Here, we propose to use a single parameter, the average energy of d-orbital levels, Md, to predict the solid solubility limit in HEAs. It is found that Md can satisfactorily describe the solid solubilities in fcc structured HEAs containing 3 d transition metals, and also in bcc structured HEAs. This finding will greatly simplify the alloys design and lends more flexibility to control the mechanical properties of HEAs. When 4 d transition metals are alloyed, Md alone cannot describe the solid solubility limit in fcc structured HEAs, due to the large increase of the bond strength that can be gauged by the bond order, Bo. The potential opportunities and challenges with applying the molecular orbital approach to HEAs are discussed.
Oxygen evolution on a SrFeO3 anode - Mechanistic considerations from molecular orbital theory
NASA Technical Reports Server (NTRS)
Mehandru, S. P.; Anderson, Alfred B.
1989-01-01
Various pathways proposed in the literature for the evolution of O2 in electrochemical oxidations are explored using the atom superposition and electron delocalization molecular orbital (ASED-MO) theory and the cluster models of the SrFeO3 surface as a prototype material. Calculations indicate that oxygen atoms can be easily formed on the (100) surface as well as on the edge cation sites of a SrFeO3 anode by the discharge of OH(-), followed by its deprotonation and electron transfer to the electrode. The O atoms can form O2 on the edge and corner sites, where the Fe(4+) is coordinated to four and three bulk oxygen anions, respectively. The calculations strongly disfavor mechanisms involving coupling of oxygen atoms adsorbed on different cations as well as a mechanism featuring an ozone intermediate.
Fragment molecular orbital calculations on red fluorescent proteins (DsRed and mFruits).
Taguchi, Naoki; Mochizuki, Yuji; Nakano, Tatsuya; Amari, Shinji; Fukuzawa, Kaori; Ishikawa, Takeshi; Sakurai, Minoru; Tanaka, Shigenori
2009-01-29
We have performed a series of fragment molecular orbital (FMO) calculations for a family of red fluorescent proteins, DsRed and mFruits. The electronic transition energies were evaluated by the method of configuration interaction singles with perturbative doubles [CIS(D)] including higher-order corrections. The calculated values were in good agreement with the corresponding experimental peak values of spectra. Additionally, the chromophore environment was systematically analyzed in terms of the interaction energies between the pigment moiety and neighboring residues. It was theoretically revealed that the electrostatic interactions play a dominant role in the DsRed chromophore, whereas the color tunings in mFruits are controlled in a more delicate fashion. PMID:19127982
AM1 and ab initio molecular orbital study of water dimer
Dannenberg, J.J.
1988-12-01
Several structures for the water dimer, including trifurcated structures similar to the optimized AM1 geometry, have been calculated by using the MP4/6-311G** level of ab initio molecular orbital theory. The relative energies of the structures become quite close at the higher levels of calculation. The best trifurcated is only 0.2 kcal/mol higher than the optimized HF/6-31G* structure and only 0.4 kcal/mol higher than the lowest energy structure found (optimized by using AM1 with the H bond constrained to be linear). It appears likely that the potential surface of the water dimer is extremely flat. The experimental geometry, which corresponds to the minimum on the free energy surface, is likely to be dominated by entropy contributions.
Hyperfine Parameters for Aluminum Hydride: An ab Initio Molecular Orbital Study
NASA Astrophysics Data System (ADS)
Gee, Myrlene; Wasylishen, Roderick E.
2001-06-01
An extensive ab initio molecular orbital study of the 27Al nuclear spin-rotation and nuclear quadrupolar coupling constants in aluminum hydride, AlH, has been performed. The 27Al nuclear spin-rotation constant (C⊥), calculated to be approximately 300 kHz, was neglected in a previous analysis of the hyperfine structure in the microwave spectrum (M. Goto and S. Saito, Astrophys. J.452, L147-148 (1995)). Unfortunately, the ab initio calculations do not provide a definitive value for the aluminum nuclear quadrupolar coupling constant, but suggest a value of -49±4 MHz. It is apparent that the microwave study of AlH should be repeated.
Hybrid RHF/MP2 Geometry Optimizations with the Effective Fragment Molecular Orbital Method
Christensen, Anders S.; Steinmann, Casper; Fedorov, Dmitri G.; Jensen, Jan H.
2014-01-01
The frozen domain effective fragment molecular orbital method is extended to allow for the treatment of a single fragment at the MP2 level of theory. The approach is applied to the conversion of chorismate to prephenate by Chorismate Mutase, where the substrate is treated at the MP2 level of theory while the rest of the system is treated at the RHF level. MP2 geometry optimization is found to lower the barrier by up to 3.5 kcal/mol compared to RHF optimzations and ONIOM energy refinement and leads to a smoother convergence with respect to the basis set for the reaction profile. For double zeta basis sets the increase in CPU time relative to RHF is roughly a factor of two. PMID:24558430
Steinmann, Casper; Fedorov, Dmitri G.; Jensen, Jan H.
2013-01-01
We extend the Effective Fragment Molecular Orbital (EFMO) method to the frozen domain approach where only the geometry of an active part is optimized, while the many-body polarization effects are considered for the whole system. The new approach efficiently mapped out the entire reaction path of chorismate mutase in less than four days using 80 cores on 20 nodes, where the whole system containing 2398 atoms is treated in the ab initio fashion without using any force fields. The reaction path is constructed automatically with the only assumption of defining the reaction coordinate a priori. We determine the reaction barrier of chorismate mutase to be kcal mol−1 for MP2/cc-pVDZ and for MP2/cc-pVTZ in an ONIOM approach using EFMO-RHF/6-31G(d) for the high and low layers, respectively. PMID:23593259
NASA Astrophysics Data System (ADS)
Safouhi, Hassan; Hoggan, Philip
2003-01-01
This review on molecular integrals for large electronic systems (MILES) places the problem of analytical integration over exponential-type orbitals (ETOs) in a historical context. After reference to the pioneering work, particularly by Barnett, Shavitt and Yoshimine, it focuses on recent progress towards rapid and accurate analytic solutions of MILES over ETOs. Software such as the hydrogenlike wavefunction package Alchemy by Yoshimine and collaborators is described. The review focuses on convergence acceleration of these highly oscillatory integrals and in particular it highlights suitable nonlinear transformations. Work by Levin and Sidi is described and applied to MILES. A step by step description of progress in the use of nonlinear transformation methods to obtain efficient codes is provided. The recent approach developed by Safouhi is also presented. The current state of the art in this field is summarized to show that ab initio analytical work over ETOs is now a viable option.
Molecular orbital studies in oxidation: Sulfate formation and metal-metal oxide adhesion
NASA Technical Reports Server (NTRS)
Anderson, A. B.
1985-01-01
The chemical mechanisms for sulfate formation from sodium chloride and sulfur trioxide, which is a product of jet fuel combustion was determined. Molten sodium sulfate leads to hot corrosion of the protective oxide layers on turbine blades. How yttrium dopants in nidkel-aluminum alloys used in turbine blades reduce the spalling rate of protective alumina films and enhance their adhesion was also determined. Two other fulfate mechanisms were deduced and structure of carbon monoxide on a clean chronium and clean platinum-titanium alloys surfaces was determined. All studies were by use of the atom superposition and electron delocalization molecular orbital (ASED-MO) theory. Seven studies were completed. Their titles and abstracts are given.
NASA Astrophysics Data System (ADS)
Chakraborty, Debajit; Karasiev, Valentin; Trickey, Samuel
Aluminum is frequently used in warm-dense matter (WDM) experiments. However, experimental diagnostic limitations make computational exploration of the Al liquid-vapor transition important. The elevated temperaure and low-density make ab initio molecular dynamics (AIMD) with Kohn-Sham (KS) density functional theory (DFT) searches for the divergent compressibility extremely time consuming. Orbital free DFT (OFDFT) in principle is a cost-effective alternative. Here we report on calculations utilizing the PROFESS@QuantumEspresso interface to explore suitable pseudo-potentials, the limitations of our wholly constraint-based VT84F non-ineracting free-energy functional as exposed in the low-density regime, and possible extensions or extrapolations via tunable non-interacting free energy functionals. Work supported by U.S. Dept. of Energy, Grant DE-SC0002139.
Yuan, Wenjuan E-mail: luojunkink@126.com; Yang, Hongping; Zhu, Jing; Luo, Jun E-mail: luojunkink@126.com
2015-11-15
Defects are capable of modulating various properties of graphene, and thus controlling defects is useful in the development of graphene-based devices. Here we present first-principles calculations, which reveal a new avenue for defect engineering of graphene: the modulation by defects on the highest occupied molecular orbital (HOMO) energy of a charged monolayer graphene quantum dot (GQD) is discriminative. When the charge of a GQD increases its HOMO energy also increases. Importantly, when the GQD contains one particular class of defects its HOMO energy is sometimes higher and sometimes lower than that of the corresponding GQD without any defects, but when the GQD contains another class of defects its HOMO energy is always higher or lower than that of the corresponding intact GQD as its excess charge reaches a critical value. This discriminative modulation could allow defect engineering to control secondary electron ejection in graphene, leading to a new way to develop graphene-based devices.
Song, Lingchun; Han, Jaebeom; Lin, Yen-lin; Xie, Wangshen; Gao, Jiali
2010-01-01
The explicit polarization (X-Pol) method has been examined using ab initio molecular orbital theory and density functional theory. The X-Pol potential was designed to provide a novel theoretical framework for developing next-generation force fields for biomolecular simulations. Importantly, the X-Pol potential is a general method, which can be employed with any level of electronic structure theory. The present study illustrates the implementation of the X-Pol method using ab initio Hartree—Fock theory and hybrid density functional theory. The computational results are illustrated by considering a set of bimolecular complexes of small organic molecules and ions with water. The computed interaction energies and hydrogen bond geometries are in good accord with CCSD(T) calculations and B3LYP/aug-cc-pVDZ optimizations. PMID:19618944
Development of the four-body corrected fragment molecular orbital (FMO4) method
NASA Astrophysics Data System (ADS)
Nakano, Tatsuya; Mochizuki, Yuji; Yamashita, Katsumi; Watanabe, Chiduru; Fukuzawa, Kaori; Segawa, Katsunori; Okiyama, Yoshio; Tsukamoto, Takayuki; Tanaka, Shigenori
2012-01-01
The four-body corrected fragment molecular orbital (FMO4) method was implemented at the second-order Møller-Plesset perturbation (MP2) level. A series of accuracy tests relative to the previous two-body and three-body treatments were performed. As expected, FMO4 provided better accuracy in total energies in comparison with the reference values by regular MO calculations. A nonconventional fragmentation by separating main and side chains in amino acid residues was examined for Ala-pentamer and Chignolin, where the four-body corrections were shown to be substantial. A large complex of HIV-1 protease (total 198 residues) with lopinavir was calculated as well. Furthermore, this new FMO scheme was successfully applied to adamantane-shaped clusters with three-dimensional bonding framework.
ERIC Educational Resources Information Center
Halkides, Christopher J.
2013-01-01
In this activity, students manipulate three-dimensional molecular models of the Ala-Ala-Ala tripeptide, where Ala is alanine. They rotate bonds to show that the pairs of dihedral angles phi = 0 degrees, psi = 180 degrees, and phi = 0 degrees, psi = 0 degrees lead to unfavorable interactions among the main chain atoms of the tripeptide. This…
Liu, Shubin; Pedersen, Lee G.
2009-01-01
An effective approach of estimating molecular pKa values from simple density functional calculations is proposed in this work. Both the molecular electrostatic potential (MEP) at the nucleus of the acidic atom and the sum of valence natural atomic orbitals are employed for three categories of compounds, amines and anilines, carbonyl acids and alcohols, and sulfonic acids and thiols. A strong correlation between experimental pKa values and each of these two quantities for each of the three categories has been discovered. Moreover, if the MEP is subtracted by the isolated atomic MEP for each category of compounds, we observe a single unique linear relationship between the resultant MEP difference and experimental pKa data of amines, anilines, carbonyl acids, alcohols, sulfonic acids, thiols, and their substituents. These results can generally be utilized to simultaneously estimate pKa values at multiple sites with a single calculation for either relatively small molecules in drug design or amino acids in proteins and macromolecules. PMID:19317439
NASA Astrophysics Data System (ADS)
Li, Xiao-Hong; Mei, Zheng; Zhang, Xian-Zhou
2014-01-01
The vibrational frequencies of acetophenone thiosemicarbazone in the ground state have been calculated using density functional method (B3LYP) with 6-31G(d), 6-31G(d,p) and 6-311++G(d,p) basis sets. The analysis of natural bond orbital was also performed. The IR spectra were obtained and interpreted by means of potential energies distributions (PEDs) using MOLVIB program. In addition, the results show that there exist Nsbnd H…N and Nsbnd H…S hydrogen bonds in the title compound, which play a major role in stabilizing the molecule and are confirmed by the natural bond orbital analysis. The predicted NLO properties show that the title compound is a good candidate as second-order NLO material. In addition, the frontier molecular orbitals were analyzed and the crystal structure obtained by molecular mechanics belongs to the Pbca space group, with lattice parameters Z = 8, a = 16.0735 Å, b = 7.1719 Å, c = 7.8725 Å, ρ = 0.808 g/cm3.
Li, Xiao-Hong; Mei, Zheng; Zhang, Xian-Zhou
2014-01-24
The vibrational frequencies of acetophenone thiosemicarbazone in the ground state have been calculated using density functional method (B3LYP) with 6-31G(d), 6-31G(d,p) and 6-311++G(d,p) basis sets. The analysis of natural bond orbital was also performed. The IR spectra were obtained and interpreted by means of potential energies distributions (PEDs) using MOLVIB program. In addition, the results show that there exist N-H…N and N-H…S hydrogen bonds in the title compound, which play a major role in stabilizing the molecule and are confirmed by the natural bond orbital analysis. The predicted NLO properties show that the title compound is a good candidate as second-order NLO material. In addition, the frontier molecular orbitals were analyzed and the crystal structure obtained by molecular mechanics belongs to the Pbca space group, with lattice parameters Z=8, a=16.0735 Å, b=7.1719 Å, c=7.8725 Å, ρ=0.808 g/cm(3). PMID:24084483
Chen, N.; Yang, R.T.
1996-11-01
An ab initio molecular orbital study is undertaken on the adsorption of N{sub 2}, O{sub 2}, and C{sub 2}H{sub 4} (adsorbate) on Ag-zeolite and Ag halides (adsorbent). Geometry optimization is performed at the HF/3-21G level, while MP2/3-21G with natural bond orbital calculations are performed to obtain energies, atomic charges, orbital energies, and orbital populations (occupancies). The bonding of adsorbate to adsorbent is discussed in the context of {sigma}-donation (i.e., overlap of the 2p orbitals of the adsorbate molecule with the 5s orbital of Ag) and d-{pi}* back donation (i.e., overlap of the 4d{sub yz} orbitals of Ag with the 2p* antibonding orbitals of the adsorbate). For all adsorbate-adsorbent pairs, the ratio of {sigma}-donation to d-{pi}* back donation is approximately 3:1. Results on occupancy analysis indicate that a considerable electron redistribution from the 4d{sub zy} orbitals to the 4d{sub yz} orbitals occurs in Ag during adsorption and that this redistribution has possibly enhanced the d-{pi}* back donation. Net charge and energy gap ({Delta}{epsilon}) analyses indicate that it is slightly easier for N{sub 2} than O{sub 2} to adsorb, whereas a comparison of N{sub 2} and O{sub 2} adsorption from calculations of the energies of adsorption is inconclusive. However, a fair agreement is obtained in comparison of theory and experiment for energy of adsorption of N{sub 2} and C{sub 2}J{sub 4} on Ag-zeolite. The dispersion energies of adsorption, based on the MP2 correlation energies, are nearly the same for all adsorption pairs, i.e,, approximately 4--5 kcal/mol.
NASA Astrophysics Data System (ADS)
Habeeb, Moustafa M.; Al-Attas, Amirah S.; Al-Raimi, Doaa S.
2015-05-01
Charge transfer (CT) interaction between 3,5-dimethylpyrazole (DMP) with the π-acceptor 2,3-dichloro-5,6-dicyano-p-benzoquinon (DDQ) has been investigated spectrophotometrically in acetonitrile (AN). Simultaneous reddish brown color has been observed upon mixing donor with acceptor solutions attributing to CT complex formation. The electronic spectra of the formed complex exhibited multi-charge transfer bands at 429, 447, 506, 542 and 589 nm, respectively. Job's method of continuous variations and spectrophotometric titration methods confirmed the formation of the studied complex in 1:2 ratio between DMP and DDQ. Benesi-Hildebrand equation has been applied to calculate the stability constant of the formed complex where it recorded high value supporting formation of stable complex. Molecular orbital calculations using MM2 method and GAMESS (General Atomic and Molecular Electronic Structure System) interface computations as a package of ChemBio3D Ultra12 software were carried out for more analysis of the formed complex in the gas phase. The computational analysis included energy minimisation, stabilisation energy, molecular geometry, Mullikan charges, molecular electrostatic potential (MEP) surfaces of reactants and complex as well as characterization of the higher occupied molecular orbitals (HOMO) and lower unoccupied molecular orbitals (LUMO) surfaces of the complex. A good consistency between experimental and theoretical results has been recorded.
NASA Technical Reports Server (NTRS)
Almloef, Jan; Taylor, Peter R.
1989-01-01
A recently proposed scheme for using natural orbitals from atomic configuration interaction (CI) wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outmost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital (ANO) sets.
NASA Technical Reports Server (NTRS)
Almlof, Jan; Taylor, Peter R.
1990-01-01
A recently proposed scheme for using natural orbitals from atomic configuration interaction wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outermost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital sets.
ERIC Educational Resources Information Center
Litofsky, Joshua; Viswanathan, Rama
2015-01-01
Matrix diagonalization, the key technique at the heart of modern computational chemistry for the numerical solution of the Schrödinger equation, can be easily introduced in the physical chemistry curriculum in a pedagogical context using simple Hückel molecular orbital theory for p bonding in molecules. We present details and results of…
Saito, N; Fuwa, A
2000-01-01
In this investigation, respective thermodynamic parameters of heats of formation, standard entropy and specific heat capacity at constant pressure for PCDDs, PCDFs, Co-PCB and PCBs as well as polychlorinated-benzenes and polychlorinated-phenols have been evaluated by quantum chemical calculation using a semi-empirical molecular orbital method with the PM3 Hamiltonian and statistical thermodynamic correlation. PMID:10665426
Stan, Raluca-Maria; Gaina, Roxana; Enachescu, Cristian E-mail: radu.tanasa@uaic.ro; Stancu, Alexandru; Tanasa, Radu E-mail: radu.tanasa@uaic.ro; Bronisz, Robert
2015-05-07
In this paper, we analyze two types of hysteresis in spin crossover molecular magnets compounds in the framework of the First Order Reversal Curve (FORC) method. The switching between the two stable states in these compounds is accompanied by hysteresis phenomena if the intermolecular interactions are higher than a threshold. We have measured the static thermal hysteresis (TH) and the kinetic light induced thermal hysteresis (LITH) major loops and FORCs for the polycrystalline Fe(II) spin crossover compound [Fe{sub 1−x}Zn{sub x}(bbtr){sub 3}](ClO{sub 4}){sub 2} (bbtr = 1,4-di(1,2,3-triazol-1-yl)butane), either in a pure state (x = 0) or doped with Zn ions (x = 0.33) considering different sweeping rates. Here, we use this method not only to infer the domains distribution but also to disentangle between kinetic and static components of the LITH and to estimate the changes in the intermolecular interactions introduced by dopants. We also determined the qualitative relationship between FORC distributions measured for TH and LITH.
NASA Astrophysics Data System (ADS)
Wyrick, Jonathan; Einstein, T. L.; Bartels, Ludwig
2015-03-01
We present a method of analyzing the results of density functional modeling of molecular adsorption in terms of an analogue of molecular orbitals. This approach permits intuitive chemical insight into the adsorption process. Applied to a set of anthracene derivates (anthracene, 9,10-anthraquinone, 9,10-dithioanthracene, and 9,10-diselenonanthracene), we follow the electronic states of the molecules that are involved in the bonding process and correlate them to both the molecular adsorption geometry and the species' diffusive behavior. We additionally provide computational code to easily repeat this analysis on any system.
Wyrick, Jonathan; Bartels, Ludwig; Einstein, T. L.
2015-03-14
We present a method of analyzing the results of density functional modeling of molecular adsorption in terms of an analogue of molecular orbitals. This approach permits intuitive chemical insight into the adsorption process. Applied to a set of anthracene derivates (anthracene, 9,10-anthraquinone, 9,10-dithioanthracene, and 9,10-diselenonanthracene), we follow the electronic states of the molecules that are involved in the bonding process and correlate them to both the molecular adsorption geometry and the species’ diffusive behavior. We additionally provide computational code to easily repeat this analysis on any system.
Analytic Gradient for Density Functional Theory Based on the Fragment Molecular Orbital Method.
Brorsen, Kurt R; Zahariev, Federico; Nakata, Hiroya; Fedorov, Dmitri G; Gordon, Mark S
2014-12-01
The equations for the response terms for the fragment molecular orbital (FMO) method interfaced with the density functional theory (DFT) gradient are derived and implemented. Compared to the previous FMO-DFT gradient, which lacks response terms, the FMO-DFT analytic gradient has improved accuracy for a variety of functionals, when compared to numerical gradients. The FMO-DFT gradient agrees with the fully ab initio DFT gradient in which no fragmentation is performed, while reducing the nonlinear scaling associated with standard DFT. Solving for the response terms requires the solution of the coupled perturbed Kohn-Sham (CPKS) equations, where the CPKS equations are solved through a decoupled Z-vector procedure called the self-consistent Z-vector method. FMO-DFT is a nonvariational method and the FMO-DFT gradient is unique compared to standard DFT gradients in that the FMO-DFT gradient requires terms from both DFT and time-dependent density functional theory (TDDFT) theories. PMID:26583213
Second order Møller-Plesset perturbation theory based upon the fragment molecular orbital method
NASA Astrophysics Data System (ADS)
Fedorov, Dmitri G.; Kitaura, Kazuo
2004-08-01
The fragment molecular orbital (FMO) method was combined with the second order Møller-Plesset (MP2) perturbation theory. The accuracy of the method using the 6-31G* basis set was tested on (H2O)n, n=16,32,64; α-helices and β-strands of alanine n-mers, n=10,20,40; as well as on (H2O)n, n=16,32,64 using the 6-31++G** basis set. Relative to the regular MP2 results that could be afforded, the FMO2-MP2 error in the correlation energy did not exceed 0.003 a.u., the error in the correlation energy gradient did not exceed 0.000 05 a.u./bohr and the error in the correlation contribution to dipole moment did not exceed 0.03 debye. An approximation reducing computational load based on fragment separation was introduced and tested. The FMO2-MP2 method demonstrated nearly linear scaling and drastically reduced the memory requirements of the regular MP2, making possible calculations with several thousands basis functions using small Pentium clusters. As an example, (H2O)64 with the 6-31++G** basis set (1920 basis functions) can be run in 1 Gbyte RAM and it took 136 s on a 40-node Pentium4 cluster.
NASA Astrophysics Data System (ADS)
Kubicki, J. D.; Heaney, P. J.
2002-12-01
Researchers recently have made the discovery that hypercoordinate Si-sorbitol complexes will readily form in biologically relevant fluids, and they have reported the first evidence for a transient organosilicon complex generated within the life cycle of an organism. These interpretations are based upon peak assignments of Si-29 NMR spectra that invoke Si-polyol complexes with Si in five- and six-fold coordination states. However, ab initio analyses of the proposed organosilicon structures do not reproduce the experimentally observed chemical shifts. We have successfully modeled one of the observed Si-29 chemical shifts with a 5-fold Si-disorbitol complex involving 5-membered ring configurations (i.e., Si-O-C-C-O), which yielded Si-29 chemical shifts that closely matched the observed values in the -100 to -102 ppm range. Likewise, Si-29 NMR peaks near -144 ppm were well fit by a model in which a 6-fold Si was complexed to three sorbitol molecules in a 5-membered ring configuration. The ability to simulate observed NMR signals using molecular orbital calculations provides strong support for the controversial role of hypercoordinate organosilicon species in the uptake and transport of silica by biological systems. The existence of such complexes in turn may explain other puzzles in Si biogeochemistry, such as the persistence of monomeric silica in concentrated biological fluids and the biofractionation of Si isotopes and Ge.
Auxiliary functions for molecular integrals with Slater-type orbitals. II. Gauss transform methods
NASA Astrophysics Data System (ADS)
Ema, I.; López, R.; Fernández, J. J.; Ramírez, G.; Rico, J. F.
The Gauss transform of Slater-type orbitals is used to express several types of molecular integrals involving these functions in terms of simple auxiliary functions. After reviewing this transform and the way it can be combined with the shift operator technique, a master formula for overlap integrals is derived and used to obtain multipolar moments associated to fragments of two-center distributions and overlaps of derivatives of Slater functions. Moreover, it is proved that integrals involving two-center distributions and irregular harmonics placed at arbitrary points (which determine the electrostatic potential, field and field gradient, as well as higher order derivatives of the potential) can be expressed in terms of auxiliary functions of the same type as those appearing in the overlap. The recurrence relations and series expansions of these functions are thoroughly studied, and algorithms for their calculation are presented. The usefulness and efficiency of this procedure are tested by developing two independent codes: one for the derivatives of the overlap integrals with respect to the centers of the functions, and another for derivatives of the potential (electrostatic field, field gradient, and so forth) at arbitrary points.0
Krah, Tim; Ben Amor, Nadia; Maynau, Daniel; Berger, J A; Robert, Vincent
2014-07-01
Based on localized molecular orbitals, the proposed method reduces large configuration interaction (CI) spaces while maintaining agreement with reference values. Our strategy concentrates the numerical effort on physically pertinent CI-contributions and is to be considered as a tool to tackle large systems including numerous open-shells. To show the efficiency of our method we consider two 4-electron parent systems. First, we illustrate our approach by describing the van der Waals interactions in the (H2)2 system. By systematically including local correlation, dispersion and charge transfer mechanisms, we show that 90% of the reference full CI dissociation energy of the H2 dimer is reproduced using only 3% of the full CI space. Second, the conformational cis/trans rotation barrier of the butadiene molecule is remarkably reproduced (97% of the reference value) with less than 1% of the reference space. This work paves the way to numerical strategies which afford the electronic structure determination of large open-shell systems avoiding the exponential limitation. At the same time, a physical analysis of the contents of the wave function is offered. PMID:24935105
Okumura, Hironori McSkimming, Brian M.; Speck, James S.; Huault, Thomas; Chaix, Catherine
2014-01-06
N-face GaN was grown on free-standing GaN (0001{sup ¯}) substrates at a growth rate of 1.5 μm/h using plasma-assisted molecular beam epitaxy. Difference in growth rate between (0001{sup ¯}) and (0001) oriented GaN depends on nitrogen plasma power, and the (0001{sup ¯}) oriented GaN had only 70% of the growth rate of the (0001) oriented GaN at 300 W. Unintentional impurity concentrations of silicon, carbon, and oxygen were 2 × 10{sup 15}, 2 × 10{sup 16}, and 7 × 10{sup 16} cm{sup −3}, respectively. A growth diagram was constructed that shows the dependence of the growth modes on the difference in the Ga and active nitrogen flux, Φ{sub Ga} − Φ{sub N*}, and the growth temperature. At high Φ{sub Ga} − Φ{sub N*} (Φ{sub Ga} ≫ Φ{sub N*}), two-dimensional (step-flow and layer-by-layer) growth modes were realized. High growth temperature (780 °C) expanded the growth window of the two-dimensional growth modes, achieving a surface with rms roughness of 0.48 nm without Ga droplets.
The Space Shuttle Orbiter molecular environment induced by the supplemental flash evaporator system
NASA Technical Reports Server (NTRS)
Ehlers, H. K. F.
1985-01-01
The water vapor environment of the Space Shuttle Orbiter induced by the supplemental flash evaporator during the on-orbit flight phase has been analyzed based on Space II model predictions and orbital flight measurements. Model data of local density, column density, and return flux are presented. Results of return flux measurements with a mass spectrometer during STS-2 and of direct flux measurements during STS-4 are discussed and compared with model predictions.
Findlater, Alexander D; Zahariev, Federico; Gordon, Mark S
2015-04-16
The local correlation "cluster-in-molecule" (CIM) method is combined with the fragment molecular orbital (FMO) method, providing a flexible, massively parallel, and near-linear scaling approach to the calculation of electron correlation energies for large molecular systems. Although the computational scaling of the CIM algorithm is already formally linear, previous knowledge of the Hartree-Fock (HF) reference wave function and subsequent localized orbitals is required; therefore, extending the CIM method to arbitrarily large systems requires the aid of low-scaling/linear-scaling approaches to HF and orbital localization. Through fragmentation, the combined FMO-CIM method linearizes the scaling, with respect to system size, of the HF reference and orbital localization calculations, achieving near-linear scaling at both the reference and electron correlation levels. For the 20-residue alanine α helix, the preliminary implementation of the FMO-CIM method captures 99.6% of the MP2 correlation energy, requiring 21% of the MP2 wall time. The new method is also applied to solvated adamantine to illustrate the multilevel capability of the FMO-CIM method. PMID:25794346
NASA Astrophysics Data System (ADS)
Philipp, Dean Michael
Methodology is discussed for mixed ab initio quantum mechanics/molecular mechanics modeling of systems where the quantum mechanics (QM) and molecular mechanics (MM) regions are within the same molecule. The ab initio QM calculations are at the restricted Hartree-Fock level using the pseudospectral method of the Jaguar program while the MM part is treated with the OPLS force fields implemented in the IMPACT program. The interface between the QM and MM regions, in particular, is elaborated upon, as it is dealt with by ``breaking'' bonds at the boundaries and using Boys-localized orbitals found from model molecules in place of the bonds. These orbitals are kept frozen during QM calculations. The mixed modeling presented here can be used for single point energy calculations and geometry optimizations. Results from tests of the method to find relative conformational energies and geometries of alanine tetrapeptides are presented along with comparisons to pure QM and pure MM calculations.
NASA Astrophysics Data System (ADS)
Drake, Kyle; Bonacum, Jason; Zhang, Guo-Ping
2014-03-01
The molecular structure of Buckminster fullerene (C60) allows for electron delocalization in all of the pi-bonding electrons of the molecule. This coupled with the symmetry of the molecule allows for the formation of super-atomic molecular orbitals (SAMOs) similar to those observed in aluminum clusters. The SAMOs behave as if the molecule that they belong to is a single atom. We compute the eigenstates of C60 compulationally using density functional theory (DFT) and a grid mesh. Using larger radii also allows us to accurately describe SAMOs and test the convergence of our data. The results are interesting because for the first time, we can show the true converged super atomic orbitals in C60. Indiana State University SURE Program, Department of Energy, Indiana State University Department of Physics, and Indiana State University Center for Student Creativity and Research.
NASA Astrophysics Data System (ADS)
Karaca, Caglar; Atac, Ahmet; Karabacak, Mehmet
2015-02-01
In this study, 5-iodosalicylic acid (5-ISA, C7H5IO3) is structurally characterized by FT-IR, FT-Raman, NMR and UV spectroscopies. There are eight conformers, Cn, n = 1-8 for this molecule therefore the molecular geometry for these eight conformers in the ground state are calculated by using the ab-initio density functional theory (DFT) B3LYP method approach with the aug-cc-pVDZ-PP basis set for iodine and the aug-cc-pVDZ basis set for the other elements. The computational results identified that the most stable conformer of 5-ISA is the C1 form. The vibrational spectra are calculated DFT method invoking the same basis sets and fundamental vibrations are assigned on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method with PQS program. Total density of state (TDOS) and partial density of state (PDOS) and also overlap population density of state (COOP or OPDOS) diagrams analysis for C1 conformer were calculated using the same method. The energy and oscillator strength are calculated by time-dependent density functional theory (TD-DFT) results complement with the experimental findings. Besides, charge transfer occurring in the molecule between HOMO and LUMO energies, frontier energy gap, molecular electrostatic potential (MEP) are calculated and presented. The NMR chemical shifts (1H and 13C) spectra are recorded and calculated using the gauge independent atomic orbital (GIAO) method. Mulliken atomic charges of the title molecule are also calculated, interpreted and compared with salicylic acid. The optimized bond lengths, bond angles and calculated NMR and UV, vibrational wavenumbers showed the best agreement with the experimental results.
Pharmacophore Modeling for Anti-Chagas Drug Design Using the Fragment Molecular Orbital Method
Ohno, Kazuki; Orita, Masaya; Inoue, Masayuki; Shiba, Tomoo; Harada, Shigeharu; Honma, Teruki; Balogun, Emmanuel Oluwadare; da Rocha, Josmar Rodrigues; Montanari, Carlos Alberto; Kita, Kiyoshi; Sekijima, Masakazu
2015-01-01
Background Chagas disease, caused by the parasite Trypanosoma cruzi, is a neglected tropical disease that causes severe human health problems. To develop a new chemotherapeutic agent for the treatment of Chagas disease, we predicted a pharmacophore model for T. cruzi dihydroorotate dehydrogenase (TcDHODH) by fragment molecular orbital (FMO) calculation for orotate, oxonate, and 43 orotate derivatives. Methodology/Principal Findings Intermolecular interactions in the complexes of TcDHODH with orotate, oxonate, and 43 orotate derivatives were analyzed by FMO calculation at the MP2/6-31G level. The results indicated that the orotate moiety, which is the base fragment of these compounds, interacts with the Lys43, Asn67, and Asn194 residues of TcDHODH and the cofactor flavin mononucleotide (FMN), whereas functional groups introduced at the orotate 5-position strongly interact with the Lys214 residue. Conclusions/Significance FMO-based interaction energy analyses revealed a pharmacophore model for TcDHODH inhibitor. Hydrogen bond acceptor pharmacophores correspond to Lys43 and Lys214, hydrogen bond donor and acceptor pharmacophores correspond to Asn67 and Asn194, and the aromatic ring pharmacophore corresponds to FMN, which shows important characteristics of compounds that inhibit TcDHODH. In addition, the Lys214 residue is not conserved between TcDHODH and human DHODH. Our analysis suggests that these orotate derivatives should preferentially bind to TcDHODH, increasing their selectivity. Our results obtained by pharmacophore modeling provides insight into the structural requirements for the design of TcDHODH inhibitors and their development as new anti-Chagas drugs. PMID:25961853
Using the fragment molecular orbital method to investigate agonist-orexin-2 receptor interactions.
Heifetz, Alexander; Aldeghi, Matteo; Chudyk, Ewa I; Fedorov, Dmitri G; Bodkin, Mike J; Biggin, Philip C
2016-04-15
The understanding of binding interactions between any protein and a small molecule plays a key role in the rationalization of affinity and selectivity and is essential for an efficient structure-based drug discovery (SBDD) process. Clearly, to begin SBDD, a structure is needed, and although there has been fantastic progress in solving G-protein-coupled receptor (GPCR) crystal structures, the process remains quite slow and is not currently feasible for every GPCR or GPCR-ligand complex. This situation significantly limits the ability of X-ray crystallography to impact the drug discovery process for GPCR targets in 'real-time' and hence there is still a need for other practical and cost-efficient alternatives. We present here an approach that integrates our previously described hierarchical GPCR modelling protocol (HGMP) and the fragment molecular orbital (FMO) quantum mechanics (QM) method to explore the interactions and selectivity of the human orexin-2 receptor (OX2R) and its recently discovered nonpeptidic agonists. HGMP generates a 3D model of GPCR structures and its complexes with small molecules by applying a set of computational methods. FMO allowsab initioapproaches to be applied to systems that conventional QM methods would find challenging. The key advantage of FMO is that it can reveal information on the individual contribution and chemical nature of each residue and water molecule to the ligand binding that normally would be difficult to detect without QM. We illustrate how the combination of both techniques provides a practical and efficient approach that can be used to analyse the existing structure-function relationships (SAR) and to drive forward SBDD in a real-world example for which there is no crystal structure of the complex available. PMID:27068972
Roemelt, Michael
2015-07-28
Spin Orbit Coupling (SOC) is introduced to molecular ab initio density matrix renormalization group (DMRG) calculations. In the presented scheme, one first approximates the electronic ground state and a number of excited states of the Born-Oppenheimer (BO) Hamiltonian with the aid of the DMRG algorithm. Owing to the spin-adaptation of the algorithm, the total spin S is a good quantum number for these states. After the non-relativistic DMRG calculation is finished, all magnetic sublevels of the calculated states are constructed explicitly, and the SOC operator is expanded in the resulting basis. To this end, spin orbit coupled energies and wavefunctions are obtained as eigenvalues and eigenfunctions of the full Hamiltonian matrix which is composed of the SOC operator matrix and the BO Hamiltonian matrix. This treatment corresponds to a quasi-degenerate perturbation theory approach and can be regarded as the molecular equivalent to atomic Russell-Saunders coupling. For the evaluation of SOC matrix elements, the full Breit-Pauli SOC Hamiltonian is approximated by the widely used spin-orbit mean field operator. This operator allows for an efficient use of the second quantized triplet replacement operators that are readily generated during the non-relativistic DMRG algorithm, together with the Wigner-Eckart theorem. With a set of spin-orbit coupled wavefunctions at hand, the molecular g-tensors are calculated following the scheme proposed by Gerloch and McMeeking. It interprets the effective molecular g-values as the slope of the energy difference between the lowest Kramers pair with respect to the strength of the applied magnetic field. Test calculations on a chemically relevant Mo complex demonstrate the capabilities of the presented method.
NASA Astrophysics Data System (ADS)
Roemelt, Michael
2015-07-01
Spin Orbit Coupling (SOC) is introduced to molecular ab initio density matrix renormalization group (DMRG) calculations. In the presented scheme, one first approximates the electronic ground state and a number of excited states of the Born-Oppenheimer (BO) Hamiltonian with the aid of the DMRG algorithm. Owing to the spin-adaptation of the algorithm, the total spin S is a good quantum number for these states. After the non-relativistic DMRG calculation is finished, all magnetic sublevels of the calculated states are constructed explicitly, and the SOC operator is expanded in the resulting basis. To this end, spin orbit coupled energies and wavefunctions are obtained as eigenvalues and eigenfunctions of the full Hamiltonian matrix which is composed of the SOC operator matrix and the BO Hamiltonian matrix. This treatment corresponds to a quasi-degenerate perturbation theory approach and can be regarded as the molecular equivalent to atomic Russell-Saunders coupling. For the evaluation of SOC matrix elements, the full Breit-Pauli SOC Hamiltonian is approximated by the widely used spin-orbit mean field operator. This operator allows for an efficient use of the second quantized triplet replacement operators that are readily generated during the non-relativistic DMRG algorithm, together with the Wigner-Eckart theorem. With a set of spin-orbit coupled wavefunctions at hand, the molecular g-tensors are calculated following the scheme proposed by Gerloch and McMeeking. It interprets the effective molecular g-values as the slope of the energy difference between the lowest Kramers pair with respect to the strength of the applied magnetic field. Test calculations on a chemically relevant Mo complex demonstrate the capabilities of the presented method.
Roemelt, Michael
2015-07-28
Spin Orbit Coupling (SOC) is introduced to molecular ab initio density matrix renormalization group (DMRG) calculations. In the presented scheme, one first approximates the electronic ground state and a number of excited states of the Born-Oppenheimer (BO) Hamiltonian with the aid of the DMRG algorithm. Owing to the spin-adaptation of the algorithm, the total spin S is a good quantum number for these states. After the non-relativistic DMRG calculation is finished, all magnetic sublevels of the calculated states are constructed explicitly, and the SOC operator is expanded in the resulting basis. To this end, spin orbit coupled energies and wavefunctions are obtained as eigenvalues and eigenfunctions of the full Hamiltonian matrix which is composed of the SOC operator matrix and the BO Hamiltonian matrix. This treatment corresponds to a quasi-degenerate perturbation theory approach and can be regarded as the molecular equivalent to atomic Russell-Saunders coupling. For the evaluation of SOC matrix elements, the full Breit-Pauli SOC Hamiltonian is approximated by the widely used spin-orbit mean field operator. This operator allows for an efficient use of the second quantized triplet replacement operators that are readily generated during the non-relativistic DMRG algorithm, together with the Wigner-Eckart theorem. With a set of spin-orbit coupled wavefunctions at hand, the molecular g-tensors are calculated following the scheme proposed by Gerloch and McMeeking. It interprets the effective molecular g-values as the slope of the energy difference between the lowest Kramers pair with respect to the strength of the applied magnetic field. Test calculations on a chemically relevant Mo complex demonstrate the capabilities of the presented method. PMID:26233112
ERIC Educational Resources Information Center
Bridgeman, Adam J.; Schmidt, Timothy W.; Young, Nigel A.
2013-01-01
The stretching modes of ML[subscript "x"] complexes have the same symmetry as the atomic orbitals on M that are used to form its s bonds. In the exercise suggested here, the atomic orbitals are used to derive the form of the stretching modes without the need for formal group theory. The analogy allows students to help understand many…
NASA Astrophysics Data System (ADS)
Seo, Won-Gap; Matsuura, Hiroyuki; Tsukihashi, Fumitaka
2006-04-01
Recently, molecular dynamics (MD) simulation has been widely employed as a very useful method for the calculation of various physicochemical properties in the molten slags and fluxes. In this study, MD simulation has been applied to calculate the structural, transport, and thermodynamic properties for the FeCl2, PbCl2, and ZnCl2 systems using the Born—Mayer—Huggins type pairwise potential with partial ionic charges. The interatomic potential parameters were determined by fitting the physicochemical properties of iron chloride, lead chloride, and zinc chloride systems with experimentally measured results. The calculated structural, transport, and thermodynamic properties of pure FeCl2, PbCl2, and ZnCl2 showed the same tendency with observed results. Especially, the calculated structural properties of molten ZnCl2 and FeCl2 show the possibility of formation of polymeric network structures based on the ionic complexes of ZnCl{4/2-}, ZnCl{3/-}, FeCl{4/2-}, and FeCl{3/-}, and these calculations have successfully reproduced the measured results. The enthalpy, entropy, and Gibbs energy of mixing for the PbCl2-ZnCl2, FeCl2-PbCl2, and FeCl2-ZnCl2 systems were calculated based on the thermodynamic and structural parameters of each binary system obtained from MD simulation. The phase diagrams of the PbCl2-ZnCl2, FeCl2-PbCl2, and FeCl2-ZnCl2 systems estimated by using the calculated Gibbs energy of mixing reproduced the experimentally measured ones reasonably well.
Pinjari, Rahul V.; Delcey, Mickaël G.; Guo, Meiyuan; Lundberg, Marcus; Odelius, Michael
2014-09-28
The metal L-edge (2p → 3d) X-ray absorption spectra are affected by a number of different interactions: electron-electron repulsion, spin-orbit coupling, and charge transfer between metal and ligands, which makes the simulation of spectra challenging. The core restricted active space (RAS) method is an accurate and flexible approach that can be used to calculate X-ray spectra of a wide range of medium-sized systems without any symmetry constraints. Here, the applicability of the method is tested in detail by simulating three ferric (3d{sup 5}) model systems with well-known electronic structure, viz., atomic Fe{sup 3+}, high-spin [FeCl{sub 6}]{sup 3−} with ligand donor bonding, and low-spin [Fe(CN){sub 6}]{sup 3−} that also has metal backbonding. For these systems, the performance of the core RAS method, which does not require any system-dependent parameters, is comparable to that of the commonly used semi-empirical charge-transfer multiplet model. It handles orbitally degenerate ground states, accurately describes metal-ligand interactions, and includes both single and multiple excitations. The results are sensitive to the choice of orbitals in the active space and this sensitivity can be used to assign spectral features. A method has also been developed to analyze the calculated X-ray spectra using a chemically intuitive molecular orbital picture.
Phipps, M J S; Fox, T; Tautermann, C S; Skylaris, C-K
2016-07-12
We report the development and implementation of an energy decomposition analysis (EDA) scheme in the ONETEP linear-scaling electronic structure package. Our approach is hybrid as it combines the localized molecular orbital EDA (Su, P.; Li, H. J. Chem. Phys., 2009, 131, 014102) and the absolutely localized molecular orbital EDA (Khaliullin, R. Z.; et al. J. Phys. Chem. A, 2007, 111, 8753-8765) to partition the intermolecular interaction energy into chemically distinct components (electrostatic, exchange, correlation, Pauli repulsion, polarization, and charge transfer). Limitations shared in EDA approaches such as the issue of basis set dependence in polarization and charge transfer are discussed, and a remedy to this problem is proposed that exploits the strictly localized property of the ONETEP orbitals. Our method is validated on a range of complexes with interactions relevant to drug design. We demonstrate the capabilities for large-scale calculations with our approach on complexes of thrombin with an inhibitor comprised of up to 4975 atoms. Given the capability of ONETEP for large-scale calculations, such as on entire proteins, we expect that our EDA scheme can be applied in a large range of biomolecular problems, especially in the context of drug design. PMID:27248370
NASA Astrophysics Data System (ADS)
Chen, Mohan; Xia, Junchao; Huang, Chen; Dieterich, Johannes M.; Hung, Linda; Shin, Ilgyou; Carter, Emily A.
2015-05-01
Orbital-free density functional theory (OFDFT) is a linear-scaling first-principles quantum mechanics method used to calculate the ground-state energy of a given system. Here we present a new version of PRinceton Orbital-Free Electronic Structure Software (PROFESS) with new features. First, PROFESS 3.0 provides a set of new kinetic energy density functionals (KEDFs) which are designed to model semiconductors or transition metals. Specifically, PROFESS 3.0 includes the Huang-Carter (HC) KEDF [1], a density decomposition method with fixed localized electronic density [2], the Wang-Govind-Carter (WGC) decomposition KEDF [3], and the Enhanced von Weizsäcker (EvW)-WGC KEDF [4]. Other major new functions are included, such as molecular dynamics with different statistical mechanical ensembles and spin-polarized density optimizers.
Liquid Water through Density-Functional Molecular Dynamics: Plane-Wave vs Atomic-Orbital Basis Sets.
Miceli, Giacomo; Hutter, Jürg; Pasquarello, Alfredo
2016-08-01
We determine and compare structural, dynamical, and electronic properties of liquid water at near ambient conditions through density-functional molecular dynamics simulations, when using either plane-wave or atomic-orbital basis sets. In both frameworks, the electronic structure and the atomic forces are self-consistently determined within the same theoretical scheme based on a nonlocal density functional accounting for van der Waals interactions. The overall properties of liquid water achieved within the two frameworks are in excellent agreement with each other. Thus, our study supports that implementations with plane-wave or atomic-orbital basis sets yield equivalent results and can be used indiscriminately in study of liquid water or aqueous solutions. PMID:27434607
NASA Technical Reports Server (NTRS)
Roberts, W. W., Jr.; Stewart, G. R.
1987-01-01
The different roles played by orbital dynamics and dissipative cloud-cloud collisions in the formation of giant molecular clouds (GMCs) in a global spiral structure are investigated. The interstellar medium (ISM) is simulated by a system of particles, representing clouds, which orbit in a spiral-perturbed, galactic gravitational field. The overall magnitude and width of the global cloud density distribution in spiral arms is very similar in the collisional and collisionless simulations. The results suggest that the assumed number density and size distribution of clouds and the details of individual cloud-cloud collisions have relatively little effect on these features. Dissipative cloud-cloud collisions play an important steadying role for the cloud system's global spiral structure. Dissipative cloud-cloud collisions also damp the relative velocity dispersion of clouds in massive associations and thereby aid in the effective assembling of GMC-like complexes.
Balachandran, V; Karunakaran, V
2014-06-01
The FT-IR (4000-400cm(-)(1)) and FT-Raman spectra (3500-100cm(-)(1)) of 2,3,4,5,6-pentafluorophenylacetic acid (PAA) have been recorded. Density functional theory calculation with LSDA/6-31+G(d,p) and B3LYP/6-31+G(d,p) basis sets have been used to determine ground state molecular geometries (bond lengths and bond angles), harmonic vibrational frequencies, infrared intensities, Raman intensities and bonding features of the title compound. The assignments of the vibrational spectra have been carried out with the help of normal coordinate analysis (NCA) following the scaled quantum mechanical force field (SQMFF) methodology. The first order hyperpolarizability (β0) and related properties (β, α0 and Δα) of PAA are calculated using B3LYP/6-31+G(d,p) method on the finite-field approach. The calculated first hyperpolarizability shows that the molecule is an attractive molecule for future applications in non-linear optics. The stability of molecule has been analyzed by using NBO analysis. The calculated HOMO and LUMO energies show that charge transfer occurs within this molecule. Mulliken population analysis on atomic charges is also calculated. Thermodynamic properties (heat capacity, enthalpy, Gibb's free energy and entropy) of the title compound at different temperatures were calculated. PMID:24662720
NASA Astrophysics Data System (ADS)
Kalfaoǧlu, Emel; Karabulut, Bünyamin
2016-03-01
Electron paramagnetic resonance (EPR) and optical absorption spectra of Cu2+ ions in cesium hydrogen oxalate single crystals have been investigated at room temperature. The spin-Hamiltonian parameters (g and A), have been determined. Crystalline field around the Cu2+ ion is almost axially symmetric. The results show a single paramagnetic site which confirms the triclinic crystal symmetry. Molecular orbital bonding coefficients are studied from the EPR and optical data. Theoretical octahedral field parameter and the tetragonal field parameters have been evaluated from the superposition model. Using these parameters, various bonding parameters are analyzed and the nature of bonding in the complex is discussed. The theoretical results are supported by experimental results.
Molecular orbital (SCF-Xα-SW) theory of metal-metal charge transfer processes in minerals
Sherman, David M.
1987-01-01
Electronic transitions between the Fe-Fe bonding and Fe-Fe antibonding orbitals results in the optically-induced intervalence charge transfer bands observed in the electronic spectra of mixed valence minerals. Such transitions are predicted to be polarized along the metal-metal bond direction, in agreement with experimental observations.
ERIC Educational Resources Information Center
Gillespie, Ronald J.; And Others
1996-01-01
Presents an alternative approach to bonding and geometry--the electron domain model--which avoids some of the problems with the conventional approach. Discusses difficulties with the orbital model at the introductory level, electron spin and the Pauli exclusion principle, electron pair domains, nonequivalent domains, multiple bonds, and origins…
Phase Equilibria Diagrams Database
National Institute of Standards and Technology Data Gateway
SRD 31 NIST/ACerS Phase Equilibria Diagrams Database (PC database for purchase) The Phase Equilibria Diagrams Database contains commentaries and more than 21,000 diagrams for non-organic systems, including those published in all 21 hard-copy volumes produced as part of the ACerS-NIST Phase Equilibria Diagrams Program (formerly titled Phase Diagrams for Ceramists): Volumes I through XIV (blue books); Annuals 91, 92, 93; High Tc Superconductors I & II; Zirconium & Zirconia Systems; and Electronic Ceramics I. Materials covered include oxides as well as non-oxide systems such as chalcogenides and pnictides, phosphates, salt systems, and mixed systems of these classes.
Martin, Daniel J; McCarthy, Brian D; Rountree, Eric S; Dempsey, Jillian L
2016-06-14
The EC' Zone Diagram, introduced by Savéant and Su over 30 years ago, has been used to classify voltammetric responses for electrocatalytic systems. With a single H2-evolving catalyst, Co(dmgBF2)2(CH3CH)2 (dmgBF2 = difluoroboryl-dimethylglyoxime), and a series of para-substituted anilinium acids, experimental conditions were carefully tuned to access to each region of the classic zone diagram. Close scrutiny revealed the extent to which the kinetic (λ) and excess (γ) factors could be experimentally controlled and used to access a variety of waveforms for this ECEC' catalytic system. It was found that most of the tunable experimental parameters (such as catalyst concentration, scan rate, and substrate concentration) predicted in the EC' Zone Diagram could be extended to a multi-electron system and produced similarly-shaped waveforms with some deviations. Tuning of a single catalyst across every region of the classic zone diagram has previously been prevented due to the seven orders of magnitude that need to be traversed across the kinetic parameter; however, the cobalt catalyst in this study provided unique control of this parameter. By varying the acids used as the proton source, the rate constants for protonation were tuned via a pKa-dependent linear free energy relationship. PMID:26998812
Danel, J-F.; Kazandjian, L.; Zerah, G.
2009-06-15
A form of the linear mixing rule involving the equality of excess pressures is tested with various mole fractions and various types of orbital-free molecular dynamics simulations. For all the cases considered, this mixing rule yields, within statistical error, the pressure of a mixture of helium and iron obtained by a direct simulation. In an attempt to interpret the robustness of the mixing rule, we show that it can be derived from thermodynamic stability if the system is regarded as a mixture of independent effective average atoms. The success of the mixing rule applied with equations of state including various degrees of approximation leads us to suggest its use in the thermodynamic domain where quantum molecular dynamics can be implemented.
Koseki, Jun; Gouda, Hiroaki; Hirono, Shuichi
2016-01-01
The natural product argadin is a cyclopentapeptide chitinase inhibitor that binds to chitinase B (ChiB) from the pathogenic bacteria Serratia marcescens. N(ω)-Acetyl-L-arginine and L-aminoadipic acid of argadin form intramolecular ionic hydrogen bonds in the aromatic hydrophobic pocket of ChiB. We performed ab initio molecular orbital and density functional theory calculations to elucidate the role of this intramolecular hydrogen bonding on intermolecular interactions between argadin and ChiB. We found that argadin accrues large stabilization energies from the van der Waals dispersion interactions, such as CH-π, π-π, and π-lone pair interactions, in the aromatic hydrophobic pocket of ChiB, although intramolecular hydrogen bonding within argadin might result in loss of entropy. The intramolecular ionic hydrogen bonding formation canceled local molecular charges and provided good van der Waals interactions with surrounding aromatic residues. PMID:27373666
Nishimoto, Yoshio; Nakata, Hiroya; Fedorov, Dmitri G; Irle, Stephan
2015-12-17
The fully analytic gradient is developed for density-functional tight-binding (DFTB) combined with the fragment molecular orbital (FMO) method (FMO-DFTB). The response terms arising from the coupling of the electronic state to the embedding potential are derived, and the gradient accuracy is demonstrated on water clusters and a polypeptide. The radial distribution functions (RDFs) obtained with FMO-DFTB are found to be similar to those from conventional DFTB, while the computational cost is greatly reduced; for 256 water molecules one molecular dynamics (MD) step takes 73.26 and 0.68 s with full DFTB and FMO-DFTB, respectively, showing a speed-up factor of 108. FMO-DFTB/MD is applied to 100 ps MD simulations of liquid hydrogen halides and is found to reproduce experimental RDFs reasonably well. PMID:26623658
Mean-Field Theory of Intra-Molecular Charge Ordering in (TTM--TTP)I3
NASA Astrophysics Data System (ADS)
Omori, Yukiko; Tsuchiizu, Masahisa; Suzumura, Yoshikazu
2011-02-01
We examine an intra-molecular charge-ordered (ICO) state in the multi-orbital molecular compound (TTM--TTP)I3 on the basis of an effective two-orbital model derived from ab initio calculations. Representing the model in terms of the fragment molecular-orbital (MO) picture, the ICO state is described as the charge disproportionation on the left and right fragment MOs. By applying the mean-field theory, the phase diagram of the ground state is obtained as a function of the inter-molecular Coulomb repulsion and the intra-molecular transfer integral. The ICO state is stabilized by large inter-fragment Coulomb interactions, and the small intra-molecular transfer energy between two fragment MOs. Furthermore, we examine the finite-temperature phase diagram. The relevance to the experimental observations in the molecular compound of (TTM--TTP)I3 is also discussed.
NASA Astrophysics Data System (ADS)
Xavier, S.; Periandy, S.; Carthigayan, K.; Sebastian, S.
2016-12-01
Vibrational spectral analysis of Diphenyl Carbonate (DPC) is carried out by using FT-IR and FT-Raman spectroscopic techniques. It is found that all vibrational modes are in the expected region. Gaussian computational calculations were performed using B3LYP method with 6-311++G (d, p) basis set. The computed geometric parameters are in good agreement with XRD data. The observation shows that the structure of the carbonate group is unsymmetrical by ∼5° due to the attachment of the two phenyl rings. The stability of the molecule arising from hyperconjugative interaction and charge delocalization are analyzed by Natural Bond Orbital (NBO) study and the results show the lone pair transition has higher stabilization energy compared to all other. The 1H and 13C NMR chemical shifts are calculated using the Gauge-Including Atomic Orbital (GIAO) method with B3LYP/6-311++G (d, p) method. The chemical shifts computed theoretically go very closer to the experimental results. A study on the electronic and optical properties; absorption wavelengths, excitation energy, dipole moment and frontier molecular orbital energies and Molecular electrostatic potential (MEP) exhibit the high reactivity nature of the molecule. The non-linear optical property of the DPC molecule predicted theoretically found to be good candidate for NLO material. TG/DTA analysis was made and decomposition of the molecule with respect to the temperature was studied. DPC having the anthelmintic activity is docked in the Hemoglobin of Fasciola hepatica protein. The DPC has been screened to antimicrobial activity and found to exhibit antibacterial effects.
Molecular quantum magnetism with strong spin-orbit coupling in inorganic solid Ba3Yb2Zn5O11
NASA Astrophysics Data System (ADS)
Park, Sang-Youn; Ji, Sungdae; Park, Jae-Hoon; Do, Seunghwan; Choi, Kwang-Yong; Jang, Dongjin; Schmidt, Burkhard; Brando, Manuel; Butch, Nicholas
The molecular magnet, assembly of finite number of spins which are isolated from environment, is a model system to study the quantum information process such as the qubit or spintronic devices. In past decades, the molecular magnet has been mostly realized in organic material, however, it has difficulty synthesizing materials or controlling their properties, meanwhile tremendous endeavors to search inorganic molecular magnet are continuing. Here, we propose Ba3Yb2Zn5O11 as a candidate of inorganic molecular magnet. This material consists of an alternating 3D-array of small and large tetrahedron containing antiferromagnetically coupled four pseudospin-1/2 Yb ions, and magnetic properties are described by an isolated tetrahedron without long-range magnetic ordering. Inelastic neutron scattering measurement with external magnetic field reveals that extraordinarily huge Dzyaloshinsky-Moriya (DM) interaction originating from strong spin-orbit coupling in Yb isospin is the key to explain energy level of tetrahedron in addition to Heisenberg exchange interaction and Zeeman effect. Magnetization measurement shows the Landau-Zener transition between avoided crossing levels caused by DM interaction.
Gravity wave transmission diagram
NASA Astrophysics Data System (ADS)
Tomikawa, Yoshihiro
2016-07-01
A possibility of gravity wave propagation from a source region to the airglow layer around the mesopause has been discussed based on the gravity wave blocking diagram taking into account the critical level filtering alone. This paper proposes a new gravity wave transmission diagram in which both the critical level filtering and turning level reflection of gravity waves are considered. It shows a significantly different distribution of gravity wave transmissivity from the blocking diagram.
Molecular Orbital Calculations on the Interaction of PF3 and CO with Ni and Cu Atoms
NASA Astrophysics Data System (ADS)
Itoh, H.; Ertl, G.
1982-04-01
The interaction of Ni and Cu atoms with PF3 and CO ligands was investigated by means of ab initio MO calculations. Coupling occurs mainly through the HOMO (8a1 for PF3 and 5 a for CO) levels with the metal 4s and 3dz2 orbitals, if the 3dn-1 4s1 electronic configuration of the metal atom is considered. The estimated M-PF3 bond lengths are 2.0 Å for Ni and 2.5 Å for Cu. Calculations with Ni(3d10) revealed for PF3 a more pronounced electron transfer to the ligand than for CO. The results are consistent with experimental UPS data for mononuclear complexes as well as corresponding adsorption systems. In particular, split-off d-states observed in UPS data for adsorbed PF3 are attributed to the pronounced lowering of the 3d-orbital energy of the metal atom upon interaction with this electropositive ligand.
Direct observation of collective modes coupled to molecular orbital-driven charge transfer
NASA Astrophysics Data System (ADS)
Ishikawa, Tadahiko; Hayes, Stuart A.; Keskin, Sercan; Corthey, Gastón; Hada, Masaki; Pichugin, Kostyantyn; Marx, Alexander; Hirscht, Julian; Shionuma, Kenta; Onda, Ken; Okimoto, Yoichi; Koshihara, Shin-ya; Yamamoto, Takashi; Cui, Hengbo; Nomura, Mitsushiro; Oshima, Yugo; Abdel-Jawad, Majed; Kato, Reizo; Miller, R. J. Dwayne
2015-12-01
Correlated electron systems can undergo ultrafast photoinduced phase transitions involving concerted transformations of electronic and lattice structure. Understanding these phenomena requires identifying the key structural modes that couple to the electronic states. We report the ultrafast photoresponse of the molecular crystal Me4P[Pt(dmit)2]2, which exhibits a photoinduced charge transfer similar to transitions between thermally accessible states, and demonstrate how femtosecond electron diffraction can be applied to directly observe the associated molecular motions. Even for such a complex system, the key large-amplitude modes can be identified by eye and involve a dimer expansion and a librational mode. The dynamics are consistent with the time-resolved optical study, revealing how the electronic, molecular, and lattice structures together facilitate ultrafast switching of the state.
Kortright, Jeffrey B; Kortright, Jeffrey B; Lincoln, Derek M; Edelstein, Ruth Shima; Epstein, Arthur J
2008-05-20
X-ray absorption spectroscopy (XAS) and magnetic circular dichroism (MCD) at the V L2,3 and C and N K edges reveal bonding/backbonding interactions in films of the 400 K magnetic semiconductor V[TCNE]x~;;2. In V spectra, dxy-like orbitals are modeled assuming V2+ in an octahedral ligand field, while dz2 and dx2-y2 orbitals involved in strong covalent bonding cannot be modeled by atomic calculations. C and N MCD, and differences in XAS from neutral TCNE molecules, reveal spin-polarized molecular orbitals in V[TCNE]x~;;2 associated with backbonding interactions that yield its novel properties.
Ohyama, Tatsuya; Hayakawa, Masato; Nishikawa, Shin; Kurita, Noriyuki
2011-06-01
Transcription mechanisms of gene information from DNA to mRNA are essentially controlled by regulatory proteins such as a lactose repressor (LacR) protein and ligand molecules. Biochemical experiments elucidated that a ligand binding to LacR drastically changes the mechanism controlled by LacR, although the effect of ligand binding has not been clarified at atomic and electronic levels. We here investigated the effect of ligand binding on the specific interactions between LacR and operator DNA by the molecular simulations combined with classical molecular mechanics and ab initio fragment molecular orbital methods. The results indicate that the binding of anti-inducer ligand strengthens the interaction between LacR and DNA, which is consistent with the fact that the binding of anti-inducer enhances the repression of gene transcription by LacR. It was also elucidated that hydrating water molecules existing between LacR and DNA contribute to the specific interactions between LacR and DNA. PMID:21328406
Stasyuk, O A; Szatylowicz, H; Krygowski, T M; Fonseca Guerra, C
2016-04-28
The substituent effect of the amino and nitro groups on the electronic system of benzene has been investigated quantum chemically using quantitative Kohn-Sham molecular orbital theory and a corresponding energy decomposition analysis (EDA). The directionality of electrophilic substitution in aniline can accurately be explained with the amount of contribution of the 2pz orbitals on the unsubstituted carbon atoms to the highest occupied π orbital. For nitrobenzene, the molecular π orbitals cannot explain the regioselectivity of electrophilic substitution as there are two almost degenerate π orbitals with nearly the same 2pz contributions on the unsubstituted carbon atoms. The Voronoi deformation density analysis has been applied to aniline and nitrobenzene to obtain an insight into the charge rearrangements due to the substituent. This analysis method identified the orbitals involved in the C-N bond formation of the π system as the cause for the π charge accumulation at the ortho and para positions in the case of the NH2 group and the largest charge depletion at these same positions for the NO2 substituent. Furthermore, we showed that it is the repulsive interaction between the πHOMO of the phenyl radical and the πHOMO of the NH2 radical that is responsible for pushing up the πHOMO of aniline and therefore activating this π orbital of the phenyl ring towards electrophilic substitution. PMID:26800159
DFT calculations of molecular excited states using an orbital-dependent nonadiabatic exchange kernel
Ipatov, A. N.
2010-02-15
A density functional method for computing molecular excitation spectra is presented that uses a frequency-dependent kernel and takes into account the nonlocality of exchange interaction. Owing to its high numerical stability and the use of a nonadiabatic (frequency-dependent) exchange kernel, the proposed approach provides a qualitatively correct description of the asymptotic behavior of charge-transfer excitation energies.
NASA Astrophysics Data System (ADS)
Pati, Ranjit; Karna, Shashi P.
2002-01-01
The dependence of electron transfer (ET) coupling element, VAB, on the length of rigid-rod-like systems consisting of bicyclo[1.1.1]pentane (BCP), cubane (CUB), and bicyclo[2.2.2]octane (BCO) monomers, has been investigated with the use of ab initio Hartree-Fock (HF) method employing Marcus-Hush two-state (TS) model. The value of VAB decreases exponentially with increase in the number of the cage units of the σ-bonded molecules. The calculated decay constant, β, shows good agreement with previously reported data. For molecular length⩾15 Å, the value of VAB becomes negligibly small, suggesting complete suppression of the through bond direct tunneling contribution to ET process.
Controlling the Interference of Multiple Molecular Orbitals in High-Harmonic Generation
Woerner, H. J.; Bertrand, J. B.; Hockett, P.; Corkum, P. B.; Villeneuve, D. M.
2010-06-11
We demonstrate a new method to investigate the origin of spectral structures in high-harmonic generation. We report detailed measurements of high-harmonic spectra in aligned nitrogen and carbon dioxide molecules. Varying the wavelength and intensity of the generating laser field, we show that the minimum in aligned N{sub 2} molecules is nearly unaffected, whereas the minimum in aligned CO{sub 2} molecules shifts over more than 15 eV. Our quantitative analysis shows that both the interference of multiple orbitals and their structural characteristics affect the position of the minimum. Our method provides a simple approach to the investigation of the high-harmonic generation process in more complex molecules.
Rings of C2H in the Molecular Disks Orbiting TW Hya and V4046 Sgr
NASA Astrophysics Data System (ADS)
Kastner, J. H.; Qi, C.; Gorti, U.; Hily-Blant, P.; Oberg, K.; Forveille, T.; Andrews, S.; Wilner, D.
2016-01-01
We have used the Submillimeter Array (SMA) to image, at ~1'' resolution, C2H(3-2) emission from the molecule-rich circumstellar disks orbiting the nearby, classical T Tauri star systems TW Hya and V4046 Sgr. The SMA imaging reveals that the C2H emission exhibits a ring-like morphology within each disk; the radius of the inner hole of the C2H ring within the V4046 Sgr disk (~70 AU) is somewhat larger than than of its counterpart within the TW Hya disk (~45 AU). We suggest that, in each case, the C2H emission likely traces irradiation of the tenuous surface layers of the outer disks by high-energy photons from the central stars.
Nishimoto, Yoshio; Fedorov, Dmitri G
2016-08-10
The energy and its analytic gradient are formulated for the fragment molecular orbital (FMO) method combined with density-functional tight-binding (DFTB) and the polarizable continuum model (PCM). The accuracy is demonstrated in comparison with unfragmented calculations and numerical gradients. The instability in the description of proteins using density functional theory (DFT) and DFTB is analyzed for both unfragmented and FMO methods. The cause of the instability is shown to be charged residues, and the problem is particularly severe in the gas phase when long-range functionals are not used. Adding solvent effects considerably increases the gap between occupied and virtual orbitals and stabilizes convergence. The pair interaction energies calculated using FMO-DFT and FMO-DFTB in solution are shown to correlate, whereas the latter method is 4840 times faster than the former for a protein consisting of 1961 atoms. The structures of five proteins (containing up to 3578 atoms) optimized using FMO-DFTB/PCM agree reasonably well with experiment. PMID:27215663
Yang, Hongping; Yam, Chi-Yung; Zhang, Aihua; Xu, Zhiping; Luo, Jun; Zhu, Jing
2015-03-21
The highest occupied molecular orbital (HOMO) energies of carbon nanotubes (CNTs) and graphene are crucial in fundamental and applied research of carbon nanomaterials, and so their modulation is desired. Our first-principles calculations reveal that the HOMO energies of CNTs and graphene can both be raised by negatively charging, and that the rate of increase of the HOMO energy of a CNT is much greater and faster than that of graphene with the same number of C atoms. This discriminative modulation holds true regardless of the number of C atoms and the CNT type, and so is universal. This work provides a new opportunity to develop all-carbon devices with CNTs and graphene as different functional elements. PMID:25692228
NASA Astrophysics Data System (ADS)
Ilieva, S.; Hadjieva, B.; Galabov, B.
1999-09-01
Ab initio molecular orbital calculations at HF/4-31G level and infrared spectroscopic data for the frequencies are applied to analyse the grouping in a series model aromatic secondary amides: formanilide; acetanilide; o-methylacetanilide; 2,6-dimethylformanilide, 2,6-dimethylacetanilide; N-benzylacetamide and N-benzylformamide. The theoretical and experimental data obtained show that the conformational state of the molecules studied is determined by the fine balance of several intramolecular factors: resonance effect between the amide group and the aromatic ring, steric interaction between various substituents around the -NH-CO- grouping in the aromatic ring, conjugation between the carbonyl bond and the nitrogen lone pair as well as direct field influences inside the amide group.
NASA Astrophysics Data System (ADS)
Tagami, Ayumu; Ishibashi, Nobuhiro; Kato, Dai-ichiro; Taguchi, Naoki; Mochizuki, Yuji; Watanabe, Hirofumi; Ito, Mika; Tanaka, Shigenori
2009-04-01
Bioluminescence spectra of firefly Luciola cruciata were theoretically analyzed on the basis of the fragment molecular orbital (FMO) method. The CIS(D) and PR-CIS(Ds) methods were employed for the calculations of emission energies of wild-type and mutant luciferase-oxyluciferin systems, and various multi-layer FMO calculations were performed changing the sizes of the luciferase protein and of the chromophore to which the excited-state calculations were applied. We have thus reproduced the experimental emission energies of wild-type and mutant luciferase systems with good accuracy, which provides useful information concerning the roles of protein environment for the color tuning of the bioluminescence spectra of firefly.
NASA Technical Reports Server (NTRS)
Defrees, D. J.; Mclean, A. D.
1986-01-01
The discovery of cyclopropenylidene in space suggests that other C3H2 isomers may be present, and a tentative detection of one such isomer, propargylene (HCCCH), has been reported. Ab initio molecular orbital theory has been used to characterize five low-lying, metastable isomers of cyclopropenylidene. Extended calculations including the electron correlation energy, show that the lowest in energy is singlet propadienylidene, followed by propargylene; the singlet and triplet of the latter are too close in energy to allow an assignment of the ground state; triplet propadienylidene is at a significantly higher energy. Rotational frequencies computed to an expected accuracy of + or - 1 -2 percent do not confirm the tentative detection of propargylene in space, although the discrepancy between theory and the observation is not so great as to unequivocally rule out this possibility.
NASA Astrophysics Data System (ADS)
Chiosi, C.; Murdin, P.
2000-11-01
The Hertzsprung-Russell diagram (HR-diagram), pioneered independently by EJNAR HERTZSPRUNG and HENRY NORRIS RUSSELL, is a plot of the star luminosity versus the surface temperature. It stems from the basic relation for an object emitting thermal radiation as a black body: ...
Jafri, Sadaf Fatima; Koumousi, Evangelia S; Sainctavit, Philippe; Juhin, Amélie; Schuler, Vivien; Bun U, Oana; Mitcov, Dmitri; Dechambenoit, Pierre; Mathonière, Corine; Clérac, Rodolphe; Otero, Edwige; Ohresser, Philippe; Cartier Dit Moulin, Christophe; Arrio, Marie-Anne
2016-07-18
Photomagnetism in three-dimensional Co/Fe Prussian blue analogues is a complex phenomenon, whose detailed mechanism is not yet fully understood. Recently, researchers have been able to prepare molecular fragments of these networks using a building block synthetic approach from mononuclear precursors. The main objective in this strategy is to isolate the smallest units that show an intramolecular electron transfer to have a better understanding of the electronic processes. A prior requirement to the development of this kind of system is to understand to what extent electronic and magnetic properties are inherited from the corresponding precursors. In this work, we investigate the electronic and magnetic properties of the FeTp precursor (N(C4H9)4)[TpFe(III)(CN)3], (Tp being tris-pyrazolyl borate) of a recently reported binuclear cyanido-bridged Fe/Co complex. X-ray absorption spectroscopy and X-ray magnetic circular dichroism measurements at the Fe L2,3 edges (2p → 3d) supported by ligand field multiplet calculations have allowed to determine the spin and orbit magnetic moments. Inaccuracy of the spin sum rule in the case of low-spin Fe(III) ion was demonstrated. An exceptionally large value of the orbital magnetic moment is found (0.9 μB at T = 2 K and B = 6.5 T) that is likely to play an important role in the magnetic and photomagnetic properties of molecular Fe/Co Prussian blue analogues. PMID:27385292
Thirman, Jonathan Head-Gordon, Martin
2015-08-28
An energy decomposition analysis (EDA) of intermolecular interactions is proposed for second-order Møller–Plesset perturbation theory (MP2) based on absolutely localized molecular orbitals (ALMOs), as an extension to a previous ALMO-based EDA for self-consistent field methods. It decomposes the canonical MP2 binding energy by dividing the double excitations that contribute to the MP2 wave function into classes based on how the excitations involve different molecules. The MP2 contribution to the binding energy is decomposed into four components: frozen interaction, polarization, charge transfer, and dispersion. Charge transfer is defined by excitations that change the number of electrons on a molecule, dispersion by intermolecular excitations that do not transfer charge, and polarization and frozen interactions by intra-molecular excitations. The final two are separated by evaluations of the frozen, isolated wave functions in the presence of the other molecules, with adjustments for orbital response. Unlike previous EDAs for electron correlation methods, this one includes components for the electrostatics, which is vital as adjustment to the electrostatic behavior of the system is in some cases the dominant effect of the treatment of electron correlation. The proposed EDA is then applied to a variety of different systems to demonstrate that all proposed components behave correctly. This includes systems with one molecule and an external electric perturbation to test the separation between polarization and frozen interactions and various bimolecular systems in the equilibrium range and beyond to test the rest of the EDA. We find that it performs well on these tests. We then apply the EDA to a halogen bonded system to investigate the nature of the halogen bond.
Yoshida, Norio
2014-06-07
The three-dimensional reference interaction site model (3D-RISM) method was efficiently implemented in the fragment molecular orbital (FMO) method. The method is referred to as the FMO/3D-RISM method, and allows us to treat electronic structure of the whole of a macromolecule, such as a protein, as well as the solvent distribution around a solute macromolecule. The formalism of the FMO/3D-RISM method, for the computationally available form and variational expressions, are proposed in detail. A major concern leading to the implementation of the method was decreasing the computational costs involved in calculating the electrostatic potential, because the electrostatic potential is calculated on numerous grid points in three-dimensional real space in the 3D-RISM method. In this article, we propose a procedure for decreasing the computational costs involved in calculating the electrostatic potential in the FMO method framework. The strategy involved in this procedure is to evaluate the electrostatic potential and the solvated Fock matrix in different manners, depending on the distance between the solute and the solvent. The electrostatic potential is evaluated directly in the vicinity of the solute molecule by integrating the molecular orbitals of monomer fragments of the solute molecule, whereas the electrostatic potential is described as the sum of multipole interactions when an analog of the fast multipole method is used. The efficiency of our method was demonstrated by applying it to a water trimer system and three biomolecular systems. The FMO/3D-RISM calculation can be performed within a reasonable computational time, retaining the accuracy of some physical properties.
Localized and Spectroscopic Orbitals: Squirrel Ears on Water.
ERIC Educational Resources Information Center
Martin, R. Bruce
1988-01-01
Reexamines the electronic structure of water considering divergent views. Discusses several aspects of molecular orbital theory using spectroscopic molecular orbitals and localized molecular orbitals. Gives examples for determining lowest energy spectroscopic orbitals. (ML)
NASA Astrophysics Data System (ADS)
Kim, Minjae; Choi, Hong Chul; Shim, Ji Hoon; Min, B. I.
2014-03-01
We have studied correlated electronic structures and the phase diagram of electron-doped hydrocarbon molecular solids, based on the dynamical mean-field theory. We have determined the phase diagram of hydrocarbon molecular solids as functions of doping and energy parameters including the Coulomb correlation, the Hund coupling, and the molecular-orbital (MO) energy level splitting. We have found that the hydrocarbon superconductors (electron-doped picene and coronene) belong to the multi-band Fermi liquid state, while non-superconducting electron-doped pentacene belongs to the single-band state in the proximity of the metal-insulator transition. The size of the MO energy level splitting plays an important role in deriving the superconductivity of electron-doped hydrocarbon solids. The multi-band nature of hydrocarbon solids from the small MO energy level splitting boosts the superconductivity through the enhanced density of states at the Fermi level.
Infrared spectroscopy of molecular ions in selected rotational and spin-orbit states
NASA Astrophysics Data System (ADS)
Jacovella, U.; Agner, J. A.; Schmutz, H.; Deiglmayr, J.; Merkt, F.
2016-07-01
First results are presented obtained with an experimental setup developed to record IR spectra of rotationally state-selected ions. The method we use is a state-selective version of a method developed by Schlemmer et al. [Int. J. Mass Spectrom. 185, 589 (1999); J. Chem. Phys. 117, 2068 (2002)] to record IR spectra of ions. Ions are produced in specific rotational levels using mass-analyzed-threshold-ionization spectroscopy. The state-selected ions generated by pulsed-field ionization of Rydberg states of high principal quantum number (n ≈ 200) are extracted toward an octupole ion guide containing a neutral target gas. Prior to entering the octupole, the ions are excited by an IR laser. The target gas is chosen so that only excited ions react to form product ions. These product ions are detected mass selectively as a function of the IR laser wavenumber. To illustrate this method, we present IR spectra of C 2 H2 + in selected rotational levels of the 2Πu,3/2 and 2Πu,1/2 spin-orbit components of the vibronic ground state.
Infrared spectroscopy of molecular ions in selected rotational and spin-orbit states.
Jacovella, U; Agner, J A; Schmutz, H; Deiglmayr, J; Merkt, F
2016-07-01
First results are presented obtained with an experimental setup developed to record IR spectra of rotationally state-selected ions. The method we use is a state-selective version of a method developed by Schlemmer et al. [Int. J. Mass Spectrom. 185, 589 (1999); J. Chem. Phys. 117, 2068 (2002)] to record IR spectra of ions. Ions are produced in specific rotational levels using mass-analyzed-threshold-ionization spectroscopy. The state-selected ions generated by pulsed-field ionization of Rydberg states of high principal quantum number (n ≈ 200) are extracted toward an octupole ion guide containing a neutral target gas. Prior to entering the octupole, the ions are excited by an IR laser. The target gas is chosen so that only excited ions react to form product ions. These product ions are detected mass selectively as a function of the IR laser wavenumber. To illustrate this method, we present IR spectra of C2H2 (+) in selected rotational levels of the (2)Πu,3/2 and (2)Πu,1/2 spin-orbit components of the vibronic ground state. PMID:27394102
Conradie, Jeanet; Ghosh, Abhik
2014-05-19
The great majority of {NiNO}(10) complexes are characterized by short Ni-N(O) distances of 1.60-1.65 Å and linear NO units. Against this backdrop, the {CuNO}(10) unit in the recently reported [Cu(CH3NO2)5(NO)](2+) cation (1) has a CuNO angle of about 120° and a very long 1.96 Å Cu-N(O) bond. According to DFT calculations, metal-NO bonding in 1 consists of a single Cu(dz(2))-NO(π*) σ-interaction and essentially no metal(dπ)-NO(π*) π-bonding, which explains both the bent CuNO geometry and the long, weak Cu-N(O) bond. This σ-interaction is strongly favored by a ligand trans to the NO; indeed such a trans ligand may be critical for the existence and stability of a {CuNO}(10) unit. By contrast, {NiNO}(10) complexes exhibit a strong avoidance of such trans ligands. Thus, a five-coordinate {NiNO}(10) complex appears to favor a trigonal-bipyramidal structure with the NO in an equatorial position, as in the case of [Ni(bipy)2(NO)](+) (6). An unusual set of Ni(d)-NO(π*) orbital interactions accounts for the strongly bent NiNO geometry for this complex. PMID:24796643
The continuous and discrete molecular orbital x-ray bands from Xeq+ (12≤q≤29) +Zn collisions
NASA Astrophysics Data System (ADS)
Guo, Yipan; Yang, Zhihu; Hu, Bitao; Wang, Xiangli; Song, Zhangyong; Xu, Qiumei; Zhang, Boli; Chen, Jing; Yang, Bian; Yang, Jie
2016-07-01
In this paper, the x-ray emissions are measured by the interaction of 1500–3500 keV Xeq+ (q = 12, 15, 17, 19, 21, 23, 26 and 29) ions with Zn target. When q < 29, we observe Ll, Lα, Lβ1, Lβ2 and Lγ characteristic x-rays from Xeq+ ions and a broad M-shell molecular orbital (MO) x-ray band from the transient quasi-molecular levels. It is found that their yields quickly increase with different rates as the incident energy increases. Besides, the widths of the broad MO x-ray bands are about 0.9–1.32 keV over the energy range studied and are proportional to v1/2 (v = projectile velocity). Most remarkably, when the projectile charge state is 29, the broad x-ray band separates into several narrow discrete spectra, which was never observed before in this field.
Pals, Justin A; Wagner, Elizabeth D; Plewa, Michael J
2016-03-15
Disinfection of drinking water protects public health against waterborne pathogens. However, during disinfection, toxic disinfection byproducts (DBPs) are formed. Exposure to DBPs was associated with increased risk of bladder cancer in humans. DBPs are generated at concentrations below their carcinogenic potencies; it is unclear how exposure leads to adverse health outcomes. We used computational estimates of the energy of the lowest unoccupied molecular orbital (ELUMO) to predict thiol reactivity and additive toxicity among soft electrophile DBPs. Bromoacetic acid (BAA) was identified as non-thiol-reactive, which was supported by in chemico and in vitro data. Bromoacetonitrile (BAN) and bromoacetamide (BAM) were thiol-reactive. Genotoxicity induced by these compounds was reduced by increasing the thiol pool with N-acetyl l-cysteine (NAC), while NAC had little effect on BAA. BAN and BAM shared depletion of glutathione (GSH) or cellular thiols as a molecular initiating event (MIE), whereas BAA induces toxicity through another pathway. Binary mixtures of BAM and BAN expressed a potentiating effect in genotoxicity. We found that soft electrophile DBPs could be an important predictor of common mechanism groups that demonstrated additive toxicity. In silico estimates of ELUMO could be used to identify the most relevant DBPs that are the forcing factors of the toxicity of finished drinking waters. PMID:26854864
Lyu, Yan; Fang, Yuan; Miao, Qingqing; Zhen, Xu; Ding, Dan; Pu, Kanyi
2016-04-26
Optical theranostic nanoagents that seamlessly and synergistically integrate light-generated signals with photothermal or photodynamic therapy can provide opportunities for cost-effective precision medicine, while the potential for clinical translation requires them to have good biocompatibility and high imaging/therapy performance. We herein report an intraparticle molecular orbital engineering approach to simultaneously enhance photoacoustic brightness and photothermal therapy efficacy of semiconducting polymer nanoparticles (SPNs) for in vivo imaging and treatment of cancer. The theranostic SPNs have a binary optical component nanostructure, wherein a near-infrared absorbing semiconducting polymer and an ultrasmall carbon dot (fullerene) interact with each other to induce photoinduced electron transfer upon light irradiation. Such an intraparticle optoelectronic interaction augments heat generation and consequently enhances the photoacoustic signal and maximum photothermal temperature of SPNs by 2.6- and 1.3-fold, respectively. With the use of the amplified SPN as the theranostic nanoagent, it permits enhanced photoacoustic imaging and photothermal ablation of tumor in living mice. Our study thus not only introduces a category of purely organic optical theranostics but also highlights a molecular guideline to amplify the effectiveness of light-intensive imaging and therapeutic nanosystems. PMID:26959505
Pals, Justin A.; Wagner, Elizabeth D.; Plewa, Michael J.
2016-01-01
Disinfection of drinking water protects public health against waterborne pathogens. However, during disinfection, toxic disinfection byproducts (DBPs) are formed. Exposure to DBPs was associated with increased risk of bladder cancer in humans. DBPs are generated at concentrations below their carcinogenic potencies; it is unclear how exposure leads to adverse health outcomes. We used computational estimates of the energy of the lowest unoccupied molecular orbital (ELUMO) to predict thiol reactivity and additive toxicity among soft electrophile DBPs. Bromoacetic acid (BAA) was identified as non-thiol-reactive, which was supported by in chemico and in vitro data. Bromoacetonitrile (BAN) and bromoacetamide (BAM) were thiol-reactive. Genotoxicity induced by these compounds was reduced by increasing the thiol pool with N-acetyl l-cysteine (NAC), while NAC had little effect on BAA. BAN and BAM shared depletion of glutathione (GSH) or cellular thiols as a molecular initiating event (MIE), whereas BAA induces toxicity through another pathway. Binary mixtures of BAM and BAN expressed a potentiating effect in genotoxicity. We found that soft electrophile DBPs could be an important predictor of common mechanism groups that demonstrated additive toxicity. In silico estimates of ELUMO could be used to identify the most relevant DBPs that are the forcing factors of the toxicity of finished drinking waters. PMID:26854864
The continuous and discrete molecular orbital x-ray bands from Xeq+ (12≤q≤29) +Zn collisions
Guo, Yipan; Yang, Zhihu; Hu, Bitao; Wang, Xiangli; Song, Zhangyong; Xu, Qiumei; Zhang, Boli; Chen, Jing; Yang, Bian; Yang, Jie
2016-01-01
In this paper, the x-ray emissions are measured by the interaction of 1500–3500 keV Xeq+ (q = 12, 15, 17, 19, 21, 23, 26 and 29) ions with Zn target. When q < 29, we observe Ll, Lα, Lβ1, Lβ2 and Lγ characteristic x-rays from Xeq+ ions and a broad M-shell molecular orbital (MO) x-ray band from the transient quasi-molecular levels. It is found that their yields quickly increase with different rates as the incident energy increases. Besides, the widths of the broad MO x-ray bands are about 0.9–1.32 keV over the energy range studied and are proportional to v1/2 (v = projectile velocity). Most remarkably, when the projectile charge state is 29, the broad x-ray band separates into several narrow discrete spectra, which was never observed before in this field. PMID:27469425
Heifetz, Alexander; Trani, Giancarlo; Aldeghi, Matteo; MacKinnon, Colin H; McEwan, Paul A; Brookfield, Frederick A; Chudyk, Ewa I; Bodkin, Mike; Pei, Zhonghua; Burch, Jason D; Ortwine, Daniel F
2016-05-12
Inhibition of inducible T-cell kinase (ITK), a nonreceptor tyrosine kinase, may represent a novel treatment for allergic asthma. In our previous reports, we described the discovery of sulfonylpyridine (SAP), benzothiazole (BZT), indazole (IND), and tetrahydroindazole (THI) series as novel ITK inhibitors and how computational tools such as dihedral scans and docking were used to support this process. X-ray crystallography and modeling were applied to provide essential insight into ITK-ligand interactions. However, "visual inspection" traditionally used for the rationalization of protein-ligand affinity cannot always explain the full complexity of the molecular interactions. The fragment molecular orbital (FMO) quantum-mechanical (QM) method provides a complete list of the interactions formed between the ligand and protein that are often omitted from traditional structure-based descriptions. FMO methodology was successfully used as part of a rational structure-based drug design effort to improve the ITK potency of high-throughput screening hits, ultimately delivering ligands with potency in the subnanomolar range. PMID:26950250
Horn, Paul R; Mao, Yuezhi; Head-Gordon, Martin
2016-08-17
An energy decomposition analysis (EDA) separates a calculated interaction energy into as many interpretable contributions as possible; for instance, permanent and induced electrostatics, Pauli repulsions, dispersion and charge transfer. The challenge is to construct satisfactory definitions of all terms in the chemically relevant regime where fragment densities overlap, rendering unique definitions impossible. Towards this goal, we present an improved EDA for Kohn-Sham density functional theory (DFT) with properties that have previously not been simultaneously attained. Building on the absolutely localized molecular orbital (ALMO)-EDA, this second generation ALMO-EDA is variational and employs valid antisymmetric electronic wavefunctions to produce all five contributions listed above. These contributions moreover all have non-trivial complete basis set limits. We apply the EDA to the water dimer, the T-shaped and parallel-displaced benzene dimer, the p-biphthalate dimer "anti-electrostatic" hydrogen bonding complex, the biologically relevant binding of adenine and thymine in stacked and hydrogen-bonded configurations, the triply hydrogen-bonded guanine-cytosine complex, the interaction of Cl(-) with s-triazine and with the 1,3-dimethyl imidazolium cation, which is relevant to the study of ionic liquids, and the water-formaldehyde-vinyl alcohol ter-molecular radical cationic complex formed in the dissociative photoionization of glycerol. PMID:27492057
Chang, Le; Ishikawa, Takeshi; Kuwata, Kazuo; Takada, Shoji
2013-05-30
Accurate computational estimate of the protein-ligand binding affinity is of central importance in rational drug design. To improve accuracy of the molecular mechanics (MM) force field (FF) for protein-ligand simulations, we use a protein-specific FF derived by the fragment molecular orbital (FMO) method and by the restrained electrostatic potential (RESP) method. Applying this FMO-RESP method to two proteins, dodecin, and lysozyme, we found that protein-specific partial charges tend to differ more significantly from the standard AMBER charges for isolated charged atoms. We did not see the dependence of partial charges on the secondary structure. Computing the binding affinities of dodecin with five ligands by MM PBSA protocol with the FMO-RESP charge set as well as with the standard AMBER charges, we found that the former gives better correlation with experimental affinities than the latter. While, for lysozyme with five ligands, both charge sets gave similar and relatively accurate estimates of binding affinities. PMID:23420697
Zarycz, M. Natalia C. Provasi, Patricio F.; Sauer, Stephan P. A.
2014-10-21
We discuss the effect of electron correlation on the unexpected differential sensitivity (UDS) in the {sup 1}J(C–H) coupling constant of CH{sub 4} using a decomposition into contributions from localized molecular orbitals and compare with the {sup 1}J(N–H) coupling constant in NH{sub 3}. In particular, we discuss the well known fact that uncorrelated coupled Hartree-Fock (CHF) calculations are not able to reproduce the UDS in methane. For this purpose we have implemented for the first time a localized molecular orbital analysis for the second order polarization propagator approximation with coupled cluster singles and doubles amplitudes—SOPPA(CCSD) in the DALTON program. Comparing the changes in the localized orbital contributions at the correlated SOPPA and SOPPA(CCSD) levels and at the uncorrelated CHF level, we find that the latter overestimates the effect of stretching the bond between the coupled atoms on the contribution to the coupling from the localized bonding orbital between these atoms. This disturbs the subtle balance between the molecular orbital contributions, which lead to the UDS in methane.
NASA Astrophysics Data System (ADS)
Suresh, S.; Gunasekaran, S.; Srinivasan, S.
The solid phase FT-IR and FT-Raman spectra of 2-[2-[2-[(2,6-dichlorophenyl)amino]phenyl]acetyl] oxyacetic acid (Aceclofenac) have been recorded in the region 4000-400 and 4000-100 cm-1 respectively. The optimized molecular geometry and fundamental vibrational frequencies are interpreted with the aid of structure optimizations and normal coordinate force field calculations based on density functional theory (DFT) method and a comparative study between Hartree Fork (HF) method 6-311++G(d,p) level basis set. The calculated harmonic vibrational frequencies were scaled and have been compared with experimental by obtained FT-IR and FT-Raman spectra. A detailed interpretation of the vibrational spectra of this compound has been made on the basis of the calculated potential energy distribution (PED). The time dependent DFT method employed to study its absorption energy and oscillator strength. The linear polarizability (α) and the first order hyper polarizability (β) values of the investigated molecule have been computed. The electronic properties, such as HOMO and LUMO energies, molecular electrostatic potential (MESP) were also performed. Stability of the molecule arising from hyper conjugative interaction, charge delocalization has been analyzed using natural bond orbital (NBO) analysis.
A Ring of C2H in the Molecular Disk Orbiting TW Hya
NASA Astrophysics Data System (ADS)
Kastner, Joel H.; Qi, Chunhua; Gorti, Uma; Hily-Blant, Pierre; Oberg, Karin; Forveille, Thierry; Andrews, Sean; Wilner, David
2015-06-01
We have used the Submillimeter Array to image, at ˜1.″5 resolution, C2H N=3\\to 2 emission from the circumstellar disk orbiting the nearby (D = 54 pc), ˜8 Myr-old, ˜0.8 {{M}⊙ } classical T Tauri star TW Hya. The SMA imaging reveals that the C2H emission exhibits a ring-like morphology. Based on a model in which the C2H column density follows a truncated radial power-law distribution, we find that the inner edge of the ring lies at ˜45 AU, and that the ring extends to at least ˜120 AU. Comparison with previous (single-dish) observations of C2H N=4\\to 3 emission indicates that the C2H molecules are subthermally excited and, hence, that the emission arises from the relatively warm (T≳ 40 K), tenuous (n\\ll {{10}7} cm-3) upper atmosphere of the disk. Based on these results and comparisons of the SMA C2H map with previous submillimeter and scattered-light imaging, we propose that the C2H emission most likely traces particularly efficient photo-destruction of small grains and/or photodesorption and photodissociation of hydrocarbons derived from grain ice mantles in the surface layers of the outer disk. The presence of a C2H ring in the TW Hya disk hence likely serves as a marker of dust grain processing and radial and vertical grain size segregation within the disk.
Moncho, Salvador; Autschbach, Jochen
2010-12-01
The NMR nuclear shielding tensors for the series LaX(3), with X = F, Cl, Br and I, have been computed using two-component relativistic density functional theory based on the zeroth-order regular approximation (ZORA). A detailed analysis of the inverse halogen dependence (IHD) of the La shielding was performed via decomposition of the shielding tensor elements into contributions from localized and delocalized molecular orbitals. Both spin-orbit and paramagnetic shielding terms are important, with the paramagnetic terms being dominant. Major contributions to the IHD can be attributed to the La-X bonding orbitals, as well as to trends associated with the La core and halogen lone pair orbitals, the latter being related to X-La π donation. An 'orbital rotation' model for the in-plane π acceptor f orbital of La helps to rationalize the significant magnitude of deshielding associated with the in-plane π donation. The IHD goes along with a large increase in the shielding tensor anisotropy as X becomes heavier, which can be associated with trends for the covalency of the La-X bonds, with a particularly effective transfer of spin-orbit coupling induced spin density from iodine to La in LaI(3). PMID:20586110
ERIC Educational Resources Information Center
Halpern, Arthur M.; Glendening, Eric D.
2013-01-01
A three-part project for students in physical chemistry, computational chemistry, or independent study is described in which they explore applications of valence bond (VB) and molecular orbital-configuration interaction (MO-CI) treatments of H[subscript 2]. Using a scientific spreadsheet, students construct potential-energy (PE) curves for several…