Catalog footprint

What is connected

60works
27topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

60 published item(s)

preprint2023arXiv

Measurement of $ϕ$-meson production in Cu$+$Au at $\sqrt{s_{_{NN}}}=200$ GeV and U$+$U at $\sqrt{s_{_{NN}}}=193$ GeV

The PHENIX experiment reports systematic measurements at the Relativistic Heavy Ion Collider of $ϕ$-meson production in asymmetric Cu$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV and in U$+$U collisions at $\sqrt{s_{_{NN}}}$=193 GeV. Measurements were performed via the $ϕ\rightarrow K^{+}K^{-}$ decay channel at midrapidity $|η|<0.35$. Features of $ϕ$-meson production measured in Cu$+$Cu, Cu$+$Au, Au$+$Au, and U$+$U collisions were found to not depend on the collision geometry, which was expected because the yields are averaged over the azimuthal angle and follow the expected scaling with nuclear-overlap size. The elliptic flow of the $ϕ$ meson in Cu$+$Au, Au$+$Au, and U$+$U collisions scales with second-order-participant eccentricity and the length scale of the nuclear-overlap region (estimated with the number of participating nucleons). At moderate $p_T$, $ϕ$-meson production measured in Cu$+$Au and U$+$U collisions is consistent with coalescence-model predictions, whereas at high $p_T$ the production is in agreement with expectations for in-medium energy loss of parent partons prior to their fragmentation. The elliptic flow for $ϕ$ mesons measured in Cu$+$Au and U$+$U collisions is well described by a (2+1)D viscous-hydrodynamic model with specific-shear viscosity $η/s=1/4π$.

preprint2022arXiv

Kinematic dependence of azimuthal anisotropies in $p$$+$Au, $d$$+$Au, $^3$He+Au at $\sqrt{s_{_{NN}}}$ = 200 GeV

There is strong evidence for the formation of small droplets of quark-gluon plasma in $p/d/^{3}$He+Au collisions at the Relativistic Heavy Ion Collider (RHIC) and in $p$+$p$/Pb collisions at the Large Hadron Collider. In particular, the analysis of data at RHIC for different geometries obtained by varying the projectile size and shape has proven insightful. In the present analysis, we find excellent agreement with the previously published PHENIX at RHIC results on elliptical and triangular flow with an independent analysis via the two-particle correlation method, which has quite different systematic uncertainties and an independent code base. In addition, the results are extended to other detector combinations with different kinematic (pseudorapidity) coverage. These results provide additional constraints on contributions from nonflow and longitudinal decorrelations.

preprint2022arXiv

Measurement of $ψ(2S)$ nuclear modification at backward and forward rapidity in $p$$+$$p$, $p$$+$Al, and $p$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV

Suppression of the $J/ψ$ nuclear-modification factor has been seen as a trademark signature of final-state effects in large collision systems for decades. In small systems, the nuclear modification was attributed to cold-nuclear-matter effects until the observation of strong differential suppression of the $ψ(2S)$ state in $p/d$$+$$A$ collisions suggested the presence of final-state effects. Results of $J/ψ$ and $ψ(2S)$ measurements in the dimuon decay channel are presented here for $p$$+$$p$, $p$$+$Al, and $p$$+$Au collision systems at $\sqrt{s_{_{NN}}}=200$ GeV. The results are predominantly shown in the form of the nuclear-modification factor, $R_{pA}$, the ratio of the $ψ(2S)$ invariant yield per nucleon-nucleon collision in collisions of proton on target nucleus to that in $p$$+$$p$ collisions. Measurements of the $J/ψ$ and $ψ(2S)$ nuclear-modification factor are compared with shadowing and transport-model predictions, as well as to complementary measurements at Large-Hadron-Collider energies.

preprint2022arXiv

Study of $ϕ$-meson production in $p$$+$Al, $p$$+$Au, $d$$+$Au, and $^3$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV

Small nuclear collisions are mainly sensitive to cold-nuclear-matter effects; however, the collective behavior observed in these collisions shows a hint of hot-nuclear-matter effects. The identified-particle spectra, especially the $ϕ$ mesons which contain strange and antistrange quarks and have a relatively small hadronic-interaction cross section, are a good tool to study these effects. The PHENIX experiment has measured $ϕ$ mesons in a specific set of small collision systems $p$$+$Al, $p$$+$Au, and $^3$He$+$Au, as well as $d$$+$Au [Phys. Rev. C {\bf 83}, 024909 (2011)], at $\sqrt{s_{_{NN}}}=200$ GeV. The transverse-momentum spectra and nuclear-modification factors are presented and compared to theoretical-model predictions. The comparisons with different calculations suggest that quark-gluon plasma may be formed in these small collision systems at $\sqrt{s_{_{NN}}}=200$ GeV. However, the volume and the lifetime of the produced medium may be insufficient for observing strangeness-enhancement and jet-quenching effects. Comparison with calculations suggests that the main production mechanisms of $ϕ$ mesons at midrapidity may be different in $p$$+$Al versus $p/d/$$^3$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. While thermal quark recombination seems to dominate in $p/d/$$^3$He$+$Au collisions, fragmentation seems to be the main production mechanism in $p$$+$Al collisions.

preprint2022arXiv

Systematic study of nuclear effects in $p$$+$Al, $p$$+$Au, $d$$+$Au, and $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV using $π^0$ production

The PHENIX collaboration presents a systematic study of $π^0$ production from $p$$+$$p$, $p$$+$Al, $p$$+$Au, $d$$+$Au, and $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Measurements were performed with different centrality selections as well as the total inelastic, 0%--100%, selection for all collision systems. For 0%--100% collisions, the nuclear modification factors, $R_{xA}$, are consistent with unity for $p_T$ above 8 GeV/$c$, but exhibit an enhancement in peripheral collisions and a suppression in central collisions. The enhancement and suppression characteristics are similar for all systems for the same centrality class. It is shown that for high-$p_T$-$π^0$ production, the nucleons in the $d$ and $^3$He interact mostly independently with the Au nucleus and that the counter intuitive centrality dependence is likely due to a physical correlation between multiplicity and the presence of a hard scattering process. These observations disfavor models where parton energy loss has a significant contribution to nuclear modifications in small systems. Nuclear modifications at lower $p_T$ resemble the Cronin effect -- an increase followed by a peak in central or inelastic collisions and a plateau in peripheral collisions. The peak height has a characteristic ordering by system size as $p$$+$Au $>$ $d$$+$Au $>$ $^{3}$He$+$Au $>$ $p$$+$Al. For collisions with Au ions, current calculations based on initial state cold nuclear matter effects result in the opposite order, suggesting the presence of other contributions to nuclear modifications, in particular at lower $p_T$.

preprint2022arXiv

