Study of the oxygen vacancy influence on magnetic properties of Fe- and Co-doped SnO2 diluted alloys

Transition-metal (TM)-doped diluted magnetic oxides (DMOs) have attracted attention from both experimental and theoretical points of view due to their potential use in spintronics towards new nanostructured devices and new technologies. In the present work, we study the magnetic properties of Sn0.96TM0.04O2 and Sn0.96TM0.04O1.98(V (O))(0.02), where TM = Fe and Co, focusing in particular in the role played by the presence of O vacancies nearby the TM. The calculated total energy as a function of the total magnetic moment per cell shows a magnetic metastability, corresponding to a ground state, respectively, with 2 and 1 mu(B)/cell, for Fe and Co. Two metastable states, with 0 and 4 mu(B)/cell were found for Fe, and a single value, 3 mu(B)/cell, for Co. The spin-crossover energies (E (S)) were calculated. The values are E (S) (0/2) = 107 meV and E (S) (4/2) = 25 meV for Fe. For Co, E (S) (3/1) = 36 meV. By creating O vacancies close to the TM site, we show that the metastablity and E (S) change. For iron, a new state appears, and the state with zero magnetic moment disappears. The ground state is 4 mu(B)/cell instead of 2 mu(B)/cell, and the energy E (S) (2/4) is 30 meV. For cobalt, the ground state is then found with 3 mu(B)/cell and the metastable state with 1 mu(B)/cell. The spin-crossover energy E (S) (1/3) is 21 meV. Our results suggest that these materials may be used in devices for spintronic applications that require different magnetization states.

Hydrogen bond and the resonance effect on the formamide-water complexes

The interaction of formamide and the two transition states of its amide group rotation with one, two, or three water molecules was studied in vacuum. Great differences between the electronic structure of formamide in its most stable form and the electronic structure of the transition states were noticed. Intermolecular interactions were intense, especially in the cases where the solvent interacted with the amide and the carbonyl groups simultaneously. In the transition states, the interaction between the lone pair of nitrogen and the water molecule becomes important. With the aid of the natural bond orbitals, natural resonance theory, and electron localization function (ELF) analyses an increase in the resonance of planar formamide with the addition of successive water molecules was observed. Such observation suggests that the hydrogen bonds in the formamidewater complexes may have some covalent character. These results are also supported by the quantitative ELF analyses. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012

Site localization of Cd impurities in sapphire

By combining first-principles electronic structure calculations and existing time-differential gamma-gamma perturbed-angular-correlation experiments we studied the site localization, the local environment, and the electronic structure of Cd impurities in sapphire (alpha-Al2O3) single crystals in different charged states. The ab initio calculations were performed with the full-potential augmented plane wave plus local orbitals method and the projector augmented wave method. Comparing the calculated electric-field-gradient tensor at the Cd nuclei in the alpha-Al2O3 host lattice and the corresponding available experimental values, we have seen that it is equally possible for Cd to replace an Al atom (in a negative charge state) or to be placed in an interstitial site (in a neutral charge state). To finally address the issue of the Cd impurity localization, we performed formation energy calculations. These results have shown that Cd placed in the substitutional Al site, in the negatively charged state, is the most probable configuration.

CdSe/CdTe interface band gaps and band offsets calculated using spin-orbit and self-energy corrections

We performed ab initio calculations of the electronic structures of bulk CdSe and CdTe, and their interface band alignments on the CdSe in-plane lattice parameters. For this, we employed the LDA-1/2 self-energy correction scheme [L.G. Ferreira, M. Marques, L.K. Teles, Phys. Rev. B 78 (2008) 125116] to obtain corrected band gaps and band offsets. Our calculations include the spin-orbit effects for the bulk cases, which have shown to be of importance for the equilibrium systems and are possibly degraded in these strained semiconductors. Therefore, the SO showed reduced importance for the band alignment of this particular system. Moreover, the electronic structure calculated along the transition region across the CdSe/CdTe interface shows an interesting non-monotonic variation of the band gap in the range 0.8-1.8 eV, which may enhance the absorption of light for corresponding frequencies at the interface between these two materials in photovoltaic applications. (C) 2012 Elsevier B.V. All rights reserved.

Theoretical study of the XP3 (X = Al, B, Ga) clusters

