Ultrafast solvation dynamics in water: Isotope effects and comparison with experimental results

A detailed theoretical study of solvation dynamics in water is presented. The motivation of the present study comes from the recent experimental observation that the dynamics of solvation of an ion in water is ultrafast and the solvation time correlation function decays with a time constant of about 55 fs. The slower decay in the long time can be described by a sum of two exponentials with time constants equal to 126 and 880 fs. The molecular theory (developed earlier) predicts a time constant equal to 52 fs for the initial Gaussian decay and time constants equal to 134 and 886 fs for the two exponential components at the long time. This nearly perfect agreement is obtained by using the most detailed dynamical information available in the literature. The present study emphasizes the importance of the intermolecular vibrational band originating from the O...O stretching mode of the O�H...O units in the initial dynamics and raises several interesting questions regarding the nature of the decay of this mode. We have also studied the effects of isotope substitution on solvation dynamics. It is predicted that a significant isotope effect may be observed in the long time. The experimental results have also been compared with the prediction of the dynamic mean spherical approximation (DMSA); the agreement is not satisfactory at the long time. It is further found that the molecular theory and the DMSA lead to virtually identical results if the translational modes of the solvent molecules are neglected in the former. DMSA has also been used to investigate the dynamics of solvation of a dipolar solute in water. It is found that the dynamics of dipolar solvation exhibit features rather different from those of ion solvation. © 1995 American Institute of Physics.

Instability, intermittency, and multiscaling in discrete growth models of kinetic roughening

We show by numerical simulations that discretized versions of commonly studied continuum nonlinear growth equations (such as the Kardar-Parisi-Zhangequation and the Lai-Das Sarma-Villain equation) and related atomistic models of epitaxial growth have a generic instability in which isolated pillars (or grooves) on an otherwise flat interface grow in time when their height (or depth) exceeds a critical value. Depending on the details of the model, the instability found in the discretized version may or may not be present in the truly continuum growth equation, indicating that the behavior of discretized nonlinear growth equations may be very different from that of their continuum counterparts. This instability can be controlled either by the introduction of higher-order nonlinear terms with appropriate coefficients or by restricting the growth of pillars (or grooves) by other means. A number of such ''controlled instability'' models are studied by simulation. For appropriate choice of the parameters used for controlling the instability, these models exhibit intermittent behavior, characterized by multiexponent scaling of height fluctuations, over the time interval during which the instability is active. The behavior found in this regime is very similar to the ''turbulent'' behavior observed in recent simulations of several one- and two-dimensional atomistic models of epitaxial growth.

Nonequilibrium phase transition in the kinetic Ising model: Critical slowing down and the specific-heat singularity

The nonequilibrium dynamic phase transition, in the kinetic Ising model in the presence of an oscillating magnetic field has been studied both by Monte Carlo simulation and by solving numerically the mean-field dynamic equation of motion for the average magnetization. In both cases, the Debye ''relaxation'' behavior of the dynamic order parameter has been observed and the ''relaxation time'' is found to diverge near the dynamic transition point. The Debye relaxation of the dynamic order parameter and the power law divergence of the relaxation time have been obtained from a very approximate solution of the mean-field dynamic equation. The temperature variation of appropriately defined ''specific heat'' is studied by the Monte Carlo simulation near the transition point. The specific heat has been observed to diverge near the dynamic transition point.

Relative performances of textural models in estimating soil moisture characteristic

The soil moisture characteristic (SMC) forms an important input to mathematical models of water and solute transport in the unsaturated-soil zone. Owing to their simplicity and ease of use, texture-based regression models are commonly used to estimate the SMC from basic soil properties. In this study, the performances of six such regression models were evaluated on three soils. Moisture characteristics generated by the regression models were statistically compared with the characteristics developed independently from laboratory and in-situ retention data of the soil profiles. Results of the statistical performance evaluation, while providing useful information on the errors involved in estimating the SMC, also highlighted the importance of the nature of the data set underlying the regression models. Among the models evaluated, the one possessing an underlying data set of in-situ measurements was found to be the best estimator of the in-situ SMC for all the soils. Considerable errors arose when a textural model based on laboratory data was used to estimate the field retention characteristics of unsaturated soils.

