Quantum-mechanical calculations of vibrationally resolved cross sections for non-dissociative charge transfer of O3+ with H-2

Charge transfer due to collisions of ground state O3+ (2s(2)2p P-2) ions with molecular hydrogen is investigated using the quantum-mechanical molecular-orbital close-coupling (MOCC) method, and electronic and vibrational state-selective cross sections along with the corresponding differential cross sections are calculated for projectile energies of 100, 500, 1000 and 5000 eV/u at the orientation angles of 25 degrees,45 degrees and 89 degrees. The adiabatic potentials and radial coupling matrix elements utilized in the QMOCC calculations were obtained with the spin-coupled valence-bond approach. The infinite order sudden approximation (IOSA) and the vibrational sudden approximation (VSA) are utilized to deal with the rotation of H-2 and the coupling between the electron and the vibration of H-2. It is found that the distribution of vibrationally resolved cross sections with the vibrational quantum number upsilon' of H-2(+) (upsilon') varies with the increment of the projectile energy; and the electronic and vibrational stateselective differential cross sections show similar behaviors: there is a highest platform within a very small scattering angle, beyond which the differential cross sections decrease as the scattering angle increases and lots of oscillating structures appear, where the scattering angle of the first structure decreases as E-P(-1/2) with the increment of the projectile energy E-P; and the structure and amplitude of the differential cross sections are sensitive to the orientation of molecule H-2, which provides a possibility to identify the orientations of molecule H-2 by the vibrational state-selective differential scattering processes.

New insights into NO-carbon and N2O-carbon reactions from quantum mechanical calculations

Density functional theory calculations were used to investigate the mechanisms of NO-carbon and N2O-carbon reactions. It was the first time that the importance of surface nitrogen groups was addressed in the kinetic behaviors of the NO-carbon reaction. It was found that the off-plane nitrogen groups that are adjacent to the zigzag edge sites and in-plane nitrogen groups that are located on the armchair sites make the bond energy of oxygen desorption even ca. 20% lower than that of the off-plane epoxy group adjacent to zigzag edge sites and in-plane o-quinone oxygen atoms on armchair sites; this may explain the reason why the experimentally obtained activation energy of the NO-carbon reaction is ca. 20% lower than that of the O-2-carbon reaction over 923 K. A higher ratio of oxygen atoms can be formed in the N2O-carbon reaction, because of the lower dissociation energy of N2O, which results in a higher ratio of off-plane epoxy oxygen atoms. The desorption energy of semiquinone with double adjacent off-plane oxygen groups is ca. 20% less than that of semiquinone with only one adjacent off-plane oxygen group. This may be the reason why the activation energy of N2O is also ca. 20% less than that of the O-2-carbon reaction. The new mechanism can also provide a good qualitative comparison for the relative reaction rates of NO-, N2O-, and O-2-carbon reactions. The anisotropic characters of these gas-carbon reactions can also be well explained.

Quantum mechanical simulation of nanosized metal-oxide-semiconductor field-effect transistor using empirical pseudopotentials: A comparison for charge density occupation methods

The atomistic pseudopotential quantum mechanical calculations are used to study the transport in million atom nanosized metal-oxide-semiconductor field-effect transistors. In the charge self-consistent calculation, the quantum mechanical eigenstates of closed systems instead of scattering states of open systems are calculated. The question of how to use these eigenstates to simulate a nonequilibrium system, and how to calculate the electric currents, is addressed. Two methods to occupy the electron eigenstates to yield the charge density in a nonequilibrium condition are tested and compared. One is a partition method and another is a quasi-Fermi level method. Two methods are also used to evaluate the current: one uses the ballistic and tunneling current approximation, another uses the drift-diffusion method. (C) 2009 American Institute of Physics. [doi:10.1063/1.3248262]

Quantum mechanical effects in nanometer field effect transistors

The atomistic pseudopotential quantum mechanical calculations for million atom nanosized metal-oxide-semiconductor field-effect transistors (MOSFETs) are presented. When compared with semiclassical Thomas-Fermi simulation results, there are significant differences in I-V curve, electron threshold voltage, and gate capacitance. In many aspects, the quantum mechanical effects exacerbate the problems encountered during device minimization, and it also presents different mechanisms in controlling the behaviors of a nanometer device than the classical one. (c) 2007 American Institute of Physics.

Accurate quantum mechanical calculation for the N+OH reaction