The lowest singlet and triplet states of AlP3, GaP3 and BP3 molecules with C-s, C-2v and C-3v symmetries were characterized using the B3LYP functional and the aug-cc-pVTZ and aug-cc-pVQZ correlated consistent basis sets. Geometrical parameters and vibrational frequencies were calculated and compared to existent experimental and theoretical data. Relative energies were obtained with single point CCSD(T) calculations using the aug-cc-pVTZ, aug-cc-pVQZ and aug-cc-pV5Z basis sets, and then extrapolating to the complete basis set (CBS) limit. (C) 2011 Elsevier B.V. All rights reserved.

A theoretical study of SnF2+, SnCl2+, and SnO2+ and their experimental search

We present a detailed theoretical study of the stability of the gas-phase diatomic dications SnF2+, SnCl2+, and SnO2+ using ab initio computer calculations. The ground states of SnF2+, SnCl2+, and SnO2+ are thermodynamically stable, respectively, with dissociation energies of 0.45, 0.30, and 0.42 eV. Whereas SnF2+ dissociates into Sn2+ + F, the long range behaviour of the potential energy curves of SnCl2+ and SnO2+ is repulsive and wide barrier heights due to avoided crossing act as a kind of effective dissociation energy. Their equilibrium internuclear distances are 4.855, 5.201, and 4.852 a(0), respectively. The double ionisation energies (T-e) to form SnF2+, SnCl2+, and SnO2+ from their respective neutral parents are 25.87, 23.71, and 25.97 eV. We combine our theoretical work with the experimental results of a search for these doubly positively charged diatomic molecules in the gas phase. SnO2+ and SnF2+ have been observed for prolonged oxygen (O-16(-)) ion beam sputtering of a tin metal foil and of tin (II) fluoride (SnF2) powder, respectively, for ion flight times of about 10(-5) s through a magnetic-sector mass spectrometer. In addition, SnCl2+ has been detected for O-16(-) ion surface bombardment of stannous (tin (II)) chloride (SnCl2) powder. To our knowledge, SnF2+ is a novel gas-phase molecule, whereas SnCl2+ had been detected previously by electron-impact ionization mass spectrometry, and SnO2+ had been observed before by spark source mass spectrometry as well as by atom probe mass spectrometry. We are not aware of any previous theoretical studies of these molecular systems. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4758475]

Digital magnetic heterostructures based on GaN using GGA-1/2 approach

We present ab-initio calculations of seven digital magnetic heterostructures, GaN delta-doped with V, Cr, Mn, Fe, Co, Ni, and Cu, forming two-dimensional systems. Only GaN delta-doped with V or Cr present a ferromagnetic ground state with high Curie temperatures. For both, to better describe the electronic properties, we used the GGA-1/2 approach. The ground state of GaN/Cr resulted in a two dimensional half-metal, with 100% spin polarization. For GaN/V, we obtained an insulating state: integer magnetic moment of 2.0 mu(B), a minority spin gap of 3.0 eV close to the gap of GaN, but a majority spin gap of 0.34 eV. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4751285]

Study of the Influence of Localized Vibrational Modes in Charge Transport Properties at Nanoscale Systems.

In molecular and atomic devices the interaction between electrons and ionic vibrations has an important role in electronic transport. The electron-phonon coupling can cause the loss of the electron's phase coherence, the opening of new conductance channels and the suppression of purely elastic ones. From the technological viewpoint phonons might restrict the efficiency of electronic devices by energy dissipation, causing heating, power loss and instability. The state of the art in electron transport calculations consists in combining ab initio calculations via Density Functional Theory (DFT) with Non-Equilibrium Green's Function formalism (NEGF). In order to include electron-phonon interactions, one needs in principle to include a self-energy scattering term in the open system Hamiltonian which takes into account the effect of the phonons over the electrons and vice versa. Nevertheless this term could be obtained approximately by perturbative methods. In the First Born Approximation one considers only the first order terms of the electronic Green's function expansion. In the Self-Consistent Born Approximation, the interaction self-energy is calculated with the perturbed electronic Green's function in a self-consistent way. In this work we describe how to incorporate the electron-phonon interaction to the SMEAGOL program (Spin and Molecular Electronics in Atomically Generated Orbital Landscapes), an ab initio code for electronic transport based on the combination of DFT + NEGF. This provides a tool for calculating the transport properties of materials' specific system, particularly in molecular electronics. Preliminary results will be presented, showing the effects produced by considering the electron-phonon interaction in nanoscale devices.

