Ab initio absorption spectrum of NiO combining molecular dynamics with the embedded cluster approach in a discrete reaction field

We developed a procedure that combines three complementary computational methodologies to improve the theoretical description of the electronic structure of nickel oxide. The starting point is a Car-Parrinello molecular dynamics simulation to incorporate vibrorotational degrees of freedom into the material model. By means ofcomplete active space self-consistent field second-order perturbation theory (CASPT2) calculations on embedded clusters extracted from the resulting trajectory, we describe localized spectroscopic phenomena on NiO with an efficient treatment of electron correlation. The inclusion of thermal motion into the theoretical description allowsus to study electronic transitions that, otherwise, would be dipole forbidden in the ideal structure and results in a natural reproduction of the band broadening. Moreover, we improved the embedded cluster model by incorporating self-consistently at the complete active space self-consistent field (CASSCF) level a discrete (or direct) reaction field (DRF) in the cluster surroundings. The DRF approach offers an efficient treatment ofelectric response effects of the crystalline embedding to the electronic transitions localized in the cluster. We offer accurate theoretical estimates of the absorption spectrum and the density of states around the Fermi level of NiO, and a comprehensive explanation of the source of the broadening and the relaxation of the charge transferstates due to the adaptation of the environment

A variational approach for calculating Franck-Condon factors including mode-mode anharmonic coupling

We have implemented our new procedure for computing Franck-Condon factors utilizing vibrational configuration interaction based on a vibrational self-consistent field reference. Both Duschinsky rotations and anharmonic three-mode coupling are taken into account. Simulations of the first ionization band of Cl O2 and C4 H4 O (furan) using up to quadruple excitations in treating anharmonicity are reported and analyzed. A developer version of the MIDASCPP code was employed to obtain the required anharmonic vibrational integrals and transition frequencies

Fronts from complex two-dimensional dispersal kernels: theory and application to Reid's paradox

Bimodal dispersal probability distributions with characteristic distances differing by several orders of magnitude have been derived and favorably compared to observations by Nathan [Nature (London) 418, 409 (2002)]. For such bimodal kernels, we show that two-dimensional molecular dynamics computer simulations are unable to yield accurate front speeds. Analytically, the usual continuous-space random walks (CSRWs) are applied to two dimensions. We also introduce discrete-space random walks and use them to check the CSRW results (because of the inefficiency of the numerical simulations). The physical results reported are shown to predict front speeds high enough to possibly explain Reid's paradox of rapid tree migration. We also show that, for a time-ordered evolution equation, fronts are always slower in two dimensions than in one dimension and that this difference is important both for unimodal and for bimodal kernels

Relaxation times of non-Markovian processes

We consider a general class of non-Markovian processes defined by stochastic differential equations with Ornstein-Uhlenbeck noise. We present a general formalism to evaluate relaxation times associated with correlation functions in the steady state. This formalism is a generalization of a previous approach for Markovian processes. The theoretical results are shown to be in satisfactory agreement both with experimental data for a cubic bistable system and also with a computer simulation of the Stratonovich model. We comment on the dynamical role of the non-Markovianicity in different situations.

Numerical algorithm for Ginzburg-Landau equations with multiplicative noise: Application to domain growth

We consider stochastic partial differential equations with multiplicative noise. We derive an algorithm for the computer simulation of these equations. The algorithm is applied to study domain growth of a model with a conserved order parameter. The numerical results corroborate previous analytical predictions obtained by linear analysis.

Spatiotemporal order out of noise

Experimental observations of self-organized behavior arising out of noise are also described, and details on the numerical algorithms needed in the computer simulation of these problems are given.

Analysis of the role of mobility-lifetime products in the performance of amorphous silicon p-i-n solar cells

An analytical model of an amorphous silicon p-i-n solar cell is presented to describe its photovoltaic behavior under short-circuit conditions. It has been developed from the analysis of numerical simulation results. These results reproduce the experimental illumination dependence of short-circuit resistance, which is the reciprocal slope of the I(V) curve at the short-circuit point. The recombination rate profiles show that recombination in the regions of charged defects near the p-i and i-n interfaces should not be overlooked. Based on the interpretation of the numerical solutions, we deduce analytical expressions for the recombination current and short-circuit resistance. These expressions are given as a function of an effective ¿¿ product, which depends on the intensity of illumination. We also study the effect of surface recombination with simple expressions that describe its influence on current loss and short-circuit resistance.

Algebraic decay of velocity fluctuations in a confined fluid

Computer simulations of a colloidal particle suspended in a fluid confined by rigid walls show that, at long times, the velocity correlation function decays with a negative algebraic tail. The exponent depends on the confining geometry, rather than the spatial dimensionality. We can account for the tail by using a simple mode-coupling theory which exploits the fact that the sound wave generated by a moving particle becomes diffusive.

Random mixtures with orientational order, and the anisotropic resistivity tensor of high Tc superconductors