Transverse-single-spin asymmetries of charged pions at midrapidity in transversely polarized $p{+}p$ collisions at $\sqrt{s}=200$ GeV

In 2015, the PHENIX collaboration has measured single-spin asymmetries for charged pions in transversely polarized proton-proton collisions at the center of mass energy of $\sqrt{s}=200$ GeV. The pions were detected at central rapidities of $|η|<0.35$. The single-spin asymmetries are consistent with zero for each charge individually, as well as consistent with the previously published neutral-pion asymmetries in the same rapidity range. However, they show a slight indication of charge-dependent differences which may suggest a flavor dependence in the underlying mechanisms that create these asymmetries.

preprint2021arXiv

Transverse momentum dependent forward neutron single spin asymmetries in transversely polarized $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV

In 2015, the PHENIX collaboration has measured very forward ($η>6.8$) single-spin asymmetries of inclusive neutrons in transversely polarized proton-proton and proton-nucleus collisions at a center of mass energy of 200 GeV. A previous publication from this data set concentrated on the nuclear dependence of such asymmetries. In this measurement the explicit transverse-momentum dependence of inclusive neutron single spin asymmetries for proton-proton collisions is extracted using a bootstrapping-unfolding technique on the transverse momenta. This explicit transverse-momentum dependence will help improve the understanding of the mechanisms that create these asymmetries.

preprint2020arXiv

$J/ψ$ and $ψ(2S)$ production at forward rapidity in $p$+$p$ collisions at $\sqrt{s}=510$ GeV

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the differential cross section, mean transverse momentum, mean transverse momentum squared of inclusive $J/ψ$ and cross-section ratio of $ψ(2S)$ to $J/ψ$ at forward rapidity in \pp collisions at \sqrts = 510 GeV via the dimuon decay channel. Comparison is made to inclusive $J/ψ$ cross sections measured at \sqrts = 200 GeV and 2.76--13 TeV. The result is also compared to leading-order nonrelativistic QCD calculations coupled to a color-glass-condensate description of the low-$x$ gluons in the proton at low transverse momentum ($p_T$) and to next-to-leading order nonrelativistic QCD calculations for the rest of the $p_T$ range. These calculations overestimate the data at low $p_T$. While consistent with the data within uncertainties above $\approx3$ GeV/$c$, the calculations are systematically below the data. The total cross section times the branching ratio is BR $dσ^{J/ψ}_{pp}/dy (1.2<|y|<2.2, 0<p_T<10~\mbox{GeV/$c$}) =$ 54.3 $\pm$ 0.5 (stat) $\pm$ 5.5 (syst) nb.

preprint2020arXiv

Evolution of pi^0 suppression in Au+Au collisions from sqrt(s_NN) = 39 to 200 GeV

Neutral-pion, pi^0, spectra were measured at midrapidity (|y|<0.35) in Au+Au collisions at sqrt(s_NN) = 39 and 62.4 GeV and compared to earlier measurements at 200 GeV in the 1<p_T<10 GeV/c transverse-momentum (p_T) range. The high-p_T tail is well described by a power law in all cases and the powers decrease significantly with decreasing center-of-mass energy. The change of powers is very similar to that observed in the corresponding p+p-collision spectra. The nuclear-modification factors (R_AA) show significant suppression and a distinct energy dependence at moderate p_T in central collisions. At high p_T, R_AA is similar for 62.4 and 200 GeV at all centralities. Perturbative-quantum-chromodynamics calculations that describe R_AA well at 200 GeV, fail to describe the 39 GeV data, raising the possibility that the relative importance of initial-state effects and soft processes increases at lower energies. A conclusion that the region where hard processes are dominant is reached only at higher p_T, is also supported by the x_T dependence of the x_T-scaling power-law exponent.

preprint2020arXiv

Measurement of $J/ψ$ at forward and backward rapidity in $p$+$p$, $p$$+A$l, $p$$+A$u, and $^3$He+Au collisions at $\sqrt{s_{_{NN}}}=200~{\rm GeV}$

Charmonium is a valuable probe in heavy-ion collisions to study the properties of the quark gluon plasma, and is also an interesting probe in small collision systems to study cold nuclear matter effects, which are also present in large collision systems. With the recent observations of collective behavior of produced particles in small system collisions, measurements of the modification of charmonium in small systems have become increasingly relevant. We present the results of $J/ψ$ measurements at forward and backward rapidity in various small collision systems, $p$$+$$p$, $p$$+$Al, $p$$+$Au and $^3$He$+$Au, at $\sqrt{s_{_{NN}}}$=200 GeV. The results are presented in the form of the observable $R_{AB}$, the nuclear modification factor, a measure of the ratio of the $J/ψ$ invariant yield compared to the scaled yield in $p$$+$$p$ collisions. We examine the rapidity, transverse momentum, and collision centrality dependence of nuclear effects on $J/ψ$ production with different projectile sizes $p$ and $^3$He, and different target sizes Al and Au. The modification is found to be strongly dependent on the target size, but to be very similar for $p$$+$Au and $^{3}$He$+$Au. However, for 0%--20% central collisions at backward rapidity, the modification for $^{3}$He$+$Au is found to be smaller than that for $p$$+$Au, with a mean fit to the ratio of $0.89\pm0.03$(stat)${\pm}0.08$(syst), possibly indicating final state effects due to the larger projectile size.

preprint2020arXiv

Measurement of charged pion double spin asymmetries at midrapidity in longitudinally polarized $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the longitudinal double spin asymmetries, $A_{LL}$, for charged pions at midrapidity ($|η|<0.35$) in longitudinally polarized $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV. These measurements are sensitive to the gluon spin contribution to the total spin of the proton in the parton momentum fraction $x$ range between 0.04 and 0.09. One can infer the sign of the gluon polarization from the ordering of pion asymmetries with charge alone. The asymmetries are found to be consistent with global quantum-chromodynamics fits of deep-inelastic scattering and data at $\sqrt{s}=200$ GeV, which show a nonzero positive contribution of gluon spin to the proton spin.

preprint2020arXiv

Nuclear-modification factor of charged hadrons at forward and backward rapidity in $p$$+$Al and $p$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV

The PHENIX experiment has studied nuclear effects in $p$$+$Al and $p$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV on charged hadron production at forward rapidity ($1.4<η<2.4$, $p$-going direction) and backward rapidity ($-2.2<η<-1.2$, $A$-going direction). Such effects are quantified by measuring nuclear modification factors as a function of transverse momentum and pseudorapidity in various collision multiplicity selections. In central $p$$+$Al and $p$$+$Au collisions, a suppression (enhancement) is observed at forward (backward) rapidity compared to the binary scaled yields in $p$+$p$ collisions. The magnitude of enhancement at backward rapidity is larger in $p$$+$Au collisions than in $p$$+$Al collisions, which have a smaller number of participating nucleons. However, the results at forward rapidity show a similar suppression within uncertainties. The results in the integrated centrality are compared with calculations using nuclear parton distribution functions, which show a reasonable agreement at the forward rapidity but fail to describe the backward rapidity enhancement.