Theoretical and experimental studies on the electronic, optical, and structural properties of poly-pyrrole-2-carboxylic acid films

A theoretical approach is used here to explain experimental results obtained from the electrosynthesis of polypyrrole-2-carboxylic acid (PPY-2-COOH) films in nonaqueous medium. An analysis of the Fukui function (reactivity index) indicates that the monomer (pyrrole-2-carboxylic acid, PY-2-COOH), and dimers and trimers are oxidized in the C4 or C5 positions of the heterocyclic ring of the PY-2-COOH structure. After calculating the heat of formation using semiempirical Austin Model 1 post-Hartree-Fock parameterization for dimer species, both C4 and C5 positions adjacent to the aromatic rings of PPY-2-COOH were considered the most susceptible ones to oxidative coupling reactions. The ZINDO-S/CI semiempirical method was used to simulate the electronic transitions typically seen in the UV-VIS-NIR range in monomer and oligomers with different conjugation lengths. The use of an electrochemical quartz crystal microbalance provides sufficient information to propose a polymerization mechanism of PY-2-COOH based on molecular modeling and experimental results.

Development of empirical potentials for amorphous silica

Development of empirical potentials for amorphous silica Amorphous silica (SiO2) is of great importance in geoscience and mineralogy as well as a raw material in glass industry. Its structure is characterized as a disordered continuous network of SiO4 tetrahedra. Many efforts have been undertaken to understand the microscopic properties of silica by classical molecular dynamics (MD) simulations. In this method the interatomic interactions are modeled by an effective potential that does not take explicitely into account the electronic degrees of freedom. In this work, we propose a new methodology to parameterize such a potential for silica using ab initio simulations, namely Car-Parrinello (CP) method [Phys. Rev. Lett. 55, 2471 (1985)]. The new potential proposed is compared to the BKS potential [Phys. Rev. Lett. 64, 1955 (1990)] that is considered as the benchmark potential for silica. First, CP simulations have been performed on a liquid silica sample at 3600 K. The structural features so obtained have been compared to the ones predicted by the classical BKS potential. Regarding the bond lengths the BKS tends to underestimate the Si-O bond whereas the Si-Si bond is overestimated. The inter-tetrahedral angular distribution functions are also not well described by the BKS potential. The corresponding mean value of theSiOSi angle is found to be ≃ 147◦, while the CP yields to aSiOSi angle centered around 135◦. Our aim is to fit a classical Born-Mayer/Coulomb pair potential using ab initio calculations. To this end, we use the force-matching method proposed by Ercolessi and Adams [Europhys. Lett. 26, 583 (1994)]. The CP configurations and their corresponding interatomic forces have been considered for a least square fitting procedure. The classical MD simulations with the resulting potential have lead to a structure that is very different from the CP one. Therefore, a different fitting criterion based on the CP partial pair correlation functions was applied. Using this approach the resulting potential shows a better agreement with the CP data than the BKS ones: pair correlation functions, angular distribution functions, structure factors, density of states and pressure/density were improved. At low temperature, the diffusion coefficients appear to be three times higher than those predicted by the BKS model, however showing a similar temperature dependence. Calculations have also been carried out on crystalline samples in order to check the transferability of the potential. The equilibrium geometry as well as the elastic constants of α-quartz at 0 K are well described by our new potential although the crystalline phases have not been considered for the parameterization. We have developed a new potential for silica which represents an improvement over the pair potentials class proposed so far. Furthermore, the fitting methodology that has been developed in this work can be applied to other network forming systems such as germania as well as mixtures of SiO2 with other oxides (e.g. Al2O3, K2O, Na2O).

Electron spectroscopy of novel charge transfer systems based on polycyclic aromatic hydrocarbons

