Computational Methods for Studies of Quantum Dots and Boson Condensates

This thesis presents ab initio studies of two kinds of physical systems, quantum dots and bosons, using two program packages of which the bosonic one has mainly been developed by the author. The implemented models, \emph{i.e.}, configuration interaction (CI) and coupled cluster (CC) take the correlated motion of the particles into account, and provide a hierarchy of computational schemes, on top of which the exact solution, within the limit of the single-particle basis set, is obtained. The theory underlying the models is presented in some detail, in order to provide insight into the approximations made and the circumstances under which they hold. Some of the computational methods are also highlighted. In the final sections the results are summarized. The CI and CC calculations on multiexciton complexes in self-assembled semiconductor quantum dots are presented and compared, along with radiative and non-radiative transition rates. Full CI calculations on quantum rings and double quantum rings are also presented. In the latter case, experimental and theoretical results from the literature are re-examined and an alternative explanation for the reported photoluminescence spectra is found. The boson program is first applied on a fictitious model system consisting of bosonic electrons in a central Coulomb field for which CI at the singles and doubles level is found to account for almost all of the correlation energy. Finally, the boson program is employed to study Bose-Einstein condensates confined in different anisotropic trap potentials. The effects of the anisotropy on the relative correlation energy is examined, as well as the effect of varying the interaction potential.}

The velocity dispersion of giant molecular clouds. II - Mathematical and numerical refinements

The energy input to giant molecular clouds is recalculated, using the proper linearized equations of motion, including the Coriolis force and allowing for changes in the guiding center. Perturbation theory yields a result in the limit of distant encounters and small initial epicyclic amplitudes. Direct integration of the motion equations allows the strong encounter regime to be studied. The present perturbation theory result differs by a factor of order unity from that of Jog and Ostriker (1988). The result of present numerical integrations for the 2D (planar) velocity dispersion is presented. The accretion rate for a molecular cloud in the Galactic disk is calculated.

Structure of Relativistic States

QCD factorization in the Bjorken limit allows to separate the long-distance physics from the hard subprocess. At leading twist, only one parton in each hadron is coherent with the hard subprocess. Higher twist effects increase as one of the active partons carries most of the longitudinal momentum of the hadron, x -> 1. In the Drell-Yan process \pi N -> \mu^- mu^+ + X, the polarization of the virtual photon is observed to change to longitudinal when the photon carries x_F > 0.6 of the pion. I define and study the Berger-Brodsky limit of Q^2 -> \infty with Q^2(1-x) fixed. A new kind of factorization holds in the Drell-Yan process in this limit, in which both pion valence quarks are coherent with the hard subprocess, the virtual photon is longitudinal rather than transverse, and the cross section is proportional to a multiparton distribution. Generalized parton distributions contain information on the longitudinal momentum and transverse position densities of partons in a hadron. Transverse charge densities are Fourier transforms of the electromagnetic form factors. I discuss the application of these methods to the QED electron, studying the form factors, charge densities and spin distributions of the leading order |e\gamma> Fock state in impact parameter and longitudinal momentum space. I show how the transverse shape of any virtual photon induced process, \gamma^*(q)+i -> f, may be measured. Qualitative arguments concerning the size of such transitions have been previously made in the literature, but without a precise analysis. Properly defined, the amplitudes and the cross section in impact parameter space provide information on the transverse shape of the transition process.

Critical Reappraisal of Colloidal Activity of Clays

Soil properties and their behavior, apart from stress history, are influenced markedly by physicochemical characteristics of the constituent clay and nonclay minerals and their relative proportions. Atterberg limits and Skempton’s colloidal activity, which are simple quantitative parameters, reflect the composite effects of the soil constituents and their interactions with pore fluid. Micromechanistic interpretations of these parameters have been provided in this paper. It has been shown that, in general, the liquid limit of fine-grained soils reflects the physicochemical potential and that each of the factors of Skempton’s colloidal activity are interdependent. It has been illustrated that property correlations with colloidal activity, as well as with Atterberg limits, result in involved interrelationships due to the interdependence of the parameters.

Dynamic Structure Factor of Sponge Phases

