Exchange-correlation energy and potential as approximate functionals of occupied and virtual Kohn-Sham orbitals: Application to dissociating H-2

The standard local density approximation and generalized gradient approximations fail to properly describe the dissociation of an electron pair bond, yielding large errors (on the order of 50 kcal/mol) at long bond distances. To remedy this failure, a self-consistent Kohn-Sham (KS) method is proposed with the exchange-correlation (xc) energy and potential depending on both occupied and virtual KS orbitals. The xc energy functional of Buijse and Baerends [Mol. Phys. 100, 401 (2002); Phys. Rev. Lett. 87, 133004 (2001)] is employed, which, based on an ansatz for the xc-hole amplitude, is able to reproduce the important dynamical and nondynamical effects of Coulomb correlation through the efficient use of virtual orbitals. Self-consistent calculations require the corresponding xc potential to be obtained, to which end the optimized effective potential (OEP) method is used within the common energy denominator approximation for the static orbital Green's function. The problem of the asymptotic divergence of the xc potential of the OEP when a finite number of virtual orbitals is used is addressed. The self-consistent calculations reproduce very well the entire H-2 potential curve, describing correctly the gradual buildup of strong left-right correlation in stretched H-2. (C) 2003 American Institute of Physics.

A theoretical color-velocity correlation for supernovae associated with gamma-ray bursts

We carry out the first multi-dimensional radiative transfer calculations to simultaneously compute synthetic spectra and light curves for models of supernovae driven by fast bipolar outflows. These allow us to make self-consistent predictions for the orientation dependence of both color evolution and spectral features. We compare models with different degrees of asphericity and metallicity and find significant observable consequences of both. In aspherical models, we find spectral and light curve features that vary systematically with observer orientation. In particular, we find that the early-phase light curves are brighter and bluer when viewed close to the polar axis but that the peak flux is highest for equatorial (off-axis) inclinations. Spectral line features also depend systematically on observer orientation, including the velocity of the Si II 6355 Å line. Consequently, our models predict a correlation between line velocity and color that could assist the identification of supernovae associated with off-axis jet-driven explosions. The amplitude and range of this correlation depends on the degree of asphericity, the metallicity, and the epoch of observation but we find that it is always present and acts in the same direction. © 2012. The American Astronomical Society. All rights reserved..

Universal behaviour of shock precursors in the presence of efficient cosmic ray acceleration

The self-consistent interaction between energetic particles and self-generated hydromagnetic waves in a cosmic ray pressure dominated plasma is considered. Using a three-dimensional hybrid magnetohydrodynamics (MHD)-kinetic code, which utilizes a spherical harmonic expansion of the Vlasov-Fokker-Planck equation, high-resolution simulations of the magnetic field growth including feedback on the cosmic rays are carried out. It is found that for shocks with high cosmic ray acceleration efficiency, the magnetic fields become highly disorganized, resulting in near isotropic diffusion, independent of the initial orientation of the ambient magnetic field. The possibility of sub-Bohm diffusion is demonstrated for parallel shocks, while the diffusion coefficient approaches the Bohm limit from below for oblique shocks. This universal behaviour suggests that Bohm diffusion in the root-mean-squared field inferred from observation may provide a realistic estimate for the maximum energy acceleration time-scale in young supernova remnants. Although disordered, the magnetic field is not self-similar suggesting a non-uniform energy-dependent behaviour of the energetic particle transport in the precursor. Possible indirect radiative signatures of cosmic ray driven magnetic field amplification are discussed.

Quantum chemistry and in situ FTir spectroscopy studies on potential-dependent properties of CO adsorbed on Pt electrodes