Organische Ladungstransfersysteme weisen eine Vielfalt von konkurrierenden Wechselwirkungen zwischen Ladungs-, Spin- und Gitterfreiheitsgraden auf. Dies führt zu interessanten physikalischen Eigenschaften, wie metallische Leitfähigkeit, Supraleitung und Magnetismus. Diese Dissertation beschäftigt sich mit der elektronischen Struktur von organischen Ladungstransfersalzen aus drei Material-Familien. Dabei kamen unterschiedliche Photoemissions- und Röntgenspektroskopietechniken zum Einsatz. Die untersuchten Moleküle wurden z.T. im MPI für Polymerforschung synthetisiert. Sie stammen aus der Familie der Coronene (Donor Hexamethoxycoronen HMC und Akzeptor Coronen-hexaon COHON) und Pyrene (Donor Tetra- und Hexamethoxypyren TMP und HMP) im Komplex mit dem klassischen starken Akzeptor Tetracyanoquinodimethan (TCNQ). Als dritte Familie wurden Ladungstransfersalze der k-(BEDT-TTF)2X Familie (X ist ein monovalentes Anion) untersucht. Diese Materialien liegen nahe bei einem Bandbreite-kontrollierten Mottübergang im Phasendiagramm.rnFür Untersuchungen mittels Ultraviolett-Photoelektronenspektroskopie (UPS) wurden UHV-deponierte dünne Filme erzeugt. Dabei kam ein neuer Doppelverdampfer zum Einsatz, welcher speziell für Milligramm-Materialmengen entwickelt wurde. Diese Methode wies im Ladungstransferkomplex im Vergleich mit der reinen Donor- und Akzeptorspezies energetische Verschiebungen von Valenzzuständen im Bereich weniger 100meV nach. Ein wichtiger Aspekt der UPS-Messungen lag im direkten Vergleich mit ab-initio Rechnungen.rnDas Problem der unvermeidbaren Oberflächenverunreinigungen von lösungsgezüchteten 3D-Kristallen wurde durch die Methode Hard-X-ray Photoelectron Spectroscopy (HAXPES) bei Photonenenergien um 6 keV (am Elektronenspeicherring PETRA III in Hamburg) überwunden. Die große mittlere freie Weglänge der Photoelektronen im Bereich von 15 nm resultiert in echter Volumensensitivität. Die ersten HAXPES Experimente an Ladungstransferkomplexen weltweit zeigten große chemische Verschiebungen (mehrere eV). In der Verbindung HMPx-TCNQy ist die N1s-Linie ein Fingerabdruck der Cyanogruppe im TCNQ und zeigt eine Aufspaltung und einen Shift zu höheren Bindungsenergien von bis zu 6 eV mit zunehmendem HMP-Gehalt. Umgekehrt ist die O1s-Linie ein Fingerabdruck der Methoxygruppe in HMP und zeigt eine markante Aufspaltung und eine Verschiebung zu geringeren Bindungsenergien (bis zu etwa 2,5eV chemischer Verschiebung), d.h. eine Größenordnung größer als die im Valenzbereich.rnAls weitere synchrotronstrahlungsbasierte Technik wurde Near-Edge-X-ray-Absorption Fine Structure (NEXAFS) Spektroskopie am Speicherring ANKA Karlsruhe intensiv genutzt. Die mittlere freie Weglänge der niederenergetischen Sekundärelektronen (um 5 nm). Starke Intensitätsvariationen von bestimmten Vorkanten-Resonanzen (als Signatur der unbesetzte Zustandsdichte) zeigen unmittelbar die Änderung der Besetzungszahlen der beteiligten Orbitale in der unmittelbaren Umgebung des angeregten Atoms. Damit war es möglich, präzise die Beteiligung spezifischer Orbitale im Ladungstransfermechanismus nachzuweisen. Im genannten Komplex wird Ladung von den Methoxy-Orbitalen 2e(Pi*) und 6a1(σ*) zu den Cyano-Orbitalen b3g und au(Pi*) und – in geringerem Maße – zum b1g und b2u(σ*) der Cyanogruppe transferiert. Zusätzlich treten kleine energetische Shifts mit unterschiedlichem Vorzeichen für die Donor- und Akzeptor-Resonanzen auf, vergleichbar mit den in UPS beobachteten Shifts.rn

Untersuchung der Kernstruktur von 3 He mittels Polarisationsobservablen

