Exact field-driven interface dynamics in the two-dimensional stochastic Ising model with helicoidal boundary conditions

We investigate the interface dynamics of the two-dimensional stochastic Ising model in an external field under helicoidal boundary conditions. At sufficiently low temperatures and fields, the dynamics of the interface is described by an exactly solvable high-spin asymmetric quantum Hamiltonian that is the infinitesimal generator of the zero range process. Generally, the critical dynamics of the interface fluctuations is in the Kardar-Parisi-Zhang universality class of critical behavior. We remark that a whole family of RSOS interface models similar to the Ising interface model investigated here can be described by exactly solvable restricted high-spin quantum XXZ-type Hamiltonians. (C) 2012 Elsevier B.V. All rights reserved.

Directed and elliptic flow of charged particles in Cu plus Cu collisions at root s(NN)=22.4 GeV

This paper reports results for directed flow v(1) and elliptic flow v(2) of charged particles in Cu + Cu collisions at root s(NN) = 22.4 GeV at the Relativistic Heavy Ion Collider. The measurements are for the 0-60% most central collisions, using charged particles observed in the STAR detector. Our measurements extend to 22.4-GeV Cu + Cu collisions the prior observation that v1 is independent of the system size at 62.4 and 200 GeV and also extend the scaling of v(1) with eta/y(beam) to this system. The measured v(2)(p(T)) in Cu + Cu collisions is similar for root s(NN) throughout the range 22.4 to 200 GeV. We also report a comparison with results from transport model (ultrarelativistic quantum molecular dynamics and multiphase transport model) calculations. The model results do not agree quantitatively with the measured v(1)(eta), v(2)(p(T)), and v(2)(eta).

Alfven eigenmode properties and dynamics in the TJ-II stellarator

Alfven eigenmodes (AEs) were studied in neutral beam injection (NBI) heated plasmas in the TJ-II stellarator using a heavy ion beam probe (HIBP) in the core, and by Langmuir (LP) and Mirnov probes (MP) at the edge. AEs were detected over the whole plasma radius by the HIBP with a spatial resolution of about 1 cm. AE-induced oscillations were detected in the plasma density n(e), electric potential phi and poloidal magnetic field B-pol with frequencies 50 kHz < f(AE) < 300 kHz. The LP, MP and HIBP data showed a high level of coherency for specific branches of AEs. Poloidal mode wave-vectors k(theta), mode numbers m (m < 8) and propagation velocities V-theta similar to 30 km s(-1) were detected for various branches of AEs, having different radial locations. When the density rose due to NBI fuelling, the AE frequency decreased as predicted by the Alfven law f(AE) similar to n(e)(-1/2). During the AE frequency decay the following new AE features were observed: (i) the poloidal wave-vector k(theta) and mode number m remained constant, (ii) the cross-phases between the oscillations in B-pol, n(e) and electric potential remained constant, having an individual value for each AE branch, (iii) V-theta decreased proportional to the AE frequency. The interaction of the AEs with the bulk (thermal) plasma resulted in clearly pronounced quasi-coherent peaks in the electrostatic turbulent particle flux spectra. Various AE branches exhibited different contributions to the particle flux: outward, inward and also zero, depending on the phase relations between the oscillations in E-pol and n(e), which are specific for each branch. A comparison with MHD mode modelling indicated that some of the more prominent frequency branches can be identified as radially extended helical AEs.

Analysis of energy dissipation in stirred suspension polymerisation reactors using computational fluid dynamics

This work evaluates the spatial distribution of normalised rates of droplet breakage and droplet coalescence in liquidliquid dispersions maintained in agitated tanks at operation conditions normally used to perform suspension polymerisation reactions. Particularly, simulations are performed with multiphase computational fluid dynamics (CFD) models to represent the flow field in liquidliquid styrene suspension polymerisation reactors for the first time. CFD tools are used first to compute the spatial distribution of the turbulent energy dissipation rates (e) inside the reaction vessel; afterwards, normalised rates of droplet breakage and particle coalescence are computed as functions of e. Surprisingly, multiphase simulations showed that the rates of energy dissipation can be very high near the free vortex surfaces, which has been completely neglected in previous works. The obtained results indicate the existence of extremely large energy dissipation gradients inside the vessel, so that particle breakage occurs primarily in very small regions that surround the impeller and the free vortex surface, while particle coalescence takes place in the liquid bulk. As a consequence, particle breakage should be regarded as an independent source term or a boundary phenomenon. Based on the obtained results, it can be very difficult to justify the use of isotropic assumptions to formulate particle population balances in similar systems, even when multiple compartment models are used to describe the fluid dynamic behaviour of the agitated vessel. (C) 2011 Canadian Society for Chemical Engineering