Accurate three-dimensional time-dependent quantum wave packet calculations for the N+OH reaction on the (3)A' potential energy surface [Guadagnini, Schatz, and Walch, J. Chem. Phys. 102, 774 (1995)] have been carried out. The calculations show for the first time that the initial state-selected reaction probabilities are dominated by resonance structures, and the lifetime of the resonance is generally in the subpicosecond time scale. The calculated reaction cross sections indicate that they are a decreasing function of the translational energy, which is in agreement qualitatively with the quasiclassical trajectory calculations. The rate constants obtained from the quantum mechanical calculations are consistent with the quasiclassical trajectory results and the experimental measurements. (C) 2003 American Institute of Physics.

Ability of the PM3 Quantum Mechanical Method to Model Intermolecular Hydrogen Bonding between Neutral Molecules

The PM3 semiempirical quantum-mechanical method was found to systematically describe intermolecular hydrogen bonding in small polar molecules. PM3 shows charge transfer from the donor to acceptor molecules on the order of 0.02-0.06 units of charge when strong hydrogen bonds are formed. The PM3 method is predictive; calculated hydrogen bond energies with an absolute magnitude greater than 2 kcal mol-' suggest that the global minimum is a hydrogen bonded complex; absolute energies less than 2 kcal mol-' imply that other van der Waals complexes are more stable. The geometries of the PM3 hydrogen bonded complexes agree with high-resolution spectroscopic observations, gas electron diffraction data, and high-level ab initio calculations. The main limitations in the PM3 method are the underestimation of hydrogen bond lengths by 0.1-0.2 for some systems and the underestimation of reliable experimental hydrogen bond energies by approximately 1-2 kcal mol-l. The PM3 method predicts that ammonia is a good hydrogen bond acceptor and a poor hydrogen donor when interacting with neutral molecules. Electronegativity differences between F, N, and 0 predict that donor strength follows the order F > 0 > N and acceptor strength follows the order N > 0 > F. In the calculations presented in this article, the PM3 method mirrors these electronegativity differences, predicting the F-H- - -N bond to be the strongest and the N-H- - -F bond the weakest. It appears that the PM3 Hamiltonian is able to model hydrogen bonding because of the reduction of two-center repulsive forces brought about by the parameterization of the Gaussian core-core interactions. The ability of the PM3 method to model intermolecular hydrogen bonding means reasonably accurate quantum-mechanical calculations can be applied to small biologic systems.

Crystal Coulomb energies: I. Organic donor-acceptor complexes--a review. II. Classical Ewald calculations of the Coulomb binding energy of four donor-acceptor crystals and Wurster's blue perchlorate. III. A rapid-convergence quantum-mechanical formalism for crystal electronic energies

Chapter I

Theories for organic donor-acceptor (DA) complexes in solution and in the solid state are reviewed, and compared with the available experimental data. As shown by McConnell et al. (Proc. Natl. Acad. Sci. U.S., 53, 46-50 (1965)), the DA crystals fall into two classes, the holoionic class with a fully or almost fully ionic ground state, and the nonionic class with little or no ionic character. If the total lattice binding energy 2ε₁ (per DA pair) gained in ionizing a DA lattice exceeds the cost 2ε_o of ionizing each DA pair, ε₁ + ε_o less than 0, then the lattice is holoionic. The charge-transfer (CT) band in crystals and in solution can be explained, following Mulliken, by a second-order mixing of states, or by any theory that makes the CT transition strongly allowed, and yet due to a small change in the ground state of the non-interacting components D and A (or D⁺ and A^-). The magnetic properties of the DA crystals are discussed.

Chapter II

A computer program, EWALD, was written to calculate by the Ewald fast-convergence method the crystal Coulomb binding energy E_C due to classical monopole-monopole interactions for crystals of any symmetry. The precision of E_C values obtained is high: the uncertainties, estimated by the effect on E_C of changing the Ewald convergence parameter η, ranged from ± 0.00002 eV to ± 0.01 eV in the worst case. The charge distribution for organic ions was idealized as fractional point charges localized at the crystallographic atomic positions: these charges were chosen from available theoretical and experimental estimates. The uncertainty in E_C due to different charge distribution models is typically ± 0.1 eV (± 3%): thus, even the simple Hückel model can give decent results.

E_C for Wurster's Blue Perchl orate is -4.1 eV/molecule: the crystal is stable under the binding provided by direct Coulomb interactions. E_C for N-Methylphenazinium Tetracyanoquino- dimethanide is 0.1 eV: exchange Coulomb interactions, which cannot be estimated classically, must provide the necessary binding.