We predict the dynamic light scattering intensity S(q,t) for the L3 phase (anomalous isotropic phase) of dilute surfactant solutions. Our results are based on a Landau-Ginzburg approach, which was previously used to explain the observed static structure factor S(q, 0). In the extreme limit of small q, we find a monoexponential decay with marginal or irrelevant hydrodynamic interactions. In most other regimes the decay of S(q,t) is strongly nonexponential; in one case, it is purely algebraic at long times.

Critical Reappraisal of Colloidal Activity of Clays

Soil properties and their behavior, apart from stress history, are influence markedly by physicochemical characteristics of the constituent clay and nonclay minerals and their relative proportions. Atterberg limits and Skempton's colloidal activity, which are simple quantitative parameters, reflect the composite effects of the soil constituents and their interactions with pore fluid. Micromechanistic interpretations of these parameters have been provided in this paper. It has been shown that, in general, the liquid limit of fine-grained soils reflects the physicochemical potential and that each of the factors of Skempton's colloidal activity are interdependent. It has been illustrated that property correlations with colloidal activity, as well as with Atterberg limits, result in involved interrelationships due to the interdependence of the parameters.

Consistency limits behavior of bentonites exposed to sea water

This article examines the changes in interparticle forces brought about on prolonged contact (1 year period) of a bentonite clay with artificial seawater. The study is undertaken with the purpose of identifying the physico-chemical factors that impart a nonswelling character to smectite clays deposited in marine environments. Results show that equilibration of the bentonite clay with artificial seawater (total pore salinity approximately 42 gL-1) for a 1 year period does not lead to any mineralogical changes in the clay specimens; however, their exchangeable cation positions become prominently dominated by magnesium ions. The consistency limits of the seawater-equilibrated bentonite was determined on stepwise leaching to lower salinities. The predominance of diffuse double-layer repulsion forces in the pore salt concentration range of 42 gL-1 to 1.1 gL-1 caused an increase in the liquid limits of the seawater-equilibrated bentonite specimens on reducing the salinity in the corresponding range (42 gL-1 to 1.1 gL-1). The attraction forces, however, prevail over the repulsion forces at salt concentrations <1.1 gL-1 and cause a decrease in liquid limit of the clay specimens with reduction in pore salinity, which is typical of nonswelling clays. The attraction forces cause aggregation of the clay unit layers into domains that break down on sodium saturation of the clay specimens. It is inferred that the physico-chemical factors responsible for the nonswelling character of the seawater-equilibrated bentonite specimens at pore salt concentrations below 1.1 gL-1 are inadequate to explain the nonswelling character of smectite-rich Ariake marine clays. The lower consistency limits of the Ariake marine clays in comparison to the nonswelling character, seawater-equilibrated bentonite specimens is attributed to a relative deficiency of interparticle forces in the Ariake marine clay.

Microbial identification by detection of ligation probes on DNA microarray

Microbes in natural and artificial environments as well as in the human body are a key part of the functional properties of these complex systems. The presence or absence of certain microbial taxa is a correlate of functional status like risk of disease or course of metabolic processes of a microbial community. As microbes are highly diverse and mostly notcultivable, molecular markers like gene sequences are a potential basis for detection and identification of key types. The goal of this thesis was to study molecular methods for identification of microbial DNA in order to develop a tool for analysis of environmental and clinical DNA samples. Particular emphasis was placed on specificity of detection which is a major challenge when analyzing complex microbial communities. The approach taken in this study was the application and optimization of enzymatic ligation of DNA probes coupled with microarray read-out for high-throughput microbial profiling. The results show that fungal phylotypes and human papillomavirus genotypes could be accurately identified from pools of PCR amplicons generated from purified sample DNA. Approximately 1 ng/μl of sample DNA was needed for representative PCR amplification as measured by comparisons between clone sequencing and microarray. A minimum of 0,25 amol/μl of PCR amplicons was detectable from amongst 5 ng/μl of background DNA, suggesting that the detection limit of the test comprising of ligation reaction followed by microarray read-out was approximately 0,04%. Detection from sample DNA directly was shown to be feasible with probes forming a circular molecule upon ligation followed by PCR amplification of the probe. In this approach, the minimum detectable relative amount of target genome was found to be 1% of all genomes in the sample as estimated from 454 deep sequencing results. Signal-to-noise of contact printed microarrays could be improved by using an internal microarray hybridization control oligonucleotide probe together with a computational algorithm. The algorithm was based on identification of a bias in the microarray data and correction of the bias as shown by simulated and real data. The results further suggest semiquantitative detection to be possible by ligation detection, allowing estimation of target abundance in a sample. However, in practise, comprehensive sequence information of full length rRNA genes is needed to support probe design with complex samples. This study shows that DNA microarray has the potential for an accurate microbial diagnostic platform to take advantage of increasing sequence data and to replace traditional, less efficient methods that still dominate routine testing in laboratories. The data suggests that ligation reaction based microarray assay can be optimized to a degree that allows good signal-tonoise and semiquantitative detection.