preprint2020arXiv

Production of $π^0$ and $η$ mesons in U$+$U collisions at $\sqrt{s_{_{NN}}}=192$ GeV

The PHENIX experiment at the Relativistic Heavy Ion Collider measured $π^0$ and $η$ mesons at midrapidity in U$+$U collisions at $\sqrt{s_{_{NN}}}=192$ GeV in a wide transverse momentum range. Measurements were performed in the $π^0(η)\rightarrowγγ$ decay modes. A strong suppression of $π^0$ and $η$ meson production at high transverse momentum was observed in central U$+$U collisions relative to binary scaled $p$$+$$p$ results. Yields of $π^0$ and $η$ mesons measured in U$+$U collisions show similar suppression pattern to the ones measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV for similar numbers of participant nucleons. The $η$/$π^0$ ratios do not show dependence on centrality or transverse momentum, and are consistent with previously measured values in hadron-hadron, hadron-nucleus, nucleus-nucleus, and $e^+e^-$ collisions.

preprint2016arXiv

Conformal Invariance of Graphene Sheets

Suspended graphene sheets exhibit correlated random deformations that can be studied under the framework of rough surfaces with a Hurst (roughness) exponent $0.72 \pm 0.01$. Here, we show that, independent of the temperature, the iso-height lines at the percolation threshold have a well-defined fractal dimension and are conformally invariant, sharing the same statistical properties as Schramm-Loewner evolution (SLE$_κ$) curves with $κ=2.24\pm0.07$. Interestingly, iso-height lines of other rough surfaces are not necessarily conformally invariant even if they have the same Hurst exponent, e.g. random Gaussian surfaces. We have found that the distribution of the modulus of the Fourier coefficients plays an important role on this property. Our results not only introduce a new universality class and place the study of suspended graphene membranes within the theory of critical phenomena, but also provide hints on the long-standing question about the origin of conformal invariance in iso-height lines of rough surfaces.

preprint2016arXiv

Crumpling Damaged Graphene

Through molecular mechanics we find that non-covalent interactions modify the fractality of crumpled damaged graphene. Pristine graphene membranes are damaged by adding random vacancies and carbon-hydrogen bonds. Crumpled membranes exhibit a fractal dimension of $ 2.71 \pm 0.02$ when all interactions between carbon atoms are considered, and $2.30 \pm 0.05$ when non-covalent interactions are suppressed. The transition between these two values, obtained by switching on/off the non-covalent interactions of equilibrium configurations, is shown to be reversible and independent on thermalisation. In order to explain this transition, we propose a theoretical model that is compatible with our numerical findings. Finally, we also compare damaged graphene membranes with other crumpled structures, as for instance, polymerised membranes and paper sheets, that share similar scaling properties.

preprint2016arXiv

Energy dissipation in flows through curved spaces

Fluid dynamics in intrinsically curved geometries is encountered in many physical systems in nature, ranging from microscopic bio-membranes all the way up to general relativity at cosmological scales. Despite the diversity of applications, all of these systems share a common feature: the free motion of particles is affected by inertial forces originating from the curvature of the embedding space. Here we reveal a fundamental process underlying fluid dynamics in curved space: the free motion of fluids, in the complete absence of solid walls or obstacles, exhibits loss of energy due exclusively to the intrinsic curvature of space. We find that local sources of curvature generate viscous stresses as a result of the inertial forces. The curvature-induced viscous forces are shown to cause hitherto unnoticed and yet appreciable energy dissipation, which might play a significant role for a variety of physical systems involving fluid dynamics in curved spaces.

preprint2016arXiv

Measurement of higher cumulants of net-charge multiplicity distributions in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=7.7-200$ GeV

We report the measurement of cumulants ($C_n, n=1\ldots4$) of the net-charge distributions measured within pseudorapidity ($|η|<0.35$) in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=7.7-200$ GeV with the PHENIX experiment at the Relativistic Heavy Ion Collider. The ratios of cumulants (e.g. $C_1/C_2$, $C_3/C_1$) of the net-charge distributions, which can be related to volume independent susceptibility ratios, are studied as a function of centrality and energy. These quantities are important to understand the quantum-chromodynamics phase diagram and possible existence of a critical end point. The measured values are very well described by expectation from negative binomial distributions. We do not observe any nonmonotonic behavior in the ratios of the cumulants as a function of collision energy. The measured values of $C_1/C_2 = μ/σ^2$ and $C_3/C_1 = Sσ^3/μ$ can be directly compared to lattice quantum-chromodynamics calculations and thus allow extraction of both the chemical freeze-out temperature and the baryon chemical potential at each center-of-mass energy.

preprint2016arXiv

Measurement of parity-violating spin asymmetries in W$^{\pm}$ production at midrapidity in longitudinally polarized $p$$+$$p$ collisions

We present measurements from the PHENIX experiment of large parity-violating single spin asymmetries of high transverse momentum electrons and positrons from $W^\pm/Z$ decays, produced in longitudinally polarized $p$$+$$p$ collisions at center of mass energies of $\sqrt{s}$=500 and 510~GeV. These asymmetries allow direct access to the anti-quark polarized parton distribution functions due to the parity-violating nature of the $W$-boson coupling to quarks and anti-quarks. The results presented are based on data collected in 2011, 2012, and 2013 with an integrated luminosity of 240 pb$^{-1}$, which exceeds previous PHENIX published results by a factor of more than 27. These high $Q^2$ data provide an important addition to our understanding of anti-quark parton helicity distribution functions.

preprint2016arXiv

Poiseuille flow in curved spaces

We investigate Poiseuille channel flow through intrinsically curved media, equipped with localized metric perturbations. To this end, we study the flux of a fluid driven through the curved channel in dependence of the spatial deformation, characterized by the parameters of the metric perturbations (amplitude, range and density). We find that the flux depends only on a specific combination of parameters, which we identify as the average metric perturbation, and derive a universal flux law for the Poiseuille flow. For the purpose of this study, we have improved and validated our recently developed lattice Boltzmann model in curved space by considerably reducing discrete lattice effects.

preprint2016arXiv

Transverse energy production and charged-particle multiplicity at midrapidity in various systems from $\sqrt{s_{NN}}=7.7$ to 200 GeV