EWALD was also used to test the McConnell classification of DA crystals. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine: 7,7,8,8-Tetracyanoquinodimethan) E_C = -4.0 eV while 2ε_o = 4.6₅ eV: clearly, exchange forces must provide the balance. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine:para-Chloranil) E_C = -4.4 eV, while 2ε_o = 5.0 eV: again EC falls short of 2ε₁. As a Gedankenexperiment, two nonionic crystals were assumed to be ionized: for (1:1)-(Hexamethyl- benzene:para-Chloranil) E_C = -4.5 eV, 2ε_o = 6.6 eV; for (1:1)- (Napthalene:Tetracyanoethylene) E_C = -4.3 eV, 2ε_o = 6.5 eV. Thus, exchange energies in these nonionic crystals must not exceed 1 eV.

Chapter III

A rapid-convergence quantum-mechanical formalism is derived to calculate the electronic energy of an arbitrary molecular (or molecular-ion) crystal: this provides estimates of crystal binding energies which include the exchange Coulomb inter- actions. Previously obtained LCAO-MO wavefunctions for the isolated molecule(s) ("unit cell spin-orbitals") provide the starting-point. Bloch's theorem is used to construct "crystal spin-orbitals". Overlap between the unit cell orbitals localized in different unit cells is neglected, or is eliminated by Löwdin orthogonalization. Then simple formulas for the total kinetic energy Q^(XT)_λ, nuclear attraction [λ/λ]^XT, direct Coulomb [λλ/λ'λ']^XT and exchange Coulomb [λλ'/λ'λ]^XT integrals are obtained, and direct-space brute-force expansions in atomic wavefunctions are given. Fourier series are obtained for [λ/λ]^XT, [λλ/λ'λ']^XT, and [λλ/λ'λ]^XT with the help of the convolution theorem; the Fourier coefficients require the evaluation of Silverstone's two-center Fourier transform integrals. If the short-range interactions are calculated by brute-force integrations in direct space, and the long-range effects are summed in Fourier space, then rapid convergence is possible for [λ/λ]^XT, [λλ/λ'λ']^XT and [λλ'/λ'λ]^XT. This is achieved, as in the Ewald method, by modifying each atomic wavefunction by a "Gaussian convergence acceleration factor", and evaluating separately in direct and in Fourier space appropriate portions of [λ/λ]^XT, etc., where some of the portions contain the Gaussian factor.

Quantum Mechanical Study on Tunnelling and Ballistic Transport of Nanometer Si MOSFETs

Using self-consistent calculations of million-atom Schrodinger-Poisson equations, we investigate the I-V characteristics of tunnelling and ballistic transport of nanometer metal oxide semiconductor field effect transistors (MOSFET) based on a full 3-D quantum mechanical simulation under nonequilibtium condition. Atomistic empirical pseudopotentials are used to describe the device Hamiltonian and the underlying bulk band structure. We find that the ballistic transport dominates the I-V characteristics, whereas the effects of tunnelling cannot be neglected with the maximal value up to 0.8mA/mu m when the channel length of MOSFET scales down to 25 nm. The effects of tunnelling transport lower the threshold voltage V-t. The ballistic current based on fully 3-D quantum mechanical simulation is relatively large and has small on-off ratio compared with results derived from the calculation methods of Luo et al.

A fully three-dimensional atomistic quantum mechanical study on random dopant-induced effects in 25-nm MOSFETs

A fully 3-D atomistic quantum mechanical simulation is presented to study the random dopant-induced effects in nanometer metal-oxide-semiconductor field-effect transistors. The empirical pseudopotential is used to represent the single particle Hamiltonian, and the linear combination of bulk band method is used to solve the million atom Schrodinger equation. The gate threshold fluctuation and lowering due to the discrete dopant configurations are studied. It is found that quantum mechanical effects increase the threshold fluctuation while decreasing the threshold lowering. The increase of threshold fluctuation is in agreement with the researchers' early study based on an approximated density gradient approach. However, the decrease in threshold lowering is in contrast with the density gradient calculations.

Effects of the transition dipole moment function on the dynamics of ozone photodissociation: an exact 3D quantum mechanical study

