Simplified theory of wideband spectra of liquid H2O and D2O (from 0 to 1000 cm(-1)) due to reorienting polar and vibrating H-bonded water molecules

A semi-phenomenological molecular model is presented, which is capable of describing with the use of analytical formulae, the wideband dielectric(1) and far-infrared spectra of ordinary and heavy water. In the model the vector of a dipole moment is presented as a sum of two components. The absolute value of the first one is constant; the second one changes harmonically with time. The key aspect of this work is consideration of FIR spectra due to the second component. In the context of the modified hybrid model presented in the work, reorientation of the dipoles in the rectangular potential well is considered, as a result of which the librational (near 700 cm (-1)) and translational (near 200 cm (-1)) absorption bands and the microwave Debye relaxation spectrum arise. It is shown that the time-dependent part of a dipole moment contributes most to the translational band, the relevant mechanism is taken to be stretching vibration of the H-bonded molecules. Previous linear-response molecular models were unsuccessful in describing this band (in heavy water) in terms of the complex dielectric permittivity. The spatial and time scales characteristic of water are estimated. (C) 2002 Elsevier Science B.V. All rights reserved.

Analytical calculation of far infrared spectra of ice in terms of a molecular model

Wideband far infrared (FIR) spectra of complex permittivity e(p) of ice are calculated in terms of a simple analytical theory based on the method of dipolar autocorrelation functions. The molecular model represents a revision of the model recently presented for liquid water in Adv. Chem. Phys. 127 (2003) 65. A composite two-fractional model is proposed. The model is characterised by three phenomenological potential wells corresponding to the three FIR bands observed in ice. The first fraction comprises dipoles reorienting in a rather narrow and deep hat-like well; these dipoles generate the librational band centred at the frequency approximate to 880 cm(-1). The second fraction comprises elastically interacting particles; they generate two nearby bands placed around frequency 200 cm(-1). For description of one of these bands the harmonic oscillator (HO) model is used, in which translational oscillations of two charged molecules along the H-bond are considered. The other band is produced by the H-bond stretch, which governs hindered rotation of a rigid dipole. Such a motion and its dielectric response are described in terms of a new cut parabolic (CP) model applicable for any vibration amplitude. The composite hat-HO-CP model results in a smooth epsilon(nu) ice spectrum, which does not resemble the noise-like spectra of ice met in the known literature. The proposed theory satisfactorily agrees with the experimental ice spectrum measured at - 7 degrees C. The calculated longitudinal optic-transverse optic (LO-TO) splitting occurring at approximate to 250 cm(-1) qualitatively agrees with the measured data. (c) 2004 Elsevier B.V. All rights reserved.

Specific and unspecific interactions in polar fluids in view of wideband dielectric far-infrared spectra

A simple molecular analytical theory of dielectric relaxation in strongly polar fluids is considered in terms of a semi- phenomenological approach. Theoretical spectra epsilon(v), a(v) of complex permittivity and absorption coefficient are fully determined by a form of intermolecular potential well, in which a dipole reorients. In a recent publication by VI. Gaiduk, O.F. Nielsen, and T.S. Perova [J. Molliq 95 (1002) 1-25] the wideband spectra of liquid H2O and D2O were described in terms of a composite model comprising the rectangular and the cosine squared potential wells. Much better results are achieved in this work, where the rectangular well is replaced by a well with a rounded bottom termed the hat-curved well. The spectrum of the auto-correlation function (ACF) is calculated for such a potential. The proposed theory of a composite model, comprising hat-curved and parabolic wells, is applied for liquid water. This model is capable for describing the Debye relaxation region, the second relaxation region in the submillimeter wavelength range, and the far infra-red (FIR) e(v), a(v) spectra, where an intense librational band and an additional weak band are placed, respectively, near 700 cm(-1) and 200 cm(-1). The latter band reflects the features of so-called specific (viz. directly related to H-bonds) interactions and the former band reflects the features of unspecific interactions. The physical mechanisms connected with these types of interactions are discussed in terms of two relevant types of water structure (types of molecular rotation). The proposed theory is also applied to a non-associated liquid in terms of one hat-curved potential well. (C) 2004 Elsevier B.V. All rights reserved.