Measurements of midrapidity charged particle multiplicity distributions, $dN_{\rm ch}/dη$, and midrapidity transverse-energy distributions, $dE_T/dη$, are presented for a variety of collision systems and energies. Included are distributions for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$, 130, 62.4, 39, 27, 19.6, 14.5, and 7.7 GeV, Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=200$ and 62.4 GeV, Cu$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, U$+$U collisions at $\sqrt{s_{_{NN}}}=193$ GeV, $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, and $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Centrality-dependent distributions at midrapidity are presented in terms of the number of nucleon participants, $N_{\rm part}$, and the number of constituent quark participants, $N_{q{\rm p}}$. For all $A$$+$$A$ collisions down to $\sqrt{s_{_{NN}}}=7.7$ GeV, it is observed that the midrapidity data are better described by scaling with $N_{q{\rm p}}$ than scaling with $N_{\rm part}$. Also presented are estimates of the Bjorken energy density, $\varepsilon_{\rm BJ}$, and the ratio of $dE_T/dη$ to $dN_{\rm ch}/dη$, the latter of which is seen to be constant as a function of centrality for all systems.

preprint2015arXiv

An Upgrade Proposal from the PHENIX Collaboration

In this document the PHENIX collaboration proposes a major upgrade to the PHENIX detector at the Relativistic Heavy Ion Collider. This upgrade, sPHENIX, enables an extremely rich jet and beauty quarkonia physics program addressing fundamental questions about the nature of the strongly coupled quark-gluon plasma (QGP), discovered experimentally at RHIC to be a perfect fluid. The startling dynamics of the QGP on fluid-like length scales is an emergent property of quantum chromodynamics (QCD), seemingly implicit in the Lagrangian but stubbornly hidden from view. QCD is an asymptotically free theory, but how QCD manifests as a strongly coupled fluid with specific shear viscosity near $T_C$, as low as allowed by the uncertainty principle, is as fundamental an issue as that of how confinement itself arises.

preprint2015arXiv

Cooling Effect of the Richtmyer-Meshkov Instability

We provide numerical evidence that the Richtmyer-Meshkov (RM) instability contributes to the cooling of a relativistic fluid. Due to the presence of jet particles traveling throughout the medium, shock waves are generated in the form of Mach cones. The interaction of multiple shock waves can trigger the RM instability, and we have found that this process leads to a down-cooling of the relativistic fluid. To confirm the cooling effect of the instability, shock tube Richtmyer-Meshkov instability simulations are performed. Additionally, in order to provide an experimental observable of the RM instability resulting from the Mach cone interaction, we measure the two particle correlation function and highlight the effects of the interaction. The simulations have been performed with an improved version of the relativistic lattice Boltzmann model, including general equations of state and external forces.

preprint2015arXiv

Coupled DEM-LBM method for the free-surface simulation of heterogeneous suspensions

The complexity of the interactions between the constituent granular and liquid phases of a suspension requires an adequate treatment of the constituents themselves. A promising way for numerical simulations of such systems is given by hybrid computational frameworks. This is naturally done, when the Lagrangian description of particle dynamics of the granular phase finds a correspondence in the fluid description. In this work we employ extensions of the Lattice-Boltzmann Method for non-Newtonian rheology, free surfaces, and moving boundaries. The models allows for a full coupling of the phases, but in a simplified way. An experimental validation is given by an example of gravity driven flow of a particle suspension.

preprint2015arXiv

Flow through randomly curved manifolds

We have found that the relation between the flow through campylotic (generically curved) media, consisting of randomly located curvature perturbations, and the average Ricci scalar of the system exhibits two distinct functional expressions (hysteresis), depending on whether the typical spatial extent of the curvature perturbation lies above or below the critical value maximizing the overall Ricci curvature. Furthermore, the flow through such systems as a function of the number of curvature perturbations presents a sublinear behavior for large concentrations due to the interference between curvature perturbations that, consequently, produces a less curved space. For the purpose of this study, we have developed and validated a lattice kinetic model capable of describing fluid flow in arbitrarily curved manifolds, which allows to deal with highly complex spaces in a very compact and efficient way.

preprint2015arXiv

Granular front formation in free-surface flow of concentrated suspensions

Granular fronts are a common yet unexplained phenomenon emerging during the gravity driven free-surface flow of concentrated suspensions. They are usually believed to be the result of fluid convection in combination with particle size segregation. However, suspensions composed of uniformly sized particles also develop a granular front. Within a large rotating drum, a stationary recirculating avalanche is generated. The flowing material is a mixture of a visco-plastic fluid obtained from a kaolin-water dispersion, with spherical ceramic particles denser than the fluid. The goal is to mimic the composition of many common granular-fluid materials, like fresh concrete or debris flow. In these materials, granular and fluid phases have the natural tendency to segregate due to particle settling. However, through the shearing caused by the rotation of the drum, a reorganization of the phases is induced, leading to the formation of a granular front. By tuning the material properties and the drum velocity, it is possible to control this phenomenon. The setting is reproduced in a numerical environment, where the fluid is solved by a Lattice-Boltzmann Method, and the particles are explicitly represented using the Discrete Element Method. The simulations confirm the findings of the experiments, and provide insight into the internal mechanisms. Comparing the time-scale of particle settling with the one of particle recirculation, a non-dimensional number is defined, and is found to be effective in predicting the formation of a granular front.

preprint2015arXiv

High-Order Kinetic Relaxation Schemes as High-Accuracy Poisson Solvers

We present a new approach to find accurate solutions to the Poisson equation, as obtained from the steady-state limit of a diffusion equation with strong source terms. For this purpose, we start from Boltzmann's kinetic theory and investigate the influence of higher order terms on the resulting macroscopic equations. By performing an appropriate expansion of the equilibrium distribution, we provide a method to remove the unnecessary terms up to a desired order and show that it is possible to find, with high level of accuracy, the steady-state solution of the diffusion equation for sizeable Knudsen numbers. In order to test our kinetic approach, we discretise the Boltzmann equation and solve the Poisson equation, spending up to six order of magnitude less computational time for a given precision than standard lattice Boltzmann methods.

preprint2015arXiv

Lattice Boltzmann Model for Electronic Structure Simulations

Recently, a new connection between density functional theory and kinetic theory has been proposed. In particular, it was shown that the Kohn-Sham (KS) equations can be reformulated as a macroscopic limit of the steady-state solution of a suitable single-particle kinetic equation. By using a discrete version of this new formalism, the exchange and correlation energies of simple atoms and the geometrical configuration of the methane molecule were calculated accurately. Here, we discuss the main ideas behind the lattice kinetic approach to electronic structure computations, offer some considerations for prospective extensions, and also show additional numerical results, namely the geometrical configuration of the water molecule.

preprint2015arXiv

Lattice Boltzmann Model for Numerical Relativity

In the Bona-Masso formulation, Einstein equations are written as a set of flux conservative first order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for Numerical Relativity. Our model is validated with well-established tests, showing good agreement with analytical solutions. Furthermore, we show that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improves. Finally, in order to show the potential of our approach a linear scaling law for parallelisation with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.

preprint2015arXiv

Lattice Boltzmann model for resistive relativistic magnetohydrodynamics

