Reformulating time-dependent density functional theory with non-orthogonal localized molecular orbitals.

Time-dependent density functional theory (TDDFT) has broad application in the study of electronic response, excitation and transport. To extend such application to large and complex systems, we develop a reformulation of TDDFT equations in terms of non-orthogonal localized molecular orbitals (NOLMOs). NOLMO is the most localized representation of electronic degrees of freedom and has been used in ground state calculations. In atomic orbital (AO) representation, the sparsity of NOLMO is transferred to the coefficient matrix of molecular orbitals (MOs). Its novel use in TDDFT here leads to a very simple form of time propagation equations which can be solved with linear-scaling effort. We have tested the method for several long-chain saturated and conjugated molecular systems within the self-consistent charge density-functional tight-binding method (SCC-DFTB) and demonstrated its accuracy. This opens up pathways for TDDFT applications to large bio- and nano-systems.

The Physical and Chemical Structure of Hot Molecular Cores

We have made self-consistent models of the density and temperature profiles of the gas and dust surrounding embedded luminous objects using a detailed radiative transfer model together with observations of the spectral energy distribution of hot molecular cores. Using these profiles we have investigated the hot core chemistry which results when grain mantles are evaporated, taking into account the different binding energies of the mantle molecules, as well a model in which we assume that all molecules are embedded in water ice and have a common binding energy. We find that most of the resulting column densities are consistent with those observed toward the hot core G34.3+0.15 at a time around 10^4 years after central luminous star formation. We have also investigated the dependence of the chemical structure on the density profile which suggests an observational possibility of constraining density profiles from determination of the source sizes of line emission from desorbed molecules.

A split-combination approach to merging knowledge bases in possibilistic logic

In this paper, we propose an adaptive approach to merging possibilistic knowledge bases that deploys multiple operators instead of a single operator in the merging process. The merging approach consists of two steps: one is called the splitting step and the other is called the combination step. The splitting step splits each knowledge base into two subbases and then in the second step, different classes of subbases are combined using different operators. Our approach is applied to knowledge bases which are self-consistent and the result of merging is also a consistent knowledge base. Two operators are proposed based on two different splitting methods. Both operators result in a possibilistic knowledge base which contains more information than that obtained by the t-conorm (such as the maximum) based merging methods. In the flat case, one of the operators provides a good alternative to syntax-based merging operators in classical logic.

A stabilization mechanism of zirconia based on oxygen vacancies only

The microscopic mechanism leading to stabilization of cubic and tetragonal forms of zirconia (ZrO2) is analyzed by means of a self-consistent tight-binding model. Using this model, energies and structures of zirconia containing different vacancy concentrations are calculated, equivalent in concentration to the charge compensating vacancies associated with dissolved yttria (Y2O3) in the tetragonal and cubic phase fields (3.2 and 14.4% mol, respectively). The model is shown to predict the large relaxations around an oxygen vacancy, and the clustering of vacancies along the 111 directions, in good agreement with experiments and first principles calculations. The vacancies alone are shown to explain the stabilization of cubic zirconia, and the mechanism is analyzed. (C) 2002 Acta Materialia Inc. Published by Elsevier Science Ltd. All rights reserved.

Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free-energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behavior above the transition temperature.

Development of complex classical force fields through force matching to ab initio data: Application to a room-temperature ionic liquid

Recent experimental neutron diffraction data and ab initio molecular dynamics simulation of the ionic liquid dimethylimidazolium chloride ([dmim]Cl) have provided a structural description of the system at the molecular level. However, partial radial distribution functions calculated from the latter, when compared to previous classical simulation results, highlight some limitations in the structural description offered by force fieldbased simulations. With the availability of ab initio data it is possible to improve the classical description of [dmim]Cl by using the force matching approach, and the strategy for fitting complex force fields in their original functional form is discussed. A self-consistent optimization method for the generation of classical potentials of general functional form is presented and applied, and a force field that better reproduces the observed first principles forces is obtained. When used in simulation, it predicts structural data which reproduces more faithfully that observed in the ab initio studies. Some possible refinements to the technique, its application, and the general suitability of common potential energy functions used within many ionic liquid force fields are discussed.

Non-linear conductance of disordered quantum wires

The self-consistent electron potential in a current-carrying disordered quantum wire is spatially inhomogeneous due to the formation of resistivity dipoles across scattering centres. In this paper it is argued that these inhomogeneities in the potential result in a suppression of the differential conductance of such a wire at finite applied voltage. A semi-classical argument allows this suppression, quadratic in the voltage, to be related directly to the amount of intrinsic defect scattering in the wire. This result is then tested against numerical calculations.

Electromigration of vacancies in copper

The total current-induced force on atoms in a Cu wire containing a vacancy are calculated using the self consistent one-electron density matrix in the presence of an electric current, without separation into electron-wind and direct forces. By integrating the total current-induced force, the change in vacancy migration energy due to the current is calculated. We use the change in migration energy with current to infer an effective electromigration driving force F-e. Finally, we calculate the proportionality constant rho* between F-e and the current density in the wire.

Current-induced forces in atomic-scale conductors

We present a self-consistent tight-binding formalism to calculate the forces on individual atoms due to the flow of electrical current in atomic-scale conductors. Simultaneously with the forces, the method yields the local current density and the local potential in the presence of current flow, allowing a direct comparison between these quantities. The method is applicable to structures of arbitrary atomic geometry and can be used to model current-induced mechanical effects in realistic nanoscale junctions and wires. The formalism is implemented within a simple Is tight-binding model and is applied to two model structures; atomic chains and a nanoscale wire containing a vacancy.

