Molecular dynamics simulations of aqueous urea solutions: Study of dimer stability and solution structure, and calculation of the total nitrogen radial distribution function GN(r

Molecular dynamics simulations have been performed in order to study the structure of two molal urea solutions in D2O. Several initial dimer configurations were considered for an adequate sampling of phase space. Eventually all of them appeared to be unstable, when system size and periodic boundary conditions are chosen properly, even after a very careful equilibration. The total nitrogen scattering function GN(r), calculated from these simulations, is in good agreement with neutron scattering experiments when both intra- and intermolecular correlations are considered and the experimental truncation ripples are introduced by a Fourier transform of GN(r) back and forth. The simple pair potential model that we used gives results in good agreement with experiments and with a much more involved potential model, recently described in the literature [J. Chem. Phys. 95, 8419 (1991)]

From wave function to crystal morphology: application to urea and alpha-glycine

In this paper the relation between the molecular electron density distribution and the crystal growth morphology is investigated. Accurate charge densities derived from ab initio quantum chemical calculations were partitioned into multipole moments, to calculate the electrostatic contribution to the intermolecular interaction energy. For urea and alpha-glycine the F-faces or connected nets were determined according to the Hartman-Perdok PBC theory. From attachment energy and critical Ising temperature calculations, theoretical growth forms were constructed using different atom-atom potential models. These were compared to the Donnay-Harker model, equilibrium form and experimental growth forms. In the case of alpha-glycine, the theoretical growth forms are in good agreement with crystals grown from aqueous solution. Crystals obtained by sublimation seem to show some faces which are not F-faces sensu stricto

Molecular dynamics simulation of amphiphilic membrane and wormlike micelles: a multi-scale modelling approach to the design of visco-elastic surfactant solutions

Bilayer membranes and wormlike micelles have been studied using molecular–dynamics simulations. The structure of the worm is analysed in terms of radial density distribution functions, and mechanical properties such as the elastic modulus are calculated. From an analysis of the fluctuation spectra of the tensionless states, we have calculated bending rigidities. Micelles consisting of coarse–grained (CG) model surfactants are studied in order to map the properties of the atomistic micelle. We optimize the CG model with respect to the structure factor S(q) of the atomistic micelle. The mechanical properties thus obtained will be used as input for a mesoscopic model of wormlike micelles where the persistence length is the smallest length–scale

Molecular-dynamics simulations of interfaces between water and crystalline urea

Molecular-dynamics simulations of several water-crystalline urea interfaces have been performed. The structure and dynamics of water close to the urea crystal surface are discussed in terms of density profiles, positional and orientational distribution functions, and diffusion coefficients. The water structure close to the interface is strongly determined by the structure of the crystal surface: the (001) and (111) interfaces reveal strong adsorption of water while the (110) and () interfaces do so to a lesser extent. Assuming that the growth rate of a specific crystal face decreases with increasing solvent adsorption, the appearance of only (111) on the urea growth form is predicted. We argue that on the other hand the dominance of (110) over (001) cannot be explained using a simple layer growth model

Foundations of Dissipative Particle Dynamics

We derive a mesoscopic modeling and simulation technique that is very close to the technique known as dissipative particle dynamics. The model is derived from molecular dynamics by means of a systematic coarse-graining procedure. Thus the rules governing our new form of dissipative particle dynamics reflect the underlying molecular dynamics; in particular all the underlying conservation laws carry over from the microscopic to the mesoscopic descriptions. Whereas previously the dissipative particles were spheres of fixed size and mass, now they are defined as cells on a Voronoi lattice with variable masses and sizes. This Voronoi lattice arises naturally from the coarse-graining procedure which may be applied iteratively and thus represents a form of renormalisation-group mapping. It enables us to select any desired local scale for the mesoscopic description of a given problem. Indeed, the method may be used to deal with situations in which several different length scales are simultaneously present. Simulations carried out with the present scheme show good agreement with theoretical predictions for the equilibrium behavior.Comment: 18 pages, 7 figure

Convex hull approach for determining rock representative elementary volume for multiple petrophysical parameters using pore-scale imaging and Lattice-Boltzmann modelling

In the last decade, the study of fluid flow in porous media has developed considerably due to the combination of X-ray Micro Computed Tomography (micro-CT) and advances in computational methods for solving complex fluid flow equations directly or indirectly on reconstructed three-dimensional pore space images. In this study, we calculate porosity and single phase permeability using micro-CT imaging and Lattice Boltzmann (LB) simulations for 8 different porous media: beadpacks (with bead sizes 50 µm and 350 µm), sandpacks (LV60 and HST95), sandstones (Berea, Clashach and Doddington) and a carbonate (Ketton). Combining the observed porosity and calculated single phase permeability, we shed new light on the existence and size of the Representative Element of Volume (REV) capturing the different scales of heterogeneity from the pore-scale imaging. Our study applies the concept of the ‘Convex Hull’ to calculate the REV by considering the two main macroscopic petrophysical parameters, porosity and single phase permeability, simultaneously. The shape of the hull can be used to identify strong correlation between the parameters or greatly differing convergence rates. To further enhance computational efficiency we note that the area of the convex hull (for well-chosen parameters such as the log of the permeability and the porosity) decays exponentially with sub-sample size so that only a few small simulations are needed to determine the system size needed to calculate the parameters to high accuracy (small convex hull area). Finally we propose using a characteristic length such as the pore size to choose an efficient absolute voxel size for the numerical rock

Search for direct production of charginos and neutralinos in events with three leptons and missing transverse momentum in √s = 7 TeV pp collisions with the ATLAS detector

A search for the direct production of charginos and neutralinos in final states with three electrons or muons and missing transverse momentum is presented. The analysis is based on 4.7 fb−1 of proton–proton collision data delivered by the Large Hadron Collider and recorded with the ATLAS detector. Observations are consistent with Standard Model expectations in three signal regions that are either depleted or enriched in Z-boson decays. Upper limits at 95% confidence level are set in R-parity conserving phenomenological minimal supersymmetric models and in simplified models, significantly extending previous results

Jet size dependence of single jet suppression in lead-lead collisions at sqrt(s(NN)) = 2.76 TeV with the ATLAS detector at the LHC

Measurements of inclusive jet suppression in heavy ion collisions at the LHC provide direct sensitivity to the physics of jet quenching. In a sample of lead-lead collisions at sqrt(s) = 2.76 TeV corresponding to an integrated luminosity of approximately 7 inverse microbarns, ATLAS has measured jets with a calorimeter over the pseudorapidity interval |eta| < 2.1 and over the transverse momentum range 38 < pT < 210 GeV. Jets were reconstructed using the anti-kt algorithm with values for the distance parameter that determines the nominal jet radius of R = 0.2, 0.3, 0.4 and 0.5. The centrality dependence of the jet yield is characterized by the jet "central-to-peripheral ratio," Rcp. Jet production is found to be suppressed by approximately a factor of two in the 10% most central collisions relative to peripheral collisions. Rcp varies smoothly with centrality as characterized by the number of participating nucleons. The observed suppression is only weakly dependent on jet radius and transverse momentum. These results provide the first direct measurement of inclusive jet suppression in heavy ion collisions and complement previous measurements of dijet transverse energy imbalance at the LHC.Comment: 15 pages plus author list (30 pages total), 8 figures, 2 tables, submitted to Physics Letters B. All figures including auxiliary figures are available at http://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/HION-2011-02

Operation and performance of the ATLAS semiconductor tracker

The semiconductor tracker is a silicon microstrip detector forming part of the inner tracking system of the ATLAS experiment at the LHC. The operation and performance of the semiconductor tracker during the first years of LHC running are described. More than 99% of the detector modules were operational during this period, with an average intrinsic hit efficiency of (99.74±0.04)%. The evolution of the noise occupancy is discussed, and measurements of the Lorentz angle, δ-ray production and energy loss presented. The alignment of the detector is found to be stable at the few-micron level over long periods of time. Radiation damage measurements, which include the evolution of detector leakage currents, are found to be consistent with predictions and are used in the verification of radiation background simulations