In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfvén waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to $σ^{-\frac{1}{2}}$, $σ$ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere.

preprint2015arXiv

Search for dark photons from neutral meson decays in $p$$+$$p$ and $d$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV

The standard model (SM) of particle physics is spectacularly successful, yet the measured value of the muon anomalous magnetic moment $(g-2)_μ$ deviates from SM calculations by 3.6$σ$. Several theoretical models attribute this to the existence of a "dark photon," an additional U(1) gauge boson, which is weakly coupled to ordinary photons. The PHENIX experiment at the Relativistic Heavy Ion Collider has searched for a dark photon, $U$, in $π^0,η\rightarrow γe^+e^-$ decays and obtained upper limits of $\mathcal{O}(2\times10^{-6})$ on $U$-$γ$ mixing at 90% CL for the mass range $30<m_U<90$ MeV/$c^2$. Combined with other experimental limits, the remaining region in the $U$-$γ$ mixing parameter space that can explain the $(g-2)_μ$ deviation from its SM value is nearly completely excluded at the 90% confidence level, with only a small region of $29<m_U<32$ MeV/$c^2$ remaining.

preprint2015arXiv

Short-lived lattice quasiparticles for strongly interacting fluids

It is shown that lattice kinetic theory based on short-lived quasiparticles proves very effective in simulating the complex dynamics of strongly interacting fluids (SIF). In particular, it is pointed out that the shear viscosity of lattice fluids is the sum of two contributions, one due to the usual interactions between particles (collision viscosity) and the other due to the interaction with the discrete lattice (propagation viscosity). Since the latter is {\it negative}, the sum may turn out to be orders of magnitude smaller than each of the two contributions separately, thus providing a mechanism to access SIF regimes at ordinary values of the collisional viscosity. This concept, as applied to quantum superfluids in one-dimensional optical lattices, is shown to reproduce shear viscosities consistent with the AdS-CFT holographic bound on the viscosity/entropy ratio. This shows that lattice kinetic theory continues to hold for strongly coupled hydrodynamic regimes where continuum kinetic theory may no longer be applicable.

preprint2014arXiv

A Gauge Condition for Studying the Origin of Intrinsic Magnetospheres

We propose an analytical model based on the solution of the magnetohydrodynamics (MHD) equations for studying the origin of intrinsic magnetospheres. For this purpose, we reveal a new gauge condition for the electromagnetic vector potential, which eases the solution of such complex system of non-linear equations. Using this model, we analyse the deformation of the terrestrial magnetic field due to the presence of the solar wind. By comparing with experimental observations, we have found that the geometrical configuration of the magnetosphere, before, and after the magnetic field of the Earth started to deviate the solar wind, has not changed notably, and that the solar wind should have had a finite conductivity. This model can also be used to perform linear stability analysis of fluid and magnetic instabilities.

preprint2014arXiv

Dean Instability in Double-Curved Channels

We study the Dean instability in curved channels using the lattice Boltzmann model for generalized metrics. For this purpose, we first improve and validate the method by measuring the critical Dean number at the transition from laminar to vortex flow for a streamwise curved rectangular channel, obtaining very good agreement with the literature values. Taking advantage of the easy implementation of arbitrary metrics within our model, we study the fluid flow through a double-curved channel, using ellipsoidal coordinates, and study the transition to vortex flow in dependence of the two perpendicular curvature radii of the channel. We observe not only transitions to 2-cell vortex flow, but also to 4-cell and even 6-cell vortex flow, and we find that the critical Dean number at the transition to 2-cell vortex flow exhibits a minimum when the two curvature radii are approximately equal.

preprint2014arXiv

Modelling the influence of photospheric turbulence on solar flare statistics

Solar flares stem from the reconnection of twisted magnetic field lines in the solar photosphere. The energy and waiting time distributions of these events follow complex patterns that have been carefully considered in the past and that bear some resemblance with earthquakes and stockmarkets. Here we explore in detail the tangling motion of interacting flux tubes anchored in the plasma and the energy ejections resulting when they recombine. The mechanism for energy accumulation and release in the flow is reminiscent of self-organized criticality. From this model we suggest the origin for two important and widely studied properties of solar flare statistics, including the time-energy correlations. We first propose that the scale-free energy distribution of solar flares is largely due to the twist exerted by the vorticity of the turbulent photosphere. Second, the long-range temporal and time-energy correlations appear to arise from the tube-tube interactions. The agreement with satellite measurements is encouraging.

preprint2014arXiv

Nuclear matter effects on $J/ψ$ production in asymmetric Cu+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV

We report on $J/ψ$ production from asymmetric Cu+Au heavy-ion collisions at $\sqrt{s_{_{NN}}}$=200 GeV at the Relativistic Heavy Ion Collider at both forward (Cu-going direction) and backward (Au-going direction) rapidities. The nuclear modification of $J/ψ$ yields in Cu$+$Au collisions in the Au-going direction is found to be comparable to that in Au$+$Au collisions when plotted as a function of the number of participating nucleons. In the Cu-going direction, $J/ψ$ production shows a stronger suppression. This difference is comparable in magnitude and has the same sign as the difference expected from shadowing effects due to stronger low-$x$ gluon suppression in the larger Au nucleus. The relative suppression is opposite to that expected from hot nuclear matter dissociation, since a higher energy density is expected in the Au-going direction.

preprint2013arXiv

Formal analogy between the Dirac equation in its Majorana form and the discrete-velocity version of the Boltzmann kinetic equation

We point out a formal analogy between the Dirac equation in Majorana form and the discrete-velocity version of the Boltzmann kinetic equation. By a systematic analysis based on the theory of operator splitting, this analogy is shown to turn into a concrete and efficient computational method, providing a unified treatment of relativistic and non-relativistic quantum mechanics. This might have potentially far-reaching implications for both classical and quantum computing, because it shows that, by splitting time along the three spatial directions, quantum information (Dirac-Majorana wavefunction) propagates in space-time as a classical statistical process (Boltzmann distribution).

preprint2013arXiv

Gaussian Quadrature and Lattice Discretization of the Fermi-Dirac Distribution for Graphene

We construct a lattice kinetic scheme to study electronic flow in graphene. For this purpose, we first derive a basis of orthogonal polynomials, using as weight function the ultrarelativistic Fermi-Dirac distribution at rest. Later, we use these polynomials to expand the respective distribution in a moving frame, for both cases, undoped and doped graphene. In order to discretize the Boltzmann equation and make feasible the numerical implementation, we reduce the number of discrete points in momentum space to 18 by applying a Gaussian quadrature, finding that the family of representative wave (2+1)-vectors, that satisfies the quadrature, reconstructs a honeycomb lattice. The procedure and discrete model are validated by solving the Riemann problem, finding excellent agreement with other numerical models. In addition, we have extended the Riemann problem to the case of different dopings, finding that by increasing the chemical potential, the electronic fluid behaves as if it increases its effective viscosity.

