Artificial molecular quantum rings under magnetic field influence

The ground states of a few electrons confined in two vertically coupled quantum rings in the presence of an external magnetic field are studied systematically within the current spin-density functional theory. Electron-electron interactions combined with inter-ring tunneling affect the electronic structure and the persistent current. For small values of the external magnetic field, we recover the zero magnetic field molecular quantum ring ground state configurations. Increasing the magnetic field many angular momentum, spin, and isospin transitions are predicted to occur in the ground state. We show that these transitions follow certain rules, which are governed by the parity of the number of electrons, the single-particle picture, Hund's rules, and many-body effects. (C) 2009 American Institute of Physics. [doi:10.1063/1.3223360]

Isotropic and anisotropic spin-spin interactions and a quantum phase transition in a dinuclear Cu(II) compound

We report electron-paramagnetic resonance (EPR) studies at similar to 9.5 GHz (X band) and similar to 34 GHz (Q band) of powder and single-crystal samples of the compound Cu(2)[TzTs](4) [N-thiazol-2-yl-toluenesulfonamidatecopper(II)], C(40)H(36)Cu(2)N(8)O(8)S(8), having copper(II) ions in dinuclear units. Our data allow determining an antiferromagnetic interaction J(0)=(-113 +/- 1) cm(-1) (H(ex)=-J(0)S(1)center dot S(2)) between Cu(II) ions in the dinuclear unit and the anisotropic contributions to the spin-spin coupling matrix D (H(ani)=S(1)center dot D center dot S(2)), a traceless symmetric matrix with principal values D/4=(0.198 +/- 0.003) cm(-1) and E/4=(0.001 +/- 0.003) cm(-1) arising from magnetic dipole-dipole and anisotropic exchange couplings within the units. In addition, the single-crystal EPR measurements allow detecting and estimating very weak exchange couplings between neighbor dinuclear units, with an estimated magnitude parallel to J(')parallel to=(0.060 +/- 0.015) cm(-1). The interactions between a dinuclear unit and the ""environment"" of similar units in the structure of the compound produce a spin dynamics that averages out the intradinuclear dipolar interactions. This coupling with the environment leads to decoherence, a quantum phase transition that collapses the dipolar interaction when the isotropic exchange coupling with neighbor dinuclear units equals the magnitude of the intradinuclear dipolar coupling. Our EPR experiments provide a new procedure to follow the classical exchange-narrowing process as a shift and collapse of the line structure (not only as a change of the resonance width), which is described with general (but otherwise simple) theories of magnetic resonance. Using complementary procedures, our EPR measurements in powder and single-crystal samples allow measuring simultaneously three types of interactions differing by more than three orders of magnitude (between 113 cm(-1) and 0.060 cm(-1)).

Entanglement in spatially inhomogeneous many-fermion systems

We investigate entanglement of strongly interacting fermions in spatially inhomogeneous environments. To quantify entanglement in the presence of spatial inhomogeneity, we propose a local-density approximation (LDA) to the entanglement entropy, and a nested LDA scheme to evaluate the entanglement entropy on inhomogeneous density profiles. These ideas are applied to models of electrons in superlattice structures with different modulation patterns, electrons in a metallic wire in the presence of impurities, and phase-separated states in harmonically confined many-fermion systems, such as electrons in quantum dots and atoms in optical traps. We find that the entanglement entropy of inhomogeneous systems is strikingly different from that of homogeneous systems.

Spin gaps and spin-flip energies in density-functional theory

Energy gaps are crucial aspects of the electronic structure of finite and extended systems. Whereas much is known about how to define and calculate charge gaps in density-functional theory (DFT), and about the relation between these gaps and derivative discontinuities of the exchange-correlation functional, much less is known about spin gaps. In this paper we give density-functional definitions of spin-conserving gaps, spin-flip gaps and the spin stiffness in terms of many-body energies and in terms of single-particle (Kohn-Sham) energies. Our definitions are as analogous as possible to those commonly made in the charge case, but important differences between spin and charge gaps emerge already on the single-particle level because unlike the fundamental charge gap spin gaps involve excited-state energies. Kohn-Sham and many-body spin gaps are predicted to differ, and the difference is related to derivative discontinuities that are similar to, but distinct from, those usually considered in the case of charge gaps. Both ensemble DFT and time-dependent DFT (TDDFT) can be used to calculate these spin discontinuities from a suitable functional. We illustrate our findings by evaluating our definitions for the Lithium atom, for which we calculate spin gaps and spin discontinuities by making use of near-exact Kohn-Sham eigenvalues and, independently, from the single-pole approximation to TDDFT. The many-body corrections to the Kohn-Sham spin gaps are found to be negative, i.e., single-particle calculations tend to overestimate spin gaps while they underestimate charge gaps.

Physical signatures of discontinuities of the time-dependent exchange-correlation potential

The exact exchange-correlation (XC) potential in time-dependent density-functional theory (TDDFT) is known to develop steps and discontinuities upon change of the particle number in spatially confined regions or isolated subsystems. We demonstrate that the self-interaction corrected adiabatic local-density approximation for the XC potential has this property, using the example of electron loss of a model quantum well system. We then study the influence of the XC potential discontinuity in a real-time simulation of a dissociation process of an asymmetric double quantum well system, and show that it dramatically affects the population of the resulting isolated single quantum wells. This indicates the importance of a proper account of the discontinuities in TDDFT descriptions of ionization, dissociation or charge transfer processes.

Density functional theory investigation of 3d, 4d, and 5d 13-atom metal clusters

The knowledge of the atomic structure of clusters composed by few atoms is a basic prerequisite to obtain insights into the mechanisms that determine their chemical and physical properties as a function of diameter, shape, surface termination, as well as to understand the mechanism of bulk formation. Due to the wide use of metal systems in our modern life, the accurate determination of the properties of 3d, 4d, and 5d metal clusters poses a huge problem for nanoscience. In this work, we report a density functional theory study of the atomic structure, binding energies, effective coordination numbers, average bond lengths, and magnetic properties of the 3d, 4d, and 5d metal (30 elements) clusters containing 13 atoms, M(13). First, a set of lowest-energy local minimum structures (as supported by vibrational analysis) were obtained by combining high-temperature first- principles molecular-dynamics simulation, structure crossover, and the selection of five well-known M(13) structures. Several new lower energy configurations were identified, e. g., Pd(13), W(13), Pt(13), etc., and previous known structures were confirmed by our calculations. Furthermore, the following trends were identified: (i) compact icosahedral-like forms at the beginning of each metal series, more opened structures such as hexagonal bilayerlike and double simple-cubic layers at the middle of each metal series, and structures with an increasing effective coordination number occur for large d states occupation. (ii) For Au(13), we found that spin-orbit coupling favors the three-dimensional (3D) structures, i.e., a 3D structure is about 0.10 eV lower in energy than the lowest energy known two-dimensional configuration. (iii) The magnetic exchange interactions play an important role for particular systems such as Fe, Cr, and Mn. (iv) The analysis of the binding energy and average bond lengths show a paraboliclike shape as a function of the occupation of the d states and hence, most of the properties can be explained by the chemistry picture of occupation of the bonding and antibonding states.

Theory of nitride oxide adsorption on transition metal (111) surfaces: a first-principles investigation

In this work, we report a density functional theory study of nitric oxide (NO) adsorption on close-packed transition metal (TM) Rh(111), Ir(111), Pd(111) and Pt(111) surfaces in terms of adsorption sites, binding mechanism and charge transfer at a coverage of Theta(NO) = 0.25, 0.50, 0.75 monolayer (ML). Based on our study, an unified picture for the interaction between NO and TM(111) and site preference is established, and valuable insights are obtained. At low coverage (0.25 ML), we find that the interaction of NO/TM(111) is determined by an electron donation and back-donation process via the interplay between NO 5 sigma/2 pi* and TM d-bands. The extent of the donation and back-donation depends critically on the coordination number (adsorption sites) and TM d-band filling, and plays an essential role for NO adsorption on TM surfaces. DFT calculations shows that for TMs with high d-band filling such as Pd and Pt, hollow-site NO is energetically the most favorable, and top-site NO prefers to tilt away from the normal direction. While for TMs with low d-band filling (Rh and Ir), top-site NO perpendicular to the surfaces is energetically most favorable. Electronic structure analysis show that irrespective of the TM and adsorption site, there is a net charge transfer from the substrate to the adsorbate due to overwhelming back-donation from the TM substrate to the adsorbed NO molecules. The adsorption-induced change of the work function with respect to bare surfaces and dipole moment is however site dependent, and the work function increases for hollow-site NO, but decreases for top-site NO, because of differences in the charge redistribution. The interplay between the energetics, lateral interaction and charge transfer, which is element dependent, rationalizes the structural evolution of NO adsorption on TM(111) surfaces in the submonolayer regime.

Transition-metal 13-atom clusters assessed with solid and surface-biased functionals

First-principles density-functional theory studies have reported open structures based on the formation of double simple-cubic (DSC) arrangements for Ru(13), Rh(13), Os(13), and Ir(13), which can be considered an unexpected result as those elements crystallize in compact bulk structures such as the face-centered cubic and hexagonal close-packed lattices. In this work, we investigated with the projected augmented wave method the dependence of the lowest-energy structure on the local and semilocal exchange-correlation (xc) energy functionals employed in density-functional theory. We found that the local-density approximation (LDA) and generalized-gradient formulations with different treatment of the electronic inhomogeneities (PBE, PBEsol, and AM05) confirm the DSC configuration as the lowest-energy structure for the studied TM(13) clusters. A good agreement in the relative total energies are obtained even for structures with small energy differences, e. g., 0.10 eV. The employed xc functionals yield the same total magnetic moment for a given structure, i.e., the differences in the bond lengths do not affect the moments, which can be attributed to the atomic character of those clusters. Thus, at least for those systems, the differences among the LDA, PBE, PBEsol, and AM05 functionals are not large enough to yield qualitatively different results. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3577999]