The electronic and vibrational properties of CO adsorbed on Pt electrodes at different potentials have been studied, by using methods of self-consistent-charge discrete variational Xa (SCC-DV-Xa) cluster calculations and in situ FTir spectroscopy. Two new models have been developed and verified to be successful: (1) using a "metallic state cluster" to imitate a metal (electrode) surface; and (2) charging the cluster and shifting its Fermi level (e{lunate}) to simulate, according to the relation of -d e{lunate}e dE, quantitatively the variation of the electrode potential (E). It is shown that the binding of PtCO is dominated by the electric charge transfer of dp ? 2p, while that of s ? Pt is less important in this binding. The electron occupancy of the 2p orbital of CO weakens the CO bond and decreases the v. Variation of E mainly influences the charge transfer process of dp ? 2p, but hardly influences that of s ? Pt. A linear potential-dependence of v has been shown and the calculated dv/dE = 35.0 cm V. All results of calculations coincide with the ir experimental data. © 1993.

Quantum chemistry studies on electronic properties of CN adsorbed on silver electrodes

The electronic properties of CN adsorbed on Ag electrodes at different potentials have been studied by using the method of self-consistent-charge discrete variational Xa (SCC-DV-Xa) cluster calculations. It is shown that the binding of NCAg is dominated by both electrostatic and polarization effects derived from the charge of CN, and that the transfer of s charge from CN to silver is the largest donation contribution. The electrode potential influences this charge transfer process and the interaction between CN adsorbate and silver electrode. The results of quantum chemistry calculations fit well with the experimental data of in situ spectroscopic studies on the CN/Ag electrode systems. © 1991.

Effect of magnetic bias on Gaussian pulse scattering by stacked nonlinear semiconductor layers

The nonlinear scattering of pulses by periodic stacks of semiconductor layers with magnetic bias has been studied in the self-consistent problem formulation, taking into account mobility of carriers. The three-wave mixing technique has been applied to the analysis of the waveform evolution in the stacks illuminated by two Gaussian pulses with different central frequencies and lengths. The effects of external magnetic bias, and stack physical and geometrical parameters on the properties of the scattered waveforms are discussed. © 2013 IEEE.

Pulsed second harmonic generation by stacks of magnetically biased semiconductor layers

The pulsed second harmonic generation (SHG) by periodic stacks of nonlinear semiconductor layers with external magnetic bias has been studied in the self-consistent problem formulation, taking into account mobility of carriers. The products of nonlinear scattering in the three-wave mixing process are examined. It is demonstrated that the waveform evolution in magnetoactive weakly nonlinear semiconductor periodic structure illuminated by Gaussian pulse is strongly affected by the magnetic bias and collision frequency of the carriers. The effect of nonreciprocity on the SHG efficiency is discussed and illustrated by the examples. © 2013 European Microwave Association.

Line-driven Disk Winds in Activ Galactic Nuclei: The Critical Importance of Ionization and Radiative Transfer

Accretion disk winds are thought to produce many of the characteristic features seen in the spectra of active galactic nuclei (AGNs) and quasi-stellar objects (QSOs). These outflows also represent a natural form of feedback between the central supermassive black hole and its host galaxy. The mechanism for driving this mass loss remains unknown, although radiation pressure mediated by spectral lines is a leading candidate. Here, we calculate the ionization state of, and emergent spectra for, the hydrodynamic simulation of a line-driven disk wind previously presented by Proga & Kallman. To achieve this, we carry out a comprehensive Monte Carlo simulation of the radiative transfer through, and energy exchange within, the predicted outflow. We find that the wind is much more ionized than originally estimated. This is in part because it is much more difficult to shield any wind regions effectively when the outflow itself is allowed to reprocess and redirect ionizing photons. As a result, the calculated spectrum that would be observed from this particular outflow solution would not contain the ultraviolet spectral lines that are observed in many AGN/QSOs. Furthermore, the wind is so highly ionized that line driving would not actually be efficient. This does not necessarily mean that line-driven winds are not viable. However, our work does illustrate that in order to arrive at a self-consistent model of line-driven disk winds in AGN/QSO, it will be critical to include a more detailed treatment of radiative transfer and ionization in the next generation of hydrodynamic simulations.