Semiclassical treatment of Wannier's theory when the exponent diverges

We consider a non-standard application of the Wannier model. A physical example is the single ionization of a hydrogenic beryllium ion with a fully stripped beryllium ion, where the ratio of the charge of the third particle to the charges of the escaping particles is 1/4; we investigate the single ionization by an electron of an atom comprising an electron and a nucleus of charge 1/4. An infinite exponent is obtained suggesting that this process is not tractable within the Wannier model. A modified version of Crothers' uniform semiclassical wavefunction for the outgoing particles has been adopted, since the Wannier exponents and are infinite for an effective charge of Z = 1/4. We use Bessel functions to describe the Peterkop functions u and u and derive a new turning point ?. Since u is well behaved at infinity, there exists only the singularity in u at infinity, thus we employ a one- (rather than two-) dimensional change of dependent variable, ensuring that a uniform solution is obtained that avoids semiclassical breakdown on the Wannier ridge. The regularized final-state asymptotic wavefunction is employed, along with a continuum-distorted-wave approximation for the initial-state wavefunction to obtain total cross sections on an absolute scale. © 2006 IOP Publishing Ltd.

Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid-fluid equilibrium, with a metastable fluid-fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76x10(6) MD steps. By progressively reducing the system temperature below the solid-fluid line, we first observe the metastable fluid-fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the predictions of classical nucleation theory. The critical nucleus is identified, and its structure is found to be predominantly fcc. Following nucleation, the solid phase grows steadily across the system, incorporating a large number of localized and extended defects. We discuss the relaxation processes taking place both during and after the crystallization stage. The relevance of our simulation for the kinetics of protein crystallization under normal experimental conditions is discussed. (C) 2002 American Institute of Physics.

Determining the electronic structure and chemical potentials of molecules in solution

A simulation scheme is proposed for determining the excess chemical potential of a substance in solution. First, a Monte Carlo simulation is performed with classical models for solute and solvent molecules. A representative sample of these configurations is then used in a hybrid quantum/classical (QM/MM) calculation, where the solute is treated quantum-mechanically, and the average electronic structure is used to construct an improved classical model. This procedure is iterated to self-consistency in the classical model, which in practice is attained in one or two steps, depending on the quality of the initial guess. The excess free energy of the molecule within the QM/MM approach is determined relative to the classical model using thermodynamic perturbation theory with a cumulant expansion. The procedure provides a method of constructing classical point charge models appropriate for the solution and gives a measure of the importance of solvent fluctuations.

Ab initio simulation of charged slabs at constant chemical potential

We present a practical scheme for performing ab initio supercell calculations of charged slabs at constant electron chemical potential mu, rather than at constant number of electrons N-e. To this end, we define the chemical potential relative to a plane (or "reference electrode") at a finite distance from the slab (the distance should reflect the particular geometry of the situation being modeled). To avoid a net charge in the supercell, and thus make possible a standard supercell calculation, we restore the electroneutrality of the periodically repeated unit by means of a compensating charge, whose contribution to the total energy and potential is subtracted afterwards. The "constant mu" mode enables one to perform supercell calculation on slabs, where the slab is kept at a fixed potential relative to the reference electrode. We expect this to be useful in modeling many experimental situations, especially in electro-chemistry. (C) 2001 American Institute of Physics.

Photodiodes array as potential diagnostic for measuring short bursts of K-alpha radiation from Ti targets irradiated with 45 fs laser pulses

We report on a study comparing absolute K-alpha yield from Ti foils measured with a calibrated system of an X-ray CCD coupled to a curved LiF Von-Hamos crystal spectrometer to the difference in the signals measured simultaneously with two similar photodiodes fitted with two different filters. Our data indicate that a combination of photodiodes with different filters could be developed into an alternative and inexpensive diagnostic for monitoring single shot pulsed emission in a narrow band of X-ray region.