Electron theoretical investigation of the stability of the B2-TiFe compound

The metastable phase diagram of the BCC-based ordering equilibria in the Ti-Fe system has been calculated using a truncated cluster expansion, through the combination of FP-LAPW and cluster variation method (CVM) in the irregular tetrahedron cluster approximation. The results are compared with phenomenological CVM assessments of the system and suggest that the value for the experimental formation enthalpy of the B2-TiFe compound should be significantly more negative than the currently assessed value. (C) 2008 Elsevier B.V. All rights reserved.

Thermodynamic self-consistency issues related to the Cluster Variation Method: The case of the BCC Cr-Fe (Chromium-Iron) system

The Cluster Variation Method (CVM), introduced over 50 years ago by Prof. Dr. Ryoichi Kikuchi, is applied to the thermodynamic modeling of the BCC Cr-Fe system in the irregular tetrahedron approximation, using experimental thermochemical data as initial input for accessing the model parameters. The results are checked against independent data on the low-temperature miscibility gap, using increasingly accurate thermodynamic models, first by the inclusion of the magnetic degrees of freedom of iron and then also by the inclusion of the magnetic degrees of freedom of chromium. It is shown that a reasonably accurate description of the phase diagram at the iron-rich side (i.e. the miscibility gap borders and the Curie line) is obtained, but only at expense of the agreement with the above mentioned thermochemical data. Reasons for these inconsistencies are discussed, especially with regard to the need of introducing vibrational degrees of freedom in the CVM model. (C) 2008 Elsevier Ltd. All rights reserved.