Spin coherence generation in negatively charged self-assembled (In,Ga)As quantum dots by pumping excited trion states

Spin coherence generation in an ensemble of negatively charged (In,Ga)As/GaAs quantum dots was investigated by picosecond time-resolved pump-probe spectroscopy measuring ellipticity. Robust coherence of the ground-state electron spins is generated by pumping excited charged exciton (trion) states. The phase of the coherent state, as evidenced by the spin ensemble precession about an external magnetic field, varies relative to spin coherence generation resonant with the ground state. The phase variation depends on the pump photon energy. It is determined by (a) pumping dominantly either singlet or triplet excited states, leading to a phase inversion, and (b) the subsequent carrier relaxation into the ground states. From the dependence of the precession phase and the measured g factors, information about the quantum dot shell splitting and the exchange energy splitting between triplet and singlet states can be extracted in the ensemble.

Anomalous centrality evolution of two-particle angular correlations from Au-Au collisions at root s(NN)=62 and 200 GeV

We present two-dimensional (2D) two-particle angular correlations measured with the STAR detector on relative pseudorapidity eta and azimuth phi for charged particles from Au-Au collisions at root s(NN) = 62 and 200 GeV with transverse momentum p(t) >= 0.15 GeV/c, vertical bar eta vertical bar <= 1, and 2 pi in azimuth. Observed correlations include a same-side (relative azimuth <pi/2) 2D peak, a closely related away-side azimuth dipole, and an azimuth quadrupole conventionally associated with elliptic flow. The same-side 2D peak and away-side dipole are explained by semihard parton scattering and fragmentation (minijets) in proton-proton and peripheral nucleus-nucleus collisions. Those structures follow N-N binary-collision scaling in Au-Au collisions until midcentrality, where a transition to a qualitatively different centrality trend occurs within one 10% centrality bin. Above the transition point the number of same-side and away-side correlated pairs increases rapidly relative to binary-collision scaling, the eta width of the same-side 2D peak also increases rapidly (eta elongation), and the phi width actually decreases significantly. Those centrality trends are in marked contrast with conventional expectations for jet quenching in a dense medium. The observed centrality trends are compared to perturbative QCD predictions computed in HIJING, which serve as a theoretical baseline, and to the expected trends for semihard parton scattering and fragmentation in a thermalized opaque medium predicted by theoretical calculations and phenomenological models. We are unable to reconcile a semihard parton scattering and fragmentation origin for the observed correlation structure and centrality trends with heavy-ion collision scenarios that invoke rapid parton thermalization. If the collision system turns out to be effectively opaque to few-GeV partons the present observations would be inconsistent with the minijet picture discussed here. DOI: 10.1103/PhysRevC.86.064902

Photodynamic Efficiency of Cationic meso-Porphyrins at Lipid Bilayers: Insights from Molecular Dynamics Simulations

Porphyrin derivatives have applications as photoactive drugs in photodynamic therapy. However, little is known about their interactions with phospholipid membranes at the molecular level. We employed molecular dynamics simulations to model the binding between a series of cationic meso-(N-methyl-4-pyridinium)phenylporphyrins and anionic phosphatidylglycerol lipid bilayers. This was done in the presence of molecular oxygen within the membrane. The ability of various porphyrins to cause photodamage was quantified in terms of their immersion depth and degree of exposition to a higher oxygen concentration inside the membrane. Simulations showed that the photodynamic efficiency could be improved as the number of hydrophobic phenyl substituents attached to the porphyrinic ring increased. In the specific case of porphyrins containing two hydrophobic and two charged substituents, the cis isomer was significantly more efficient than the trans. These results correlate well with previous experimental observations. They highlight the importance of both the total charge and amphiphilicity of the photosensitizer for its performance in photodynamic therapy.