preprint2013arXiv

Hydrodynamic Model for Conductivity in Graphene

Based on the recently developed picture of an electronic ideal relativistic fluid at the Dirac point, we present an analytical model for the conductivity in graphene that is able to describe the linear dependence on the carrier density and the existence of a minimum conductivity. The model treats impurities as submerged rigid obstacles, forming a disordered medium through which graphene electrons flow, in close analogy with classical fluid dynamics. To describe the minimum conductivity, we take into account the additional carrier density induced by the impurities in the sample. The model, which predicts the conductivity as a function of the impurity fraction of the sample, is supported by extensive simulations for different values of ${\cal E}$, the dimensionless strength of the electric field, and provides excellent agreement with experimental data.

preprint2013arXiv

Kinetic Formulation of the Kohn-Sham Equations for ab initio Electronic Structure Calculations

We introduce a new approach to density functional theory based on kinetic theory, showing that the Kohn-Sham equations can be derived as a macroscopic limit of a suitable Boltzmann kinetic equation in the limit of small mean free path versus the typical scale of density gradients (Chapman-Enskog expansion). To derive the approach, we first write the Schrödinger equation as a special case of a Boltzmann equation for a gas of quasi-particles, with the potential playing the role of an external source that generates and destroys particles, so as to drive the system towards the ground state. The ions are treated as classical particles, using the Born-Oppenheimer dynamics, or by imposing concurrent evolution with the electronic orbitals. In order to provide quantitative support to our approach, we implement a discrete (lattice) model and compute, the exchange and correlation energies of simple atoms, and the geometrical configuration of the methane molecule. Excellent agreement with values in the literature is found.

preprint2013arXiv

Lattice Boltzmann model for ultra-relativistic flows

We develop a relativistic lattice Boltzmann model capable of describing relativistic fluid dynamics at ultra-high velocities, with Lorentz factors up to $γ\sim 10$. To this purpose, we first build a new lattice kinetic scheme by expanding the Maxwell-Jüttner distribution function in an orthogonal basis of polynomials and applying an appropriate quadrature, providing the discrete versions of the relativistic Boltzmann equation and the equilibrium distribution. To achieve ultra-high velocities, we include a flux limiter scheme, and introduce the bulk viscosity by a suitable extension of the discrete relativistic Boltzmann equation. The model is validated by performing simulations of shock waves in viscous quark-gluon plasmas and comparing with existing models, finding very good agreement. To the best of our knowledge, we for the first time successfully simulate viscous shock waves in the highly relativistic regime. Moreover, we show that our model can also be used for near-inviscid flows even at very high velocities. Finally, as an astrophysical application, we simulate a relativistic shock wave, generated by, say, a supernova explosion, colliding with a massive interstellar cloud, e.g. molecular gas.

preprint2013arXiv

Relativistic effects on the Richtmyer-Meshkov instability

Theoretical and numerical analysis of the relativistic effects on the Richtmyer-Meshkov (RM) instability reveals new and potentially very useful effects. We find that, in contrast with the non- relativistic case, the growth rate of the RM instability depends strongly on the equation of state of the fluid, opening up the possibility to infer equations of state from experimental observations of the RM instability. As opposed to the non-relativistic case, we also discover that, above a critical value of the fluid velocity, the growth rate of the instability counter-intuitively decreases due to the Lorentz's factor, and vanishes in the ultrarelativistic limit, as the speed of the particles approaches the speed of light. Both effects might prove very useful for leading-edge applications, such as the study of the equation of state of quark-gluon matter, and the design of fast ignition inertial confinement fusion (ICF) schemes. We perform a linear stability analysis to characterize the instability, for an arbitrary equation of state, and implement numerical simulations to study the instability in the non-linear regime, using the equation of state of an ideal gas. Furthermore, based on the numerical results, we propose a general expression that characterizes the long term evolution of the instability.

preprint2013arXiv

Relativistic Lattice Boltzmann Model with Improved Dissipation

We develop a relativistic lattice Boltzmann (LB) model, providing a more accurate description of dissipative phenomena in relativistic hydrodynamics than previously available with existing LB schemes. The procedure applies to the ultra-relativistic regime, in which the kinetic energy (temperature) far exceeds the rest mass energy, although the extension to massive particles and/or low temperatures is conceptually straightforward. In order to improve the description of dissipative effects, the Maxwell-Juettner distribution is expanded in a basis of orthonormal polynomials, so as to correctly recover the third order moment of the distribution function. In addition, a time dilatation is also applied, in order to preserve the compatibility of the scheme with a cartesian cubic lattice. To the purpose of comparing the present LB model with previous ones, the time transformation is also applied to a lattice model which recovers terms up to second order, namely up to energy-momentum tensor. The approach is validated through quantitative comparison between the second and third order schemes with BAMPS (the solution of the full relativistic Boltzmann equation), for moderately high viscosity and velocities, and also with previous LB models in the literature. Excellent agreement with BAMPS and more accurate results than previous relativistic lattice Boltzmann models are reported.

preprint2013arXiv

Ultrarelativistic Transport Coefficients in Two Dimensions

We compute the shear and bulk viscosities, as well as the thermal conductivity of an ultrarelativistic fluid obeying the relativistic Boltzmann equation in 2+1 space-time dimensions. The relativistic Boltzmann equation is taken in the single relaxation time approximation, based on two approaches, the first, due to Marle and using the Eckart decomposition, and the second, proposed by Anderson and Witting and using the Landau-Lifshitz decomposition. In both cases, the local equilibrium is given by a Maxwell-Juettner distribution. It is shown that, apart from slightly different numerical prefactors, the two models lead to a different dependence of the transport coefficients on the fluid temperature, quadratic and linear, for the case of Marle and Anderson-Witting, respectively. However, by modifying the Marle model according to the prescriptions given in Ref.[1], it is found that the temperature dependence becomes the same as for the Anderson-Witting model.

preprint2012arXiv

J/psi suppression at forward rapidity in Au+Au collisions at sqrt(s_NN)=39 and 62.4 GeV

We present measurements of the J/psi invariant yields in sqrt(s_NN)=39 and 62.4 GeV Au+Au collisions at forward rapidity (1.2<|y|<2.2). Invariant yields are presented as a function of both collision centrality and transverse momentum. Nuclear modifications are obtained for central relative to peripheral Au+Au collisions (R_CP) and for various centrality selections in Au+Au relative to scaled p+p cross sections obtained from other measurements (R_AA). The observed suppression patterns at 39 and 62.4 GeV are quite similar to those previously measured at 200 GeV. This similar suppression presents a challenge to theoretical models that contain various competing mechanisms with different energy dependencies, some of which cause suppression and others enhancement.