Screw dynamo and the generation of nonaxisymmetric magnetic fields

A mechanism is presented here for the amplification of large-scale nonaxisymmetric magnetic fields as a manifestation of the dynamo effect. We generalize a result on restrictions of dynamo actions due to laminar flow originally derived by Zeldovich, Ruzmaikin, and Sokolov [Magnetic Fields in Astrophysics (Gordon and Breach, New York, 1983)]. We show how a screwlike motion having phi and z components of velocity can help to grow a magnetic field. This model postulates a large-scale flow having phi and z components with radial dependences (helical flow). Shear in the radial field, because of a near-flux-freezing condition, causes amplification of the phi component of the magnetic field. The radial and axial components grow due to the presence of turbulent diffusion. The shear in the large scale flow induces an indefinite growth of magnetic field without the a effect; nevertheless, turbulent diffusion forms an important part in the overall mechanism.

Structural sensitivity as a measure of redundancy

The conventional definition of redundancy is applicable to skeletal structural systems only, whereas the concept of redundancy has never been discussed in the context of a continuum. Generally, structures in civil engineering constitute a combination of both skeletal and continuum segments. Hence, this gaper presents a generalized definition of redundancy that has been defined in terms of structural response sensitivity, which is applicable to both continuum and discrete structures. In contrast to the conventional definition of redundancy, which is assumed to be fixed for a given structure and is believed to be independent of loading and material properties, the new definition would depend on strength and response of the structure at a given stage of its service life. The redundancy measure proposed in this paper is linked to the structural response sensitivities. Thus, the structure can have different degrees of redundancy during its lifetime, depending on the response sensitivity under consideration It is believed that this new redundancy measure would be more relevant in structural evaluation, damage assessment, and reliability analysis of structures at large.

Metal-insulator transition in a degenerate Hubbard model

We have studied the metal-insulator transition at integer fillings in a triply degenerate Hubbard model using the Lanczos method. The critical Coulomb interaction strength U-c, is found to depend strongly on the band filling, with U-c similar to root 3 W (W is the bandwidth) at half filling for this case with threefold degeneracy. We discuss the implications of our results on metal-insulator transitions in strongly correlated systems in general, and on the unusual electronic ground state of the alkali-metal-doped fullerenes, in particular. [S0163-1829(99)11003-8].

Micellar structures of dimeric surfactants with phosphate head groups and wettable spacers: A small-angle neutron scattering study

Dimeric or gemini surfactants consist of two hydrophobic chains and two hydrophilic head groups covalently connected by a hydrophobic or hydrophilic spacer. This paper reports the small-angle neutron scattering (SANS) measurements from aqueous micellar solutions of two different recently developed types of dimeric surfactants: (i) bis-anionic C16H33PO4--(CH2)(m)-PO4-C16H33,2Na(+) dimeric surfactants composed of phosphate head groups and a hydrophobic polymethylene spacer, referred to as 16-m-16,2Na(+), for spacer lengths m = 2, 4, 6, and 10, (ii) bis-cationic C16H33N+(CH3)(2)-CH2-(CH2-O-CH2)(p)-CH2-N+ (CH3)(2)C16H33,2Br(-) dimeric surfactants composed of dimethylammonium head groups and a wettable polyethylene oxide spacer, referred to as 16-CH2-p-CH2-16,2Br(-), for spacer lengths p = 1 - 3. The micellar structures of these surfactants are compared with the earlier studied bis-cationic C16H33N+ (CH3)(2)-(CH2)(m)-N+ (CH3)(2)C16H33,2Br(-) dimeric surfactants composed of dimethylammonium head groups and a hydrophobic polymethylene spacer, referred to as 16-m-16,2Br(-). It is found that 16-m-16,2Na(+), similar to 16-m-16,2Br(-), form various micellar structures depending on the spacer length. Micelles an disklike for rn = 2, rodlike for m = 4, and prolate ellipsoidal fur m = 6 and 10. The micelles of 16-CH2-p-CH2-16,2Br(-) are prolate ellipsoidal for all the values of p = 1 - 3. It is also found that micelles of 16-m-16,2Na(+) and 16-CH2-p-CH2-16,2Br(-) are large in comparison to those of 16-in-16,2Br(-) for similar spacer lengths. This is connected with the fact that both in 16-m-16,2Na(+) and 16-CH2-p-CH2-16,2Br(-), the head group or the spacer is more hydrated as compared to that in the 16-m-16,2Br(-). An increase in the hydration of the spacer or the head group increases the screening of the Coulomb repulsion between the charged head groups. This effect has been found to be more pronounced in the dimeric surfactants having wettable spacers. [S1063-651X(99)00303-7].