Controlling chaos in wave-particle interactions

We analyze the behavior of a relativistic particle moving under the influence of a uniform magnetic field and a stationary electrostatic wave. We work with a set of pulsed waves that allows us to obtain an exact map for the system. We also use a method of control for near-integrable Hamiltonians that consists of the addition of a small and simple control term to the system. This control term creates invariant tori in phase space that prevent chaos from spreading to large regions, making the controlled dynamics more regular. We show numerically that the control term just slightly modifies the system but is able to drastically reduce chaos with a low additional cost of energy. Moreover, we discuss how the control of chaos and the consequent recovery of regular trajectories in phase space are useful to improve regular particle acceleration.

Inclusive charged hadron elliptic flow in Au+Au collisions at root s(NN)=7.7-39 GeV

A systematic study is presented for centrality, transverse momentum (p(T)), and pseudorapidity (eta) dependence of the inclusive charged hadron elliptic flow (v(2)) at midrapidity (vertical bar eta vertical bar < 1.0) in Au + Au collisions at root s(NN) = 7.7, 11.5, 19.6, 27, and 39 GeV. The results obtained with different methods, including correlations with the event plane reconstructed in a region separated by a large pseudorapidity gap and four-particle cumulants (v(2){4}), are presented to investigate nonflow correlations and v(2) fluctuations. We observe that the difference between v(2){2} and v(2){4} is smaller at the lower collision energies. Values of v(2), scaled by the initial coordinate space eccentricity, v(2)/epsilon, as a function of p(T) are larger in more central collisions, suggesting stronger collective flow develops in more central collisions, similar to the results at higher collision energies. These results are compared to measurements at higher energies at the Relativistic Heavy Ion Collider (root s(NN) = 62.4 and 200 GeV) and at the Large Hadron Collider (Pb + Pb collisions at root s(NN) = 2.76 TeV). The v(2)(pT) values for fixed pT rise with increasing collision energy within the pT range studied (<2 GeV/c). A comparison to viscous hydrodynamic simulations is made to potentially help understand the energy dependence of v(2)(pT). We also compare the v(2) results to UrQMD and AMPT transport model calculations, and physics implications on the dominance of partonic versus hadronic phases in the system created at beam energy scan energies are discussed.

Hydrodynamic analysis of particle collection efficiency: comparing downflow and upflow filtration

Models of the filtration phenomenon describe the mass balance in bed filtration in terms of particle removal mechanisms, and allow for the determination of global particle removal efficiencies. These models are defined in terms of the geometry and characteristic elements of granule collectors, particles and fluid, and also the composition of the balance of forces that act in the particle collector system. This work analyzes particles collection efficiency comparing downflow and upflow direct filtration, taking into account the contribution of the gravitational factor of the settling removal efficiency in future proposal of initial collection efficiency models for upflow filtration. A qualitative analysis is also made of the proposal for the collection efficiency models for particle removal in direct downflow and upflow filtration using a Computational Fluid Dynamics (CFD) tool. This analysis showed a strong influence of gravitational factor in initial collection efficiency (t = 0) of particles, as well as the reasons of their values to be smaller for upflow filtration in comparison with the downflow filtration.

Particle Acceleration in Turbulence and Weakly Stochastic Reconnection

In this Letter we analyze the energy distribution evolution of test particles injected in three dimensional (3D) magnetohydrodynamic (MHD) simulations of different magnetic reconnection configurations. When considering a single Sweet-Parker topology, the particles accelerate predominantly through a first-order Fermi process, as predicted in [3] and demonstrated numerically in [8]. When turbulence is included within the current sheet, the acceleration rate is highly enhanced, because reconnection becomes fast and independent of resistivity [4,11] and allows the formation of a thick volume filled with multiple simultaneously reconnecting magnetic fluxes. Charged particles trapped within this volume suffer several head-on scatterings with the contracting magnetic fluctuations, which significantly increase the acceleration rate and results in a first-order Fermi process. For comparison, we also tested acceleration in MHD turbulence, where particles suffer collisions with approaching and receding magnetic irregularities, resulting in a reduced acceleration rate. We argue that the dominant acceleration mechanism approaches a second order Fermi process in this case.