N-mesityl-C-acylketenimines: 1,5-sigmatropic shifts and electrocyclization to quinolines

Flash vacuum thermolysis (FVT) of triazoles 6a-c generates alpha-oxoketenimines 10, the ester 10a being isolable. FVT of pyrroledione 8 generates the isomeric imidoylketene 9a. Ketenes 9 and ketenimines 10 undergo thermal interconversion by 1,3-shifts of methoxy and dimethylamino groups under mild FVT conditions (ca. 350-400 degrees C). Both 9 and 10 are directly observable by IR spectroscopy at either 77 K or on Ar matrix isolation at 12 K. On FVT at temperatures above ca. 400 degrees C, the ketenimines 10 undergo a 1,5-H shift to o-quinoid imines 12/13, followed by electrocyclization to dihydroquinolines 14 (unobserved) and 15 (observed by NMR). The latter are easily oxidized to alkylquinoline-3-carboxylates or quinoline-3-carboxamides 16 by atmospheric oxygen. Ab initio calculations on model compounds 18-23 predict an energy barrier of ca. 38 kcal mol(-1) (161 kJ mol(-1)) for the 1,5-H shift in N-(o-methylphenyl)ketenimines via the transition state TS19 followed by an electrocyclization barrier to dihydroquinoline 23a via TS22a of ca. 16 kcal mol(-1).

Transport properties of strongly correlated metals: A dynamical mean-field approach

The temperature dependence of the transport properties of the metallic phase of a frustrated Hubbard model on the hypercubic lattice at half-filling is calculated. Dynamical mean-held theory, which maps the Hubbard model onto a single impurity,Anderson model that is solved self-consistently, and becomes exact in the limit of large dimensionality, is used. As the temperature increases there is a smooth crossover from coherent Fermi liquid excitations at low temperatures to incoherent excitations at high temperatures. This crossover leads to a nonmonotonic temperature dependence for the resistance, thermopower, and Hall coefficient, unlike in conventional metals. The resistance smoothly increases from a quadratic temperature dependence at low temperatures to large values which can exceed the Mott-Ioffe-Regel value ha/e(2) (where a is a lattice constant) associated with mean free paths less than a lattice constant. Further signatures of the thermal destruction of quasiparticle excitations are a peak in the thermopower and the absence of a Drude peak in the optical conductivity. The results presented here are relevant to a wide range of strongly correlated metals, including transition metal oxides, strontium ruthenates, and organic metals.