Path integral description of polymers using fractional Brownian walks

The statistical properties of fractional Brownian walks are used to construct a path integral representation of the conformations of polymers with different degrees of bond correlation. We specifically derive an expression for the distribution function of the chains’ end‐to‐end distance, and evaluate it by several independent methods, including direct evaluation of the discrete limit of the path integral, decomposition into normal modes, and solution of a partial differential equation. The distribution function is found to be Gaussian in the spatial coordinates of the monomer positions, as in the random walk description of the chain, but the contour variables, which specify the location of the monomer along the chain backbone, now depend on an index h, the degree of correlation of the fractional Brownian walk. The special case of h=1/2 corresponds to the random walk. In constructing the normal mode picture of the chain, we conjecture the existence of a theorem regarding the zeros of the Bessel function.

Ultrafast underdamped solvation: Agreement between computer simulation and various theories of solvation dynamics

A theoretical analysis of the three currently popular microscopic theories of solvation dynamics, namely, the dynamic mean spherical approximation (DMSA), the molecular hydrodynamic theory (MHT), and the memory function theory (MFT) is carried out. It is shown that in the underdamped limit of momentum relaxation, all three theories lead to nearly identical results when the translational motions of both the solute ion and the solvent molecules are neglected. In this limit, the theoretical prediction is in almost perfect agreement with the computer simulation results of solvation dynamics in the model Stockmayer liquid. However, the situation changes significantly in the presence of the translational motion of the solvent molecules. In this case, DMSA breaks down but the other two theories correctly predict the acceleration of solvation in agreement with the simulation results. We find that the translational motion of a light solute ion can play an important role in its own solvation. None of the existing theories describe this aspect. A generalization of the extended hydrodynamic theory is presented which, for the first time, includes the contribution of solute motion towards its own solvation dynamics. The extended theory gives excellent agreement with the simulations where solute motion is allowed. It is further shown that in the absence of translation, the memory function theory of Fried and Mukamel can be recovered from the hydrodynamic equations if the wave vector dependent dissipative kernel in the hydrodynamic description is replaced by its long wavelength value. We suggest a convenient memory kernel which is superior to the limiting forms used in earlier descriptions. We also present an alternate, quite general, statistical mechanical expression for the time dependent solvation energy of an ion. This expression has remarkable similarity with that for the translational dielectric friction on a moving ion.

Single-molecule thermodynamics: the heat distribution function of a charged particle in a static magnetic field

Interest in the applicability of fluctuation theorems to the thermodynamics of single molecules in external potentials has recently led to calculations of the work and total entropy distributions of Brownian oscillators in static and time-dependent electromagnetic fields. These calculations, which are based on solutions to a Smoluchowski equation, are not easily extended to a consideration of the other thermodynamic quantity of interest in such systems-the heat exchanges of the particle alone-because of the nonlinear dependence of the heat on a particle's stochastic trajectory. In this paper, we show that a path integral approach provides an exact expression for the distribution of the heat fluctuations of a charged Brownian oscillator in a static magnetic field. This approach is an extension of a similar path integral approach applied earlier by our group to the calculation of the heat distribution function of a trapped Brownian particle, which was found, in the limit of long times, to be consistent with experimental data on the thermal interactions of single micron-sized colloids in a viscous solvent.