Interactions of 1-Methylimidazole with UO2(CH3CO2)2 and UO2(NO3)2: Structural, Spectroscopic and Theoretical Evidence for Imidazole Binding to the Uranyl Ion

The first definitive high-resolution single-crystal X-ray structure for the coordination of the 1-methylimidazole (Meimid) ligand to UO2(Ac)2 (Ac = CH3CO2) is reported. The crystal structure evidence is confirmed by IR, Raman, and UV-vis spectroscopic data. Direct participation of the nitrogen atom of the Meimid ligand in binding to the uranium center is confirmed. Structural analysis at the DFT (B3LYP) level of theory showed a conformational difference of the Meimid ligand in the free gas-phase complex versus the solid state due to small energetic differences and crystal packing effects. Energetic analysis at the MP2 level in the gas phase supported stronger Meimid binding over H2O binding to both UO2(Ac)2 and UO2(NO3)2. In addition, self-consistent reaction field COSMO calculations were used to assess the aqueous phase energetics of combination and displacement reactions involving H2O and Meimid ligands to UO2R2 (R = Ac, NO3). For both UO2(NO3)2 and UO2(Ac)2, the displacement of H2O by Meimid was predicted to be energetically favorable, consistent with experimental results that suggest Meimid may bind uranyl at physiological pH. Also, log(Knitrate/KAc) calculations supported experimental evidence that the binding stoichiometry of the Meimid ligand is dependent upon the nature of the reactant uranyl complex. These results clearly demonstrate that imidazole binds to uranyl and suggest that binding of histidine residues to uranyl could occur under normal biological conditions.

Magnetic tight binding and the iron-chromium enthalpy anomaly

We describe a self-consistent magnetic tight-binding theory based in an expansion of the Hohenberg-Kohn density functional to second order, about a non-spin-polarized reference density. We show how a first order expansion about a density having a trial input magnetic moment leads to a fixed moment model. We employ a simple set of tight-binding parameters that accurately describes electronic structure and energetics, and show these to be transferable between first row transition metals and their alloys. We make a number of calculations of the electronic structure of dilute Cr impurities in Fe, which we compare with results using the local spin density approximation. The fixed moment model provides a powerful means for interpreting complex magnetic configurations in alloys; using this approach, we are able to advance a simple and readily understood explanation for the observed anomaly in the enthalpy of mixing.

Effect of quantization of vibrations on the structural properties of crystals

We study the structural effects produced by the quantization of vibrational degrees of freedom in periodic crystals at zero temperature. To this end we introduce a methodology based on mapping a suitable subspace of the vibrational manifold and solving the Schrödinger equation in it. A number of increasingly accurate approximations ranging from the quasiharmonic approximation (QHA) to the vibrational self-consistent field (VSCF) method and the exact solution are described. A thorough analysis of the approximations is presented for model monatomic and hydrogen-bonded chains, and results are presented for a linear H-F chain where the potential-energy surface is obtained via first-principles electronic structure calculations. We focus on quantum nuclear effects on the lattice constant and show that the VSCF is an excellent approximation, meaning that correlation between modes is not extremely important. The QHA is excellent for covalently bonded mildly anharmonic systems, but it fails for hydrogen-bonded ones. In the latter, the zero-point energy exhibits a nonanalytic behavior at the lattice constant where the H atoms center, which leads to a spurious secondary minimum in the quantum-corrected energy curve. An inexpensive anharmonic approximation of noninteracting modes appears to produce rather good results for hydrogen-bonded chains for small system sizes. However, it converges to the incorrect QHA results for increasing size. Isotope effects are studied for the first-principles H-F chain. We show how the lattice constant and the H-F distance increase with decreasing mass and how the QHA proves to be insufficient to reproduce this behavior.

Ring currents in azulene

We propose a self consistent polarisable ion tight binding theory for the study of push-pull processes in aromatic molecules. We find that the method quantitatively reproduces ab initio calculations of dipole moments and polarisability. We apply the scheme in a simulation which solves the time dependent Schroedinger equation to follow the relaxation of azulene from the second excited to the ground states. We observe rather spectacular oscillating ring currents which we explain in terms of interference between the HOMO and LUMO states.

Coherent X-ray production via pulse reflection from laser-driven dense electron sheets

A scheme to obtain brilliant x-ray sources by coherent reflection of a counter-propagating pulse from laser-driven dense electron sheets is theoretically and numerically investigated in a self-consistent manner. A radiation pressure acceleration model for the dynamics of the electron sheets blown out from laser-irradiated ultrathin foils is developed and verified by PIC simulations. The first multidimensional and integral demonstration of the scheme by 2D PIC simulations is presented. It is found that the reflected pulse undergoes Doppler-upshift by a factor 4?z2, where ?z = (1- vz2/c2)-1/2 is the effective Lorentz factor of the electron sheet al ong its normal direction. Meanwhile the pulse electric field is intensified by a factor depending on the electron density of the sheet in its moving frame ne/?, where ? is the full Lorentz factor.

Structural changes in quasi-one-dimensional many-electron systems: From linear to zigzag and beyond

Many-electron systems confined to a quasi-one-dimensional geometry by a cylindrical distribution of positive charge have been investigated by density functional computations in the unrestricted local spin density approximation. Our investigations have been focused on the low-density regime, in which electrons are localized. The results reveal a wide variety of different charge and spin configurations, including linear and zig-zag chains, single-and double-strand helices, and twisted chains of dimers. The spin-spin coupling turns from weakly antiferromagnetic at relatively high density, to weakly ferromagnetic at the lowest densities considered in our computations. The stability of linear chains of localized charge has been investigated by analyzing the radial dependence of the self-consistent potential and by computing the dispersion relation of low-energy harmonic excitations.