preprint2012arXiv

Klein Tunneling in the presence of random impurities

In this paper, we study Klein tunneling in random media. To this purpose, we simulate the propagation of a relativistic Gaussian wavepacket through a graphene sample with randomly distributed potential barriers (impurities). The simulations, based on a relativistic quantum lattice Boltzmann method, permit to compute the transmission coefficient across the sample, thereby providing an estimate for the conductivity as a function of impurity concentration and strength of the potentials. It is found that the conductivity loss due to impurities is significantly higher for wave-packets of massive particles, as compared to massless ones. A general expression for the loss of conductivity as a function of the impurity percentage is presented and successfully compared with the Kozeny-Carman law for disordered media in classical fluid dynamics.

preprint2012arXiv

Medium modification of jet fragmentation in Au+Au collisions at sqrt(s_NN)=200 GeV measured in direct photon-hadron correlations

The jet fragmentation function is measured with direct photon-hadron correlations in p+p and Au+Au collisions at sqrt(s_NN)=200 GeV. The p_T of the photon is an excellent approximation to the initial p_T of the jet and the ratio z_T=p_T^h/p_T^γis used as a proxy for the jet fragmentation function. A statistical subtraction is used to extract the direct photon-hadron yields in Au+Au collisions while a photon isolation cut is applied in p+p. I_ AA, the ratio of jet fragment yield in Au+Au to that in p+p, indicates modification of the jet fragmentation function. Suppression, most likely due to energy loss in the medium, is seen at high z_T. The fragment yield at low z_T is enhanced at large angles. Such a trend is expected from redistribution of the lost energy into increased production of low-momentum particles.

preprint2012arXiv

Pre-Turbulent Regimes in Graphene Flows

We provide numerical evidence that electronic pre-turbulent phenomena in graphene could be observed, under current experimental conditions, through detectable current fluctuations, echoing the detachment of vortices past localized micron-sized impurities. Vortex generation, due to micron sized constriction, is also explored with special focus on the effects of relativistic corrections to the normal Navier-Stokes equations. These corrections are found to cause a delay in the stability breakout of the fluid as well as a small shift in the vortex shedding frequency. Finally, a relation between the Strouhal number, a dimensionless measure of the vortex shedding frequency, and the Reynolds number is provided under conditions of interest for future experiments.

preprint2012arXiv

Single to Double Hump Transition in the Equilibrium Distribution Function of Relativistic Particles

We unveil a transition from single peaked to bimodal velocity distribution in a relativistic fluid under increasing temperature, in contrast with a non-relativistic gas, where only a monotonic broadening of the bell-shaped distribution is observed. Such transition results from the interplay between the raise in thermal energy and the constraint of maximum velocity imposed by the speed of light. We study the Bose-Einstein, the Fermi-Dirac, and the Maxwell-Jüttner distributions, all exhibiting the same qualitative behavior. We characterize the nature of the transition in the framework of critical phenomena and show that it is either continuous or discontinuous, depending on the group velocity. We analyze the transition in one, two, and three dimensions, with special emphasis on two-dimensions, for which a possible experiment in graphene, based on the measurement of the Johnson-Nyquist noise, is proposed.

preprint2012arXiv

sPHENIX: An Upgrade Concept from the PHENIX Collaboration

The PHENIX collaboration presents a concept for a major upgrade to the PHENIX detector at the Relativistic Heavy Ion Collider (RHIC). This upgrade, referred to as sPHENIX, brings exciting new capability to the RHIC program by opening new and important channels for experimental investigation and utilizing fully the luminosity of the recently upgraded RHIC facility. sPHENIX enables a compelling jet physics program that will address fundamental questions about the nature of the strongly coupled quark-gluon plasma discovered experimentally at RHIC to be a perfect fluid. The upgrade concept addresses specific questions whose answers are necessary to advance our understanding of the quark-gluon plasma: (1) How to reconcile the observed strongly coupled quark-gluon plasma with the asymptotically free theory of quarks and gluons? (2) What are the dynamical changes to the quark-gluon plasma in terms of quasiparticles and excitations as a function of temperature? (3) How sharp is the transition of the quark-gluon plasma from the most strongly coupled regime near Tc to a weakly coupled system of partons known to emerge at asymptotically high temperatures? In three Appendices, we detail the additional physics capabilities gained through further upgrades: (A) two midrapidity detector additions, (B) a forward rapidity upgrade, and (C) an evolution to an ePHENIX detector suitable for a future Electron Ion Collider at RHIC.

preprint2012arXiv

Taylor-Couette Instability in General Manifolds: A Lattice Kinetic Approach

We present a new lattice kinetic method to simulate fluid dynamics in curvilinear geometries. A suitable discrete Boltzmann equation is solved in contravariant coordinates, and the equilibrium distribution function is obtained by a Hermite polynomials expansion of the Maxwell-Boltzmann distribution, expressed in terms of the contravariant coordinates and the metric tensor. To validate the model, we calculate the critical Reynolds number for the onset of the Taylor-Couette instability between two concentric cylinders, obtaining excellent agreement with the theory. In order to extend this study to more general geometries, we also calculate the critical Reynolds number for the case of two concentric spheres, finding good agreement with experimental data. In the case of two concentric tori, we have found that the critical Reynolds is about 10% larger than the respective value for the two concentric cylinders.

preprint2011arXiv

A fully relativistic lattice Boltzmann algorithm

Starting from the Maxwell-Juettner equilibrium distribution, we develop a relativistic lattice Boltzmann (LB) algorithm capable of handling ultrarelativistic systems with flat, but expanding, spacetimes. The algorithm is validated through simulations of quark-gluon plasma, yielding excellent agreement with hydrodynamic simulations. The present scheme opens the possibility of transferring the recognized computational advantages of lattice kinetic theory to the context of both weakly and ultra-relativistic systems.

preprint2011arXiv

On the relativistic Lattice Boltzmann method for quark-gluon plasma simulations

In this paper, we investigate the recently developed lattice Boltzmann model for relativistic hydrodynamics. To this purpose, we perform simulations of shock waves in quark-gluon plasma in the low and high viscosities regime, using three different computational models, the relativistic lattice Boltzmann (RLB), the Boltzmann Approach Multi-Parton Scattering (BAMPS), and the viscous sharp and smooth transport algorithm (vSHASTA). From the results, we conclude that the RLB model departs from BAMPS in the case of high speeds and high temperature(viscosities), the departure being due to the fact that the RLB is based on a quadratic approximation of the Maxwell-Jüttner distribution, which is only valid for sufficiently low temperature and velocity. Furthermore, we have investigated the influence of the lattice speed on the results, and shown that inclusion of quadratic terms in the equilibrium distribution improves the stability of the method within its domain of applicability. Finally, we assess the viability of the RLB model in the various parameter regimes relevant to ultra-relativistic fluid dynamics.