By generalizing effective-medium theory to the case of orientationally ordered but positionally disordered two component mixtures, it is shown that the anisotropic dielectric tensor of oxide superconductors can be extracted from microwave measurements on oriented crystallites of YBa2Cu3O7¿x embedded in epoxy. Surprisingly, this technique appears to be the only one which can access the resistivity perpendicular to the copper¿oxide planes in crystallites that are too small for depositing electrodes. This possibility arises in part because the real part of the dielectric constant of oxide superconductors has a large magnitude. The validity of the effective-medium approach for orientationally ordered mixtures is corroborated by simulations on two¿dimensional anisotropic random resistor networks. Analysis of the experimental data suggests that the zero-temperature limit of the finite frequency resistivity does not vanish along the c axis, a result which would simply the existence of states at the Fermi surface, even in the superconducting state

Finite-size scaling investigation of the liquid-liquid critical point in ST2 water and its stability with respect to crystallization.

The liquid-liquid critical point scenario of water hypothesizes the existence of two metastable liq- uid phases low-density liquid (LDL) and high-density liquid (HDL) deep within the supercooled region. The hypothesis originates from computer simulations of the ST2 water model, but the stabil- ity of the LDL phase with respect to the crystal is still being debated. We simulate supercooled ST2 water at constant pressure, constant temperature, and constant number of molecules N for N ≤ 729 and times up to 1 μs. We observe clear differences between the two liquids, both structural and dynamical. Using several methods, including finite-size scaling, we confirm the presence of a liquid-liquid phase transition ending in a critical point. We find that the LDL is stable with respect to the crystal in 98% of our runs (we perform 372 runs for LDL or LDL-like states), and in 100% of our runs for the two largest system sizes (N = 512 and 729, for which we perform 136 runs for LDL or LDL-like states). In all these runs, tiny crystallites grow and then melt within 1 μs. Only for N ≤ 343 we observe six events (over 236 runs for LDL or LDL-like states) of spontaneous crystal- lization after crystallites reach an estimated critical size of about 70 ± 10 molecules.

Morphology and segregation in two-component diffusion-limited aggregation

A diffusion-limited-aggregation (DLA) model with two components (A and B species) is presented to investigate the structure of the composite deposits. The sticking probability PAB (=PBA) between the different species is introduced into the original DLA model. By using computer simulation it is shown that various patterns are produced with varying the sticking probabilities PAB (=PBA) and PAA (= PBB), where PAA (=PBB) is the sticking probability between the same species. Segregated patterns can be analyzed under the condition PAB < PAA, assumed throughout the paper. With decreasing sticking probability PAB, a clustering of the same species occurs. With sufficiently small values of both sticking probabilities PAB and PAA, the deposit becomes dense and the segregated patterns of the composite deposit show a striped structure. The effect of the concentration on the pattern morphology is also shown.

Spatiotemporal order out of noise

Experimental observations of self-organized behavior arising out of noise are also described, and details on the numerical algorithms needed in the computer simulation of these problems are given.

Spatial behavior in groups: an agent-based approach

We present an agent-based model with the aim of studying how macro-level dynamics of spatial distances among interacting individuals in a closed space emerge from micro-level dyadic and local interactions. Our agents moved on a lattice (referred to as a room) using a model implemented in a computer program called P-Space in order to minimize their dissatisfaction, defined as a function of the discrepancy between the real distance and the ideal, or desired, distance between agents. Ideal distances evolved in accordance with the agent's personal and social space, which changed throughout the dynamics of the interactions among the agents. In the first set of simulations we studied the effects of the parameters of the function that generated ideal distances, and in a second set we explored how group macrolevel behavior depended on model parameters and other variables. We learned that certain parameter values yielded consistent patterns in the agents' personal and social spaces, which in turn led to avoidance and approaching behaviors in the agents. We also found that the spatial behavior of the group of agents as a whole was influenced by the values of the model parameters, as well as by other variables such as the number of agents. Our work demonstrates that the bottom-up approach is a useful way of explaining macro-level spatial behavior. The proposed model is also shown to be a powerful tool for simulating the spatial behavior of groups of interacting individuals.

Short-time dynamics of colloidal suspensions in confined geometries

We analyze the short-time dynamical behavior of a colloidal suspension in a confined geometry. We analyze the relevant dynamical response of the solvent, and derive the temporal behavior of the velocity autocorrelation function, which exhibits an asymptotic negative algebraic decay. We are able to compare quantitatively with theoretical expressions, and analyze the effects of confinement on the diffusive behavior of the suspension.

Effect of drift on segregation in two-component diffusion-limited aggregation

An effect of drift is investigated on the segregation pattern in diffusion-limited aggregation (DLA) with two components (A and B species). The sticking probability PAB (=PBA) between the different species is introduced into the DLA model with drift, where the sticking probability PAA (=PBB) between the same species equals 1. By using computer simulation it is found that the drift has an important effect on not only the morphology but also the segregation pattern. Under the drift and the small sticking probability, a characteristic pattern appears where elongated clusters of A species and of B species are periodically dispersed. The period decreases with increasing drift. The periodic structure of the deposits is characterized by an autocorrelation function. The shape of the cluster consisting of only A species (or B species) shows a vertically elongated filamentlike structure. Each cluster becomes vertically longer with decreasing sticking probability PAB. The segregation pattern is distinctly different from that with no drift and a small sticking probability PAA. The effect of the concentration on the segregation pattern is also shown.