The effects of the transition dipole moment function (TDMF) on the dynamics Of O-3 photodissociation in the Hartley band have been exploited by means of exact 3D time-dependent wavepacket method using the SW potential energy surface [J. Chem. Phys. 78 (1983) 7191]. The calculations show that the explicit inclusion of the TDMF results in slight uniform reductions for the intensities of recurrence peaks of the autocorrelation function and a slight broadening of the absorption spectrum, in comparison with the result where the TDMF is assumed to be constant. The pattern of recurrence structures of the autocorrelation function is essentially unaffected. (C) 2001 Elsevier Science B.V. All rights reserved.

Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

Febrifugine derivative antimalarial activity: quantum mechanical predictors

O desenvolvimento de linhagens resistentes de Plasmodium falciparum tem encorajado a busca por novas drogas antimalariais. A febrifugina é uma substância natural com alta atividade contra o P. falciparum que apresenta propriedade emética e toxicidade para o fígado tal que não permitem o seu uso clínico. A busca por análogos que possam ter uma performance clínica melhor é um tema de pesquisa atual. Nosso objetivo é investigar a estrutura eletrônica teórica de uma família de derivados da febrifugina empregando cálculos semi-empíricos de orbitais moleculares, procurando por índices eletrônicos que possam ajudar a modelar novos derivados mais eficientes. Os resultados teóricos mostram que para as moléculas mais seletivas existe um agrupamento dos valores de determinados índices em intervalos bem definidos. O modelo proposto para se obter alta seletividade foi testado com sucesso.

Quantum mechanical study of YTiO3 to the investigation of piezoelectricity

In previous articles we reported through theoretical studies the piezoelectric effect in BaTiO3, SmTiO3, and YFeO3. In this paper, we used the Douglas-Kroll-Hess (DKH) second-order scalar relativistic method to investigate the piezoelectricity in YTiO3. In the calculations we used the [6s4p] and [10s5p4d] Gaussian basis sets for the O (3P) and Ti (5S) atoms, respectively, from the literature in combination with the (30s21p16d)/[15s9p6d] basis set for the Y ( 3D) atom, obtained by generator coordinate Hartree-Fock (GCHF) method, and they had their quality evaluated using calculations of total energy and orbital energies (HOMO and HOMO-1) of the 2TiO+1 and 1YO+1 fragments. The dipole moment, the total energy, and the total atomic charges in YTiO3 in C s space group were calculated. When we analyze those properties we verify that it is reasonable to believe that YTiO3 does not present piezoelectric properties. Copyright © 2011 Raimundo Dirceu de Paula Ferreira et al.

Quantum mechanical modeling of excited electronic states and their relationship to cathodoluminescence of BaZrO3

First-principles calculations set the comprehension over performance of novel cathodoluminescence (CL) properties of BaZrO3 prepared through microwave-assisted hydrothermal. Ground (singlet, s*) and excited (singlet s** and triplet t**) electronic states were built from zirconium displacement of 0.2 Å in {001} direction. Each ground and excited states were characterized by the correlation of their corresponding geometry with electronic structures and Raman vibrational frequencies which were also identified experimentally. A kind of optical polarization switching was identified by the redistribution of 4dz2 and 4dxz (Zr) orbitals and 2pz O orbital. As a consequence, asymmetric bending and stretching modes theoretically obtained reveal a direct dependence with their polyhedral intracluster and/or extracluster ZrO6 distortions with electronic structure. Then, CL of the as-synthesized BaZrO3 can be interpreted as a result of stable triplet excited states, which are able to trap electrons, delaying the emission process due to spin multiplicity changes. © 2013 AIP Publishing LLC.

Further Quantum Mechanical Evidence that Difluorotoluene does not Hydrogen Bond

Calculations were run on the methylated DNA base pairs adenine:thymine and adenine:difluorotoluene to further investigate the hydrogen-bonding properties of difluorotoluene (F). Geometries were optimized using hybrid density functional theory. Single-point calculations at the MP2(full) level were performed to obtain more rigorous energies. The functional counterpoise method was used to correct for the basis set superposition error (BSSE), and the interaction energies were also corrected for fragment relaxation. These corrections brought the B3LYP and MP2 interaction energies into excellent agreement. In the gas phase, the Gibbs free energies calculated at the B3LYP and MP2 levels of theory predict that A and T will spontaneously form an A:T pair while A:F spontaneously dissociates into A and F. Solvation effects on the pairing of the bases were explored using implicit solvent models for water and chloroform. In aqueous solution, both A:T and A:F are predicted to dissociate into their component monomers. Semiempirical calculations were performed on small sections of B-form DNA containing the two pairs, and the results provide support for the concept that base stacking is more important than hydrogen bonding for the stability of the A:F pair within a DNA helix.