Universal tight binding model for chemical reactions in solution and at surfaces. II. Water

A revised water model intended for use in condensed phase simulations in the framework of the self consistent polarizable ion tight binding theory is constructed. The model is applied to water monomer, dimer, hexamers, ice, and liquid, where it demonstrates good agreement with theoretical results obtained by more accurate methods, such as DFT and CCSD(T), and with experiment. In particular, the temperature dependence of the self diffusion coefficient in liquid water predicted by the model, closely reproduces experimental curves in the temperature interval between 230 K and 350 K. In addition, and in contrast to standard DFT, the model properly orders the relative densities of liquid water and ice. A notable, but inevitable, shortcoming of the model is underestimation of the static dielectric constant by a factor of two. We demonstrate that the description of inter and intramolecular forces embodied in the tight binding approximation in quantum mechanics leads to a number of valuable insights which can be missing from ab initio quantum chemistry and classical force fields. These include a discussion of the origin of the enhanced molecular electric dipole moment in the condensed phases, and a detailed explanation for the increase of coordination number in liquid water as a function of temperature and compared with ice-leading to insights into the anomalous expansion on freezing. The theory holds out the prospect of an understanding of the currently unexplained density maximum of water near the freezing point. 

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules

As is now well established, a first order expansion of the Hohenberg-Kohn total energy density functional about a trial input density, namely, the Harris-Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis.

Nonreciprocal scattering by stacked nonlinear magneto-active semiconductor layers

The combinatorial frequency generation by the periodic stacks of magnetically biased semiconductor layers has been modelled in a self-consistent problem formulation, taking into account the nonlinear dynamics of carriers. It is shown that magnetic bias not only renders nonreciprocity of the three-wave mixing process but also significantly enhances the nonlinear interactions in the stacks, especially at the frequencies close to the intrinsic magneto-plasma resonances of the constituent layers. The main mechanisms and properties of the combinatorial frequency generation and emission from the stacks are illustrated by the simulation results, and the effects of the individual layer parameters and the structure arrangement on the stack nonlinear and nonreciprocal response are discussed. © 2014 Elsevier B.V. All rights reserved.

Gaussian pulse mixing by stacked semiconductor layers

The nonlinear scattering of two Gaussian pulses with different central frequencies incident at slant angles on the periodic stack of binary semiconductor layers has been modelled in the self-consistent problem formulation taking into account the dynamics of charges. The effects of the pump pulse length and central frequencies, and the stack physical and geometrical parameters on the properties of the emitted combinatorial frequency waveforms are analysed and discussed.

Pulse mixing by Fibonacci and Thue-Morse ordered semiconductor layers

The properties of combinatorial frequency generation by two-tone Gaussian pulses incident at oblique angles on quasiperiodic (Fibonacci and Thue-Morse) stacks of binary semiconductor layers are discussed. The analysis has been performed using the self-consistent model taking into account the nonlinear dynamics of mobile charges in the layers. The effects of the stack arrangements and constituent layer parameters on the combinatorial frequency waveforms are presented for the specific structures of both types

Modeling relativistic soliton interactions in overdense plasmas: A perturbed nonlinear Schrodinger equation framework

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrodinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrodinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude

Relativistic theory for localized electrostatic excitations in degenerate electron-ion plasmas

A self-consistent relativistic two-fluid model is proposed for electron-ion plasma dynamics. A one-dimensional geometry is adopted. Electrons are treated as a relativistically degenerate fluid, governed by an appropriate equation of state. The ion fluid is also allowed to be relativistic, but is cold, nondegenerate, and subject only to an electrostatic potential. Exact stationary-profile solutions are sought, at the ionic scale, via the Sagdeev pseudopotential method. The analysis provides the pulse existence region, in terms of characteristic relativistic parameters, associated with the (ultrahigh) particle density.