Prediction of reaeration rates in square, stirred tanks

Experiments were conducted on the oxygen transfer coefficient, k(L)a(20), through surface aeration in geometrically similar square tanks, with a rotor of diameter D fitted with six flat blades. An optimal geometric similarity of various linear dimensions, which produced maximum k(L)a(20) for any rotational speed of rotor N by an earlier study, was maintained. A simulation equation uniquely correlating k = k(L)a(20)(nu/g(2))(1/3) (nu and g are kinematic viscosity of water and gravitational constant, respectively), and a parameter governing the theoretical power per unit volume, X = (ND2)-D-3/(g(4/3)nu(1/3)), is developed. Such a simulation equation can be used to predict maximum k for any N in any size of such geometrically similar square tanks. An example illustrating the application of results is presented. Also, it has been established that neither the Reynolds criterion nor the Froude criterion is singularly valid to simulate either k or K = k(L)a(20)/N, simultaneously in all the sizes of tanks, even through they are geometrically similar. Occurrence of "scale effects" due to the Reynolds and the Froude laws of similitude on both k and K are also evaluated.

Density-matrix renormalization-group study of low-lying excitations of polyacene within a Pariser-Parr-Pople model

We have carried out symmetrized density-matrix renormalization-group calculations to study the nature of excited states of long polyacene oligomers within a Pariser-Parr-Pople Hamiltonian. We have used the C-2 symmetry, the electron-hole symmetry, and the spin parity of the system in our calculations. We find that there is a crossover in the lowest dipole forbidden two-photon state and the lowest dipole allowed excited state with size of the oligomer. In the long system limit, the two-photon state lies below the lowest dipole allowed excited state. The triplet state lies well below the two-photon state and energetically does not correspond to its description as being made up of two triplets. These results are in agreement with the general trends in linear conjugated polymers. However, unlike in linear polyenes wherein the two-photon state is a localized excitation, we find that in polyacenes, the two-photon excitation is spread out over the system. We have doped the systems with a hole and an electron and have calculated the charge excitation gap. Using the charge gap and the optical gap, we estimate the binding energy of the 1(1)B(-) exciton to be 2.09 eV. We have also studied doubly doped polyacenes and find that the bipolaron in these systems, to be composed of two separated polarons, as indicated by the calculated charge-density profile and charge-charge correlation function. We have studied bond orders in various states in order to get an idea of the excited state geometry of the system. We find that the ground state, the triplet state, the dipole allowed state, and the polaron excitations correspond to lengthening of the rung bonds in the interior of the oligomer while the two-photon excitation corresponds to the rung bond lengths having two maxima in the system.

Effect of two-level systems on the spectral density of universal conductance fluctuations in doped Si