Durch die Möglichkeit, gleichzeitig mehrere Polarisationsfreiheitsgradernin der quasi-elastischen Elektronstreuung an $^3mathrm{He}$ zurnmessen, bietet sich ein neuer experimenteller Zugang zu kleinen, aber rnwichtigen Partialwellenbeiträgen ($S'$, $D$-Welle) desrn$^3mathrm{He}$-Grundzustands. Dies ermöglicht nicht nur ein tieferesrnVerständnis des Drei-Körper-Systems, sondern bietet auch diernMöglichkeit, Erkenntnisse über die $^3mathrm{He}$-Struktur undrnDynamik zu erlangen. Mit Hilfe dieser Informationen lassen sich abrninitio Rechnungen testen, sowie Korrekturen berechnen, die für anderernExperimente (z.B. Messung von $G_{en}$) benötigt werden. rnrnModerne Faddeev-Rechnungen liefern nicht nur eine quantitativernBeschreibung des $^3mathrm{He}$-Grundzustands, sondern geben auchrneinen Einblick in die sogenannten spinabhängigenrnImpulsverteilungen. Eine gründliche experimentelle Untersuchung ist in rndiesem Zusammenhang nötig, um eine solide Basis für die Üperprüfungrnder theoretische Modelle zu liefern. EinrnDreifach-Polarisationsexperiment liefert hier zum einen wichtigernDaten, zum anderen kann damit untersucht werden, ob mit der Methoderndes glqq Deuteron-Tagginggrqq polarisiertes $^3mathrm{He}$ alsrneffektives polarisiertes Protonentarget verwendet werden kann. rnrnDas hier vorgestellte Experiment kombiniert erstmals Strahl- undrnTargetpolarisation sowie die Messung der Polarisation des auslaufendenrnProtons. Das Experiment wurde im Sommer 2007 an derrnDrei-Spektrometer-Anlage der A1-Kollaboration an MAMI rndurchgeführt. Dabei wurde mit einer Strahlenergie vonrn$E=855,mathrm{MeV}$ bei $q^2=-0,14,(mathrm{GeV/c})^2$rn$(omega=0,13,mathrm{GeV}$, $q=0,4,mathrm{GeV/c})$ gemessen.rnrnDie bestimmten Wirkungsquerschnitte, sowie die Strahl-Target- und diernDreifach-Asymmetrie werden mit theoretischen Modellrechnungen vonrnJ. Golak (Plane Wave Impuls Approximation PWIA, sowie ein Modell mitrnEndzustandswechselwirkung) verglichen. Zudem wurde das Modell von dernForest verwendet, welches den Wirkungsquerschnitt über eine gemessenernSpektralfunktion berechnet. Der Vergleich mit den Modellrechnungenrnzeigt, dass sowohl der Wirkungsquerschnitt, als auch die Doppel- undrnDreifach-Asymmetrie gut mit den theoretischen Rechnungenrnübereinstimmen. rnrnDie Ergebnisse dieser Arbeit bestätigen, dass polarisiertesrn$^3mathrm{He}$ nicht nur als polarisiertes Neutronentarget, sondernrndurch Nachweis des Deuterons ebenfalls als polarisiertesrnProtonentarget verwendet werden kann.

Electronic structure and physical properties of Heusler compounds for thermoelectric and spintronic applications

This thesis focuses on synthesis as well as investigations of the electronic structure and properties of Heusler compounds for spintronic and thermoelectric applications.rnThe first part reports on the electronic and crystal structure as well as the mechanical, magnetic, and transport properties of the polycrystalline Heusler compound Co2MnGe. The crystalline structure was examined in detail by extended X-ray absorption fine structure spectroscopy and anomalous X-ray diffraction. The low-temperature magnetic moment agrees well with the Slater-Pauling rule and indicates a half-metallic ferromagnetic state of the compound, as is predicted by ab-initio calculations. Transport measurements and hard X-ray photoelectron spectroscopy (HAXPES) were performed to explain the electronic structure of the compound.rnA major part of the thesis deals with a systematical investigation of Heusler compounds for thermoelectric applications. Few studies have been reported on thermoelectric properties of p-type Heusler compounds. Therefore, this thesis focuses on the search for new p-type Heusler compounds with high thermoelectric efficiency. The substitutional series NiTi1−xMxSn and CoTi1−xMxSb (where M = Sc, V and 0 ≤ x ≤ 0.2) were synthesized and investigated theoretically and experimentally with respect to electronic structure and transport properties. The results show the possibility to create n-type and p-type thermoelectrics within one Heusler compound. The pure compounds showed n-type behavior, while under Sc substitution the system switched to p-type behavior. A maximum Seebeck coefficient of +230 μV/K (at 350 K) was obtained for NiTi0.26Sc0.04Zr0.35Hf0.35Sn, which is one of the highest values for p-type thermoelectric compounds based on Heusler alloys up to now. HAXPES valence band measurement show massive in gap states for the parent compounds NiTiSn, CoTiSb and NiTi0.3Zr0.35Hf0.35Sn. This proves that the electronic states close to the Fermi energy play a key role for the behavior of the transport properties. Furthermore, the electronic structure of the gapless Heusler compounds PtYSb, PtLaBi and PtLuSb were investigated by bulk sensitive HAXPES. The linear behavior of the spectra close to εF proves the bulk origin of Dirac-cone type density of states. Furthermore, a systematic study on the optical and transport properties of PtYSb is presented. The compound exhibits promising thermoelectric properties with a high figure of merit (ZT = 0.2) and a Hall mobility μh of 300 cm2/Vs at 350 K.rnThe last part of this thesis describes the linear dichroism in angular-resolved photoemission from the valence band of NiTi0.9Sc0.1Sn and NiMnSb. High resolution photoelectron spectroscopy was performed with an excitation energy of hν = 7.938 keV. The linear polarization of the photons was changed using an in-vacuum diamond phase retarder. Noticeable linear dichroism is found in the valence bands and this allows for a symmetry analysis of the contributing states. The differences in the spectra are found to be caused by symmetry dependent angular asymmetry parameters, and these occur even in polycrystalline samples without preferential crystallographic orientation.rnIn summary, Heusler compounds with 1:1:1 and 2:1:1 stoichiometry were synthesized and examined by chemical and physical methods. Overall, this thesis shows that the combination of first-principle calculations, transport measurements and high resolution high energy photoelectron spectroscopy analysis is a very powerful tool for the design and development of new materials for a wide range of applications from spintronic applications to thermoelectric applications.rn

Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)n where n = 1–10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.