preprint2010arXiv

Adhesive Penetration in Beech Wood Part I: Experiments

A study with synchrotron radiation X-ray tomographic microscopy (SRXTM) of PUR, PVAC, and UF adhesive bond lines in beech wood, bonded under various growth ring angles is presented. After determining the hardening characteristics of the adhesives, we evaluate the bond line morphologies, and the adhesive penetration into the wood structure. We find distinct bond line imperfections for the different adhesive systems. To describe the adhesive distribution inside the bond line we propose the saturation of the pore space instead of the commonly used maximum penetration depth. The results are the basis for a penetration model of hardening fluids into hardwood (part II).

preprint2010arXiv

Adhesive penetration in Beech wood Part II: Penetration Model

We propose an analytical model to predict the adhesives penetration into hard wood. Penetration of hard wood is dominated by the vessel network which prohibits porous medium approximations. Our model considers two scales: a one dimensional capillary fluid transport of a hardening adhesive through a single, straight vessel with diffusion of solvent through its walls and a mesoscopic scale based on topological characteristics of the vessel network, where results from the single vessel scale are mapped onto a periodic network. Given an initial amount of adhesive and applied bonding pressure, we calculate the portion of the filled structure. The model is applied to beech wood samples joined with three different types of adhesive (PUR, UF, PVAC) under various growth ring angles. We evaluate adhesive properties and bond line morphologies described in part I of this work. The model contains one free parameter that can be adjusted in order to fit the experimental data.

preprint2010arXiv

Cross Section and Parity Violating Spin Asymmetries of W^+/- Boson Production in Polarized p+p Collisions at sqrt(s)=500 GeV

Large parity violating longitudinal single spin asymmetries A^{e^-}_L= -0.86^{+0.14}_{-0.30} and A^{e^+}_L= 0.88^{+0.12}_{-0.71} are observed for inclusive high transverse momentum electrons and positrons in polarized pp collisions at a center of mass energy of \sqrt{s}=500\ GeV with the PHENIX detector at RHIC. These e^{+/-} come mainly from the decay of W^{+/-} and Z^0 bosons, and the asymmetries directly demonstrate parity violation in the couplings of the W^{\pm} to the light quarks. The observed electron and positron yields were used to estimate W^\pm boson production cross sections equal to σ(pp \to W^+ X) \times BR(W^ \to ν_e)= 144.1+/-21.2(stat)^{+3.4}_{-10.3}(syst) +/- 15%(norm) pb, and σ(pp \to W^{-}X) \times BR(W^\to e^-\bar{ν_e}) = 31.7+/-12.1(stat)^{+10.1}_{-8.2}(syst)+/-15%(norm) pb.

preprint2010arXiv

Derivation of the Lattice Boltzmann Model for Relativistic Hydrodynamics

A detailed derivation of the Lattice Boltzmann (LB) scheme for relativistic fluids recently proposed in Ref. [1], is presented. The method is numerically validated and applied to the case of two quite different relativistic fluid dynamic problems, namely shock-wave propagation in quark-gluon plasmas and the impact of a supernova blast-wave on massive interstellar clouds. Close to second order convergence with the grid resolution, as well as linear dependence of computational time on the number of grid points and time-steps, are reported.

preprint2010arXiv

Three Dimensional Lattice-Boltzmann Model for Electrodynamics

In this paper we introduce a novel 3D Lattice-Boltzmann model that recovers in the continuous limit the Maxwell equations in materials. In order to build conservation equations with antisymmetric tensors, like the Faraday law, the model assigns four auxiliary vectors to each velocity vector. These auxiliary vectors, when combined with the distribution functions, give the electromagnetic fields. The evolution is driven by the usual BGK collision rule, but with a different form for the equilibrium distribution functions. This LBGK model allows us to consider for both dielectrics and conductors with realistic parameters, and therefore it is adequate to simulate the most diverse electromagnetic problems, like the propagation of electromagnetic waves (both in dielectric media and in waveguides), the skin effect, the radiation pattern of a small dipole antenna and the natural frequencies of a resonant cavity, all with 2% accuracy. Actually, it shows to be one order of magnitude faster than the original FDTD formulation by Yee to reach the same accuracy. It is, therefore, a valuable alternative to simulate electromagnetic fields and opens lattice Boltzmann for a broad spectrum of new applications in electrodynamics. In this paper we develop a 3D Lattice-Boltzmann model that recovers in the continuous limit the Maxwell equations for macroscopic mediums. The model can sucessfully reproduces the propagation of the electromagnetic waves in dielectric mediums and waveguide, the skin effect, the electrical dipole radiation and the electromagnetic response of a resonant cavity.

preprint2009arXiv

Lattice Boltzmann scheme for relativistic fluids

A Lattice Boltzmann formulation for relativistic fluids is presented and numerically verified through quantitative comparison with recent hydrodynamic simulations of relativistic shock-wave propagation in viscous quark-gluon plasmas. This formulation opens up the possibility of exporting the main advantages of Lattice Boltzmann methods to the relativistic context, which seems particularly useful for the simulation of relativistic fluids in complicated geometries.

preprint2009arXiv

Simulation of Flow of Mixtures Through Anisotropic Porous Media using a Lattice Boltzmann Model

We propose a description for transient penetration simulations of miscible and immiscible fluid mixtures into anisotropic porous media, using the lattice Boltzmann (LB) method. Our model incorporates hydrodynamic flow, diffusion, surface tension, and the possibility for global and local viscosity variations to consider various types of hardening fluids. The miscible mixture consists of two fluids, one governed by the hydrodynamic equations and one by diffusion equations. We validate our model on standard problems like Poiseuille flow, the collision of a drop with an impermeable, hydrophobic interface and the deformation of the fluid due to surface tension forces. To demonstrate the applicability to complex geometries, we simulate the invasion process of mixtures into wood spruce samples.

preprint2006arXiv

3D Lattice-Boltzmann Model for Magnetic Reconnection

In this paper we develop a 3D Lattice-Boltzmann model that recovers in the continuous limit the two-fluids theory for plasmas, and consecuently includes the generalizated Ohm's law. The model reproduces the magnetic reconnection process just by given the right initial equlibrium conditions in the magnetotail, without any assumption on the resistivity in the diffusive region. In this model, the plasma is handled like two fluids with an interaction term, each one with distribution functions associated to a cubic lattice with 19 velocities (D3Q19). The electromagnetic fields are considered like a third fluid with an external force on a cubic lattice with 13 velocities (D3Q13). The model can simulate either viscous fluids in the incompressible limit or non-viscous compressible fluids, and sucessfully reproduces both the Hartmann flow and the magnetic reconnection in the magnetotail. The reconnection rate obtained with this model is R=0.109, which is in excellent agreement with the observations.