Model for the Dynamics of the Cytoskeleton.

Particle tracking of microbeads attached to the cytoskeleton (CSK) reveals an intermittent dynamic. The mean squared displacement (MSD) is subdiffusive for small Δt and superdiffusive for large Δt, which are associated with periods of traps and periods of jumps respectively. The analysis of the displacements has shown a non-Gaussian behavior, what is indicative of an active motion, classifying the cells as a far from equilibrium material. Using Langevin dynamics, we reconstruct the dynamic of the CSK. The model is based on the bundles of actin filaments that link themself with the bead RGD coating, trapping it in an harmonic potential. We consider a one- dimensional motion of a particle, neglecting inertial effects (over-damped Langevin dynamics). The resultant force is decomposed in friction force, elastic force and random force, which is used as white noise representing the effect due to molecular agitation. These description until now shows a static situation where the bead performed a random walk in an elastic potential. In order to modeling the active remodeling of the CSK, we vary the equilibrium position of the potential. Inserting a motion in the well center, we change the equilibrium position linearly with time with constant velocity. The result found exhibits a MSD versus time ’tau’ with three regimes. The first regime is when ‘tau’ < ‘tau IND 0’, where ‘tau IND 0’ is the relaxation time, representing the thermal motion. At this regime the particle can diffuse freely. The second regime is a plateau, ‘tau IND 0’ < ‘tau’ < ‘tau IND 1’, representing the particle caged in the potential. Here, ‘tau IND 1’ is a characteristic time that limit the confinement period. And the third regime, ‘tau’ > ‘tau IND 1’, is when the particles are in the superdiffusive behavior. This is where most of the experiments are performed, under 20 frames per second (FPS), thus there is no experimental evidence that support the first regime. We are currently performing experiments with high frequency, up to 100 FPS, attempting to visualize this diffusive behavior. Beside the first regime, our simple model can reproduce MSD curves similar to what has been found experimentally, which can be helpful to understanding CSK structure and properties.

New insights into POC dynamics in the subtropical northeast Atlantic Ocean

[EN]The capacity of the ocean to sequester atmospheric carbon (CO2) depends to a large extent on the dynamics of biogenic carbon in the water column. However, most current global and regional estimates of carbon balances are solely based on particles collected with drifting and moored sediment traps. As a consequence, construction of ocean carbon budgets has long been guided by the simplification introduced by sediment traps, which give a 1D vision of the whole picture. In this thesis we have assessed a quantitative analysis of the flux magnitude and the mechanisms of transport of the whole particle spectrum (suspended, slowly-sinking and sinking particles).

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

Vortex dynamics in the high-temperature superconductors YBa 2 Cu 3 O 7 and Bi 2 Sr 2 CaCu 2 O 8 + delta in low- and high-dissipative transport measurements

This thesis investigates phenomena of vortex dynamics in type II superconductors depending on the dimensionality of the flux-line system and the strength of the driving force. In the low dissipative regime of Bi_2Sr_2CaCu_2O_{8+delta} (BSCCO) the influence of oxygen stoichiometry on flux-line tension was examined. An entanglement crossover of the vortex system at low magnetic fields was identified and a comprehensive B-T phase diagram of solid and fluid phases derived.In YBa_2Cu_3O_7 (YBCO) extremely long (>100 mm) high-quality measurement bridges allowed to extend the electric-field window in transport measurements by up to three orders of magnitude. Complementing analyses of the data conclusively produced dynamic exponents of the glass transition z~9 considerably higher than theoretically predicted and previously reported. In high-dissipative measurements a voltage instability appearing in the current-voltage characteristics of type II superconductors was observed for the first time in BSCCO and shown to result from a Larkin-Ovchinnikov flux-flow vortex instability under the influence of quasi-particle heating. However, in an analogous investigation of YBCO the instability was found to appear only in the temperature and magnetic-field regime of the vortex-glass state. Rapid-pulse measurements fully confirmed this correlation of vortex glass and instability in YBCO and revealed a constant rise time (~µs).