Hydrogen bonding of nucleotide base pairs: Application of the PM3 method

The ability of the pm3 semiempirical quantum mechanical method to reproduce hydrogen bonding in nucleotide base pairs was assessed. Results of pm3 calculations on the nucleotides 2′-deoxyadenosine 5′-monophosphate (pdA), 2′-deoxyguanosine 5′-monophosphate (pdG), 2′-deoxycytidine 5′-monophosphate (pdC), and 2′-deoxythymidine 5′-monophosphate (pdT) and the base pairs pdA–pdT, pdG–pdC, and pdG(syn)–pdC are presented and discussed. The pm3 method is the first of the parameterized nddo quantum mechanical models with any ability to reproduce hydrogen bonding between nucleotide base pairs. Intermolecular hydrogen bond lengths between nucleotides displaying Watson–Crick base pairing are 0.1–0.2 Å less than experimental results. Nucleotide bond distances, bond angles, and torsion angles about the glycosyl bond (χ), the C4′C5′ bond (γ), and the C5′O5′ bond (β) agree with experimental results. There are many possible conformations of nucleotides. pm3 calculations reveal that many of the most stable conformations are stabilized by intramolecular CHO hydrogen bonds. These interactions disrupt the usual sugar puckering. The stacking interactions of a dT–pdA duplex are examined at different levels of gradient optimization. The intramolecular hydrogen bonds found in the nucleotide base pairs disappear in the duplex, as a result of the additional constraints on the phosphate group when part of a DNA backbone. Sugar puckering is reproduced by the pm3 method for the four bases in the dT–pdA duplex. pm3 underestimates the attractive stacking interactions of base pairs in a B-DNA helical conformation. The performance of the pm3 method implemented in SPARTAN is contrasted with that implemented in MOPAC. At present, accurate ab initio calculations are too timeconsuming to be of practical use, and molecular mechanics methods cannot be used to determine quantum mechanical properties such as reaction-path calculations, transition-state structures, and activation energies. The pm3 method should be used with extreme caution for examination of small DNA systems. Future parameterizations of semiempirical methods should incorporate base stacking interactions into the parameterization data set to enhance the ability of these methods.