We have studied the power spectral density [S(f) = gamma/f(alpha)] of universal conductance fluctuations (UCF's) in heavily doped single crystals of Si, when the scatterers themselves act as the primary source of dephasing. We observed that the scatterers, with internal dynamics like two-level-systems, produce a significant, temperature-dependent reduction in the spectral slope alpha when T less than or similar to 10 K, as compared to the bare 1/f (alphaapproximate to1) spectrum at higher temperatures. It is further shown that an upper cutoff frequency (f(m)) in the UCF spectrum is necessary in order to restrict the magnitude of conductance fluctuations, [(deltaG(phi))(2)], per phase coherent region (L-phi(3)) to [(deltaGphi)(2)](1/2) less than or similar to e(2)/h. We find that f(m) approximate to tau(D)(-1), where tau(D) = L-2/D, is the time scale of the diffusive motion of the electron along the active length (L) of the sample (D is the electron diffusivity).

K(2)O(2): The most stable oxide of K

We have analyzed the stability of various oxides of K and find that K(2)O(2) is the most stable one. The additional stability is traced to the presence of oxygen dimers in K(2)O(2) which interact to form molecular orbitals. Other oxides such as KO(2) and KO(3) which also have dimers/trimers of oxygens are found to be less stable. This is traced to the shorter O-O bonds that one finds in them which gives rise to a significant coulomb repulsion between the electrons on the oxygen atoms making up the dimer/trimer, making them less stable.

Two-channel Kondo physics in odd impurity chains

We study odd-membered chains of spin-1/2 impurities, with each end connected to its own metallic lead. For antiferromagnetic exchange coupling, universal two-channel Kondo (2CK) physics is shown to arise at low energies. Two overscreening mechanisms are found to occur depending on coupling strength, with distinct signatures in physical properties. For strong interimpurity coupling, a residual chain spin-1/2 moment experiences a renormalized effective coupling to the leads, while in the weak-coupling regime, Kondo coupling is mediated via incipient single-channel Kondo singlet formation. We also investigate models in which the leads are tunnel-coupled to the impurity chain, permitting variable dot filling under applied gate voltages. Effective low-energy models for each regime of filling are derived, and for even fillings where the chain ground state is a spin singlet, an orbital 2CK effect is found to be operative. Provided mirror symmetry is preserved, 2CK physics is shown to be wholly robust to variable dot filling; in particular, the single-particle spectrum at the Fermi level, and hence the low-temperature zero-bias conductance, is always pinned to half-unitarity. We derive a Friedel-Luttinger sum rule and from it show that, in contrast to a Fermi liquid, the Luttinger integral is nonzero and determined solely by the ``excess'' dot charge as controlled by gate voltage. The relevance of the work to real quantum dot devices, where interlead charge-transfer processes fatal to 2CK physics are present, is also discussed. Physical arguments and numerical renormalization-group techniques are used to obtain a detailed understanding of these problems.

Model of resistances in systems of Tomonaga-Luttinger liquid wires

In a recent paper, we combined the technique of bosonization with the concept of a Rayleigh dissipation function to develop a model for resistances in one-dimensional systems of interacting spinless electrons Europhys. Lett. 93, 57007 (2011)]. We also studied the conductance of a system of three wires by using a current splitting matrix M at the junction. In this paper, we extend our earlier work in several ways. The power dissipated in a three-wire system is calculated as a function of M and the voltages applied in the leads. By combining two junctions of three wires, we examine a system consisting of two parallel resistances. We study the conductance of this system as a function of the M matrices and the two resistances; we find that the total resistance is generally quite different from what one expects for a classical system of parallel resistances. We do a sum over paths to compute the conductance of this system when one of the two resistances is taken to be infinitely large. We study the conductance of a three-wire system of interacting spin-1/2 electrons, and show that the charge and spin conductances can generally be different from each other. Finally, we consider a system of two wires that are coupled by a dissipation function, and we show that this leads to a current in one wire when a voltage bias is applied across the other wire.

Modeling reservoir irrigation in uncertain hydrologic environment

In a detailed model for reservoir irrigation taking into account the soil moisture dynamics in the root zone of the crops, the data set for reservoir inflow and rainfall in the command will usually be of sufficient length to enable their variations to be described by probability distributions. However, the potential evapotranspiration of the crop itself depends on the characteristics of the crop and the reference evaporation, the quantification of both being associated with a high degree of uncertainty. The main purpose of this paper is to propose a mathematical programming model to determine the annual relative yield of crops and to determine its reliability, for a single reservoir meant for irrigation of multiple crops, incorporating variations in inflow, rainfall in the command area, and crop consumptive use. The inflow to the reservoir and rainfall in the reservoir command area are treated as random variables, whereas potential evapotranspiration is modeled as a fuzzy set. The model's application is illustrated with reference to an existing single-reservoir system in Southern India.