Enhanced sensitivity for the differential pulse polarographic determination of molybdenum(VI) based on adsorption catalytic current

A differential pulse polarographic (DPP) method based on the adsorption catalytic current in a medium containing chlorate and 8-hydroxyquinoline (oxine) is suggested for the determination of molybdenum(VI). Experimental conditions such as pH and the composition of supporting electrolyte have been optimized to get a linear calibration graph at trace levels of Mo(VI). The sensitivity for molybdenum can be considerably enhanced by this method. The influence of possible interferences on the catalytic current has been investigated. The sensitivity of the method is compared with those obtained for other DPP methods for molybdenum. A detection limit of 1.0 x 10(-8) mol/L has been found.

Vortex sheet simulation of a plane �canonical� mixing layer

Discrete vortex simulations of the mixing layer carried out in the past have usually involved large induced velocity fluctuations, and thus demanded rather long time-averaging to obtain satisfactory values of Reynolds stresses and third-order moments. This difficulty has been traced here, in part, to the use of discrete vortices to model what in actuality are continuous vortex sheets. We propose here a novel two-dimensional vortex sheet technique for computing mixing layer flow in the limit of infinite Reynolds number. The method divides the vortex sheet into constant-strength linear elements, whose motions are computed using the Biot-Savart law. The downstream far-field is modelled by a steady vorticity distribution derived by application of conical similarity from the solution obtained in a finite computational domain. The boundary condition on the splitter plate is satisfied rigorously using a doublet sheet. The computed large-scale roll-up of the vortex sheet is qualitatively similar to experimentally obtained shadow-graphs of the plane turbulent mixing layer. The mean streamwise velocity profile and the growth rate agree well with experimental data. The presently computed Reynolds stresses and third-order moments are comparable with experimental and previous vortex-dynamical results, without using any external parameter (such as the vortex core-size) of the kind often used in the latter. The computed autocorrelations are qualitatively similar to experimental results along the top and bottom edges of the mixing layer, and show a well-defined periodicity along the centreline. The accuracy of the present computation is independently established by demonstrating negligibly small changes in the five invariants (including the Hamiltonian) in vortex dynamics.

Characteristic water contents of a fine-grained soil-water system

This paper critically appraises the limitations of the liquid-limit water content of clayey soils determined conventionally either by percussion cup or by the cone penetration method. It is shown that the conventional liquid limit and plastic limit are arbitrary, strength-based water contents and that they cannot represent the plasticity limits, and that the state of the soil-water system at the conventional liquid limit does not correspond to a stress-free reference state. The present investigation identifies three characteristic limiting water contents for a soil-water system which have well-defined, unique mechanisms controlling them, namely the free swell limit, settling limit and shrinkage limit. Simple procedures for the determination of the free swell limit and settling limit of natural soils are presented. The settling limit is shown to be the 'real liquid limit' of any clayey soil. With a number of experimental illustrations, it is clearly shown that the settling limit represents the maximum water-holding capacity of clayey soils and that it corresponds to the stress-free reference state.

Magnetization properties of some quantum spin ladders

The experimental realization of various spin ladder systems has prompted their detailed theoretical investigations. Hen we study the evolution of ground-state magnetization with an external magnetic field for two different antiferromagnetic systems: a three-legged spin-1/2 ladder, and a two-legged spin-1/2 ladder with an additional diagonal interaction. The finite system density-matrix renormalization-group method is employed for numerical studies of the three-chain system, and an effective low-energy Hamiltonian is used in the limit of strong interchain coupling to study the two- and three-chain systems. The three-chain system has a magnetization plateau at one-third of the saturation magnetization. The two-chain system has a plateau at zero magnetization due to a gap above the singlet ground state. It also has a plateau at half of the saturation magnetization for a certain range of values of the couplings. We study the regions of transitions between plateaus numerically and analytically, and find that they are described, at first order in a strong-coupling expansion, by an XXZ spin-1/2 chain in a magnetic field; the second-order terms give corrections to the XXZ model, We also study numerically some low-temperature properties of the three-chain system, such as the magnetization, magnetic susceptibility and specific heat. [S0163-1829(99)303001-5].