Potential role of intensity-modulated radiotherapy in the treatment of tumors of the maxillary sinus.

To assess 3-dimensional conformal radiotherapy (3D-CRT) and intensity-modulated radiotherapy (IMRT) techniques to see whether doses to critical structures could be reduced while maintaining planning target volume (PTV) coverage in patients receiving conventional radiotherapy (RT) for carcinoma of the maxillary sinus because of the risk of radiation-induced complications, particularly visual loss. Six patients who had recently received conventional RT for carcinoma of the maxillary sinus were studied. Conventional RT, 3D-CRT, and step-and-shoot IMRT plans were prepared using the same 2-field arrangement. The effect of reducing the number of segments in the IMRT beams was investigated. 3D-CRT and IMRT reduced the brain and ipsilateral parotid gland doses compared with the conventional plans. IMRT reduced doses to both optic nerves; for the contralateral optic nerve, 15-segment IMRT plans delivered an average maximal dose of 56.4 Gy (range 53.9–59.3) compared with 65.7 Gy (range 65.3–65.9) and 64.2 Gy (range 61.4–65.6) for conventional RT and 3D-CRT, respectively. IMRT also gave improved PTV homogeneity and improved coverage, with an average of 8.5% (range 7.0–11.7%) of the volume receiving

Structural and chemical embrittlement of grain boundaries by impurities: A general theory and first-principles calculations for copper

First-principles calculations of the Sigma 5(310)[001] symmetric tilt grain boundary in Cu with Bi, Na, and Ag substitutional impurities provide evidence that in the phenomenon of Bi embrittlement of Cu grain boundaries electronic effects do not play a major role; on the contrary, the embrittlement is mostly a structural or "size" effect. Na is predicted to be nearly as good an embrittler as Bi, whereas Ag does not embrittle the boundary in agreement with experiment. While we reject the prevailing view that "electronic" effects (i.e., charge transfer) are responsible for embrittlement, we do not exclude the role of chemistry. However, numerical results show a striking equivalence between the alkali metal Na and the semimetal Bi, small differences being accounted for by their contrasting "size" and "softness" (defined here). In order to separate structural and chemical effects unambiguously if not uniquely, we model the embrittlement process by taking the system of grain boundary and free surfaces through a sequence of precisely defined gedanken processes; each of these representing a putative mechanism. We thereby identify three mechanisms of embrittlement by substitutional impurities, two of which survive in the case of embrittlement or cohesion enhancement by interstitials. Two of the three are purely structural and the third contains both structural and chemical elements that by their very nature cannot be further unraveled. We are able to take the systems we study through each of these stages by explicit computer simulations and assess the contribution of each to the net reduction in intergranular cohesion. The conclusion we reach is that embrittlement by both Bi and Na is almost exclusively structural in origin; that is, the embrittlement is a size effect.

Single-ionization of helium at Ti:sapphire wavelengths: rates and scaling laws

We present a numerical and theoretical study of intense-field single-electron ionization of helium at 390 nm and 780 nm. Accurate ionization rates (over an intensity range of (0.175-34) X10^14 W/ cm^2 at 390 nm, and (0.275 - 14.4) X 10^14 W /cm^2 at 780 nm) are obtained from full-dimensionality integrations of the time-dependent helium-laser Schroedinger equation. We show that the power law of lowest order perturbation theory, modified with a ponderomotive-shifted ionization potential, is capable of modelling the ionization rates over an intensity range that extends up to two orders of magnitude higher than that applicable to perturbation theory alone. Writing the modified perturbation theory in terms of scaled wavelength and intensity variables, we obtain to first approximation a single ionization law for both the 390 nm and 780 nm cases. To model the data in the high intensity limit as well as in the low, a new function is introduced for the rate. This function has, in part, a resemblance to that derived from tunnelling theory but, importantly, retains the correct frequency-dependence and scaling behaviour derived from the perturbative-like models at lower intensities. Comparison with the predictions of classical ADK tunnelling theory confirms that ADK performs poorly in the frequency and intensity domain treated here.