Cyclotron effective masses in layered metals

Many layered metals such as quasi-two-dimensional organic molecular crystals show properties consistent with a Fermi-liquid description at low temperatures. The effective masses extracted from the temperature dependence of the magnetic oscillations observed in these materials are in the range, m(c)*/m(e) similar to 1 - 7, suggesting that these systems are strongly correlated. However, the ratio m(c)*/m(e) contains both the renormalization due to the electron-electron interaction and the periodic potential of the lattice. We show that for any quasi-two-dimensional band structure, the cyclotron mass is proportional to the density-of-states at the Fermi energy. Due to Luttinger's theorem, this result is also valid in the presence of interactions. We then evaluate m(c) for several model band structures for the beta, kappa, and theta families of (BEDT-TTF)(2)X, where BEDT-TTF is bis-(ethylenedithia-tetrathiafulvalene) and X is an anion. We find that for kappa-(BEDT-TTF)(2)X, the cyclotron mass of the beta orbit, m(c)*(beta) is close to 2 m(c)*(alpha), where m(c)*(alpha) is the effective mass of the alpha orbit. This result is fairly insensitive to the band-structure details. For a wide range of materials we compare values of the cyclotron mass deduced from band-structure calculations to values deduced from measurements of magnetic oscillations and the specific-heat coefficient gamma.

A Computational Study of Tetrafluoro-[2.2]Cyclophanes

A computational study of the isomers of tetrafluorinated [2.2]cyclophanes persubstituted in one ring, namely F-4-[2.2]paracyclophane (4), F-4-anti-[2.2]metacyclophane (5a), F-4-syn-[2.2]metacyclophane (5b), and F-4-[2.2]metaparacyclophane (6a and 6b), was carried out. The effects of fluorination on the geometries, relative energies, local and global aromaticity, and strain energies of the bridges and rings were investigated. An analysis of the electron density by B3PW91/6-31+G(d,p), B3LYP/6-31+G(d,p), and MP2/6-31+G(d,p) was carried out using the natural bond orbitals (NBO), natural steric analysis (NSA), and atoms in molecules (AIM) methods. The analysis of frontier molecular orbitals (MOs) was also employed. The results indicated that the molecular structure of [2.2]paracyclophane is the most affected by the fluorination. Isodesmic reactions showed that the fluorinated rings are more strained than the nonfluorinated ones. The NICS, HOMA, and PDI criteria evidenced that the fluorination affects the aromaticity of both the fluorinated and the nonfluorinated rings. The NBO and NSA analyses gave an indication that the fluorination increases not only the number of through-space interactions but also their magnitude. The AIM analysis suggested that the through-space interactions are restricted to the F-4-[2.2]metacyclophanes. In addition, the atomic properties, computed over the atomic basins, shave evidence that not only the substitution, but also the position of the bridges could affect the atomic charges. the first atomic moments, and the atomic volumes.

Evaluation of the enthalpy of formation, proton affinity, and gas-phase basicity of gamma-butyrolactone and 2-pyrrolidinone by isodesmic reactions

The knowledge of thermochemical parameters such as the enthalpy of formation, gas-phase basicity, and proton affinity may be the key to understanding molecular reactivity. The obtention of these thermochemical parameters by theoretical chemical models may be advantageous when experimental measurements are difficult to accomplish. The development of ab initio composite models represents a major advance in the obtention of these thermochemical parameters,. but these methods do not always lead to accurate values. Aiming at achieving a comparison between the ab initio models and the hybrid models based on the density functional theory (DFT), we have studied gamma-butyrolactone and 2-pyrrolidinone with a goal of obtaining high-quality thermochemical parameters using the composite chemical models G2, G2MP2, MP2, G3, CBS-Q, CBS-4, and CBS-QB3; the DFT methods B3LYP, B3P86, PW91PW91, mPW1PW, and B98; and the basis sets 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d), 6-311+G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ. Values obtained for the enthalpies of formation, proton affinity, and gas-phase basicity of the two target molecules were compared to the experimental data reported in the literature. The best results were achieved with the use of DFT models, and the B3LYP method led to the most accurate data.