We measure triangular flow relative to the reaction plane at 3 GeV center-of-mass energy in Au+Au collisions at the BNL Relativistic Heavy Ion Collider. A significant $v_3$ signal for protons is observed, which increases for higher rapidity, higher transverse momentum, and more peripheral collisions. The triangular flow is essentially rapidity-odd with a slope at mid-rapidity, $dv_3/dy|_{(y=0)}$, opposite in sign compared to the slope for directed flow. No significant $v_3$ signal is observed for charged pions and kaons. Comparisons with models suggest that a mean field potential is required to describe these results, and that the triangular shape of the participant nucleons is the result of stopping and nuclear geometry.

A comprehensive survey of event-by-event fluctuations of charged hadron multiplicity in relativistic heavy ions is presented. The survey covers Au+Au collisions at sqrt(s_NN) = 62.4 and 200 GeV, and Cu+Cu collisions sqrt(s_NN) = 22.5, 62.4, and 200 GeV. Fluctuations are measured as a function of collision centrality, transverse momentum range, and charge sign. After correcting for non-dynamical fluctuations due to fluctuations in the collision geometry within a centrality bin, the remaining dynamical fluctuations expressed as the variance normalized by the mean tend to decrease with increasing centrality. The dynamical fluctuations are consistent with or below the expectation from a superposition of participant nucleon-nucleon collisions based upon p+p data, indicating that this dataset does not exhibit evidence of critical behavior in terms of the compressibility of the system. An analysis of Negative Binomial Distribution fits to the multiplicity distributions demonstrates that the heavy ion data exhibit weak clustering properties.

The PHENIX experiment at RHIC has measured transverse energy and charged particle multiplicity at mid-rapidity in Au+Au collisions at sqrt(s_NN) = 19.6, 130 and 200 GeV as a function of centrality. The presented results are compared to measurements from other RHIC experiments, and experiments at lower energies. The sqrt(s_NN) dependence of dE_T/deta and dN_ch/deta per pair of participants is consistent with logarithmic scaling for the most central events. The centrality dependence of dE_T/deta and dN_ch/deta is similar at all measured incident energies. At RHIC energies the ratio of transverse energy per charged particle was found independent of centrality and growing slowly with sqrt(s_NN). A survey of comparisons between the data and available theoretical models is also presented.

We report here the first observation of directed flow ($v_1$) of the hypernuclei $^3_{\Lambda}$H and $^4_{\Lambda}$H in mid-central Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV at RHIC. These data are taken as part of the beam energy scan program carried out by the STAR experiment. From 165 $\times$ 10$^{6}$ events in 5%-40% centrality, about 8400 $^3_{\Lambda}$H and 5200 $^4_{\Lambda}$H candidates are reconstructed through two- and three-body decay channels. We observe that these hypernuclei exhibit significant directed flow. Comparing to that of light nuclei, it is found that the midrapidity $v_1$ slopes of $^3_{\Lambda}$H and $^4_{\Lambda}$H follow baryon number scaling, implying that the coalescence is the dominant mechanism for these hypernuclei production in such collisions.

We report a measurement of cumulants and correlation functions of event-by-event proton multiplicity distributions from fixed-target Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV measured by the STAR experiment. Protons are identified within the rapidity ($y$) and transverse momentum ($p_{\rm T}$) region $-0.9 < y<0$ and $0.4 < p_{\rm T} <2.0 $ GeV/$c$ in the center-of-mass frame. A systematic analysis of the proton cumulants and correlation functions up to sixth-order as well as the corresponding ratios as a function of the collision centrality, $p_{\rm T}$, and $y$ are presented. The effect of pileup and initial volume fluctuations on these observables and the respective corrections are discussed in detail. The results are compared to calculations from the hadronic transport UrQMD model as well as a hydrodynamic model. In the most central 5% collisions, the value of proton cumulant ratio $C_4/C_2$ is negative, drastically different from the values observed in Au+Au collisions at higher energies. Compared to model calculations including Lattice QCD, a hadronic transport model, and a hydrodynamic model, the strong suppression in the ratio of $C_4/C_2$ at 3 GeV Au+Au collisions indicates an energy regime dominated by hadronic interactions.

The Beam Energy Scan (BES) program at the Relativistic Heavy Ion Collider (RHIC) was extended to energies below $\sqrt{\textit{s}_{NN}}$ = 7.7 GeV in 2015 by successful implementation of the fixed-target mode of operation in the STAR (Solenoidal Track At RHIC) experiment. In the fixed-target mode, ions circulate in one ring of the collider and interact with a stationary target at the entrance of the STAR Time Projection Chamber. The first results for Au+Au collisions at $\sqrt{\textit{s}_{NN}}$ = 4.5 GeV are presented, including directed and elliptic flow of identified hadrons, and radii from pion femtoscopy. The proton flow and pion femtoscopy results agree quantitatively with earlier measurements by Alternating Gradient Synchrotron experiments at similar energies. This validates running the STAR experiment in the fixed-target configuration. Pion directed and elliptic flow are presented for the first time at this beam energy. Pion and proton elliptic flow show behavior which hints at constituent quark scaling, but large error bars preclude reliable conclusions. The ongoing second phase of BES (BES-II) will provide fixed-target data sets with 100 times more events at each of several energies down to $\sqrt{\textit{s}_{NN}}$ = 3.0 GeV.

Measurements of mass and $\Lambda$ binding energy of $\rm ^4_{\Lambda}H$ and $\rm ^4_{\Lambda}He$ in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}=3$ GeV are presented, with an aim to address the charge symmetry breaking (CSB) problem in hypernuclei systems with atomic number A = 4. The $\Lambda$ binding energies are measured to be $\rm 2.22\pm0.06(stat.) \pm0.14(syst.)$ MeV and $\rm 2.38\pm0.13(stat.) \pm0.12(syst.)$ MeV for $\rm ^4_{\Lambda}H$ and $\rm ^4_{\Lambda}He$, respectively. The measured $\Lambda$ binding-energy difference is $\rm 0.16\pm0.14(stat.)\pm0.10(syst.)$ MeV for ground states. Combined with the $\gamma$-ray transition energies, the binding-energy difference for excited states is $\rm -0.16\pm0.14(stat.)\pm0.10(syst.)$ MeV, which is negative and comparable to the value of the ground states within uncertainties. These new measurements on the $\Lambda$ binding-energy difference in A = 4 hypernuclei systems are consistent with the theoretical calculations that result in $\rm \Delta B_{\Lambda}^4(1_{exc}^{+})\approx -\Delta B_{\Lambda}^4(0_{g.s.}^{+})<0$ and present a new method for the study of CSB effect using relativistic heavy-ion collisions.

We report the first multi-differential measurements of strange hadrons of $K^{-}$, $\phi$ and $\Xi^{-}$ yields as well as the ratios of $\phi/K^-$ and $\phi/\Xi^-$ in Au+Au collisions at ${\sqrt{s_{\rm NN}} = \rm{3\,GeV}}$ with the STAR experiment fixed target configuration at RHIC. The $\phi$ mesons and $\Xi^{-}$ hyperons are measured through hadronic decay channels, $\phi\rightarrow K^+K^-$ and $\Xi^-\rightarrow \Lambda\pi^-$. Collision centrality and rapidity dependence of the transverse momentum spectra for these strange hadrons are presented. The $4\pi$ yields and ratios are compared to thermal model and hadronic transport model predictions. At this collision energy, thermal model with grand canonical ensemble (GCE) under-predicts the $\phi/K^-$ and $\phi/\Xi^-$ ratios while the result of canonical ensemble (CE) calculations reproduce $\phi/K^-$, with the correlation length $r_c \sim 2.7$ fm, and $\phi/\Xi^-$, $r_c \sim 4.2$ fm, for the 0-10% central collisions. Hadronic transport models including high mass resonance decays could also describe the ratios. While thermal calculations with GCE work well for strangeness production in high energy collisions, the change to CE at $\rm{3\,GeV}$ implies a rather different medium property at high baryon density.

We report cumulants of the proton multiplicity distribution from dedicated fixed-target Au+Au collisions at 3.0 GeV, measured by the STAR experiment in the kinematic acceptance of rapidity ($y$) and transverse momentum ($p_{\rm T}$) within $-0.5 < y<0$ and $0.4 < p_{\rm T} <2.0 $ GeV/$c$. In the most central 0--5% collisions, a proton cumulant ratio is measured to be $C_4/C_2=-0.85 \pm 0.09 ~(\rm stat.) \pm 0.82 ~(\rm syst.)$, which is less than unity, the Poisson baseline. The hadronic transport UrQMD model reproduces our $C_4/C_2$ in the measured acceptance. Compared to higher energy results and the transport model calculations, the suppression in $C_4/C_2$ is consistent with fluctuations driven by baryon number conservation and indicates an energy regime dominated by hadronic interactions. These data imply that the QCD critical region, if created in heavy-ion collisions, could only exist at energies higher than 3 GeV.

We report precision measurements of hypernuclei ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ lifetimes obtained from Au+Au collisions at \snn = 3.0 GeV and 7.2 GeV collected by the STAR experiment at RHIC, and the first measurement of ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ mid-rapidity yields in Au+Au collisions at \snn = 3.0 GeV. ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$, being the two simplest bound states composed of hyperons and nucleons, are cornerstones in the field of hypernuclear physics. Their lifetimes are measured to be $221\pm15(\rm stat.)\pm19(\rm syst.)$ ps for ${}^3_\Lambda \rm{H}$ and $218\pm6(\rm stat.)\pm13(\rm syst.)$ ps for ${}^4_\Lambda \rm{H}$. The $p_T$-integrated yields of ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ are presented in different centrality and rapidity intervals. It is observed that the shape of the rapidity distribution of ${}^4_\Lambda \rm{H}$ is different for 0--10% and 10--50% centrality collisions. Thermal model calculations, using the canonical ensemble for strangeness, describes the ${}^3_\Lambda \rm{H}$ yield well, while underestimating the ${}^4_\Lambda \rm{H}$ yield. Transport models, combining baryonic mean-field and coalescence (JAM) or utilizing dynamical cluster formation via baryonic interactions (PHQMD) for light nuclei and hypernuclei production, approximately describe the measured ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ yields. Our measurements provide means to precisely assess our understanding of the fundamental baryonic interactions with strange quarks, which can impact our understanding of more complicated systems involving hyperons, such as the interior of neutron stars or exotic hypernuclei.

We report on the measurements of directed flow $v_1$ and elliptic flow $v_2$ for hadrons ($\pi^{\pm}$, $K^{\pm}$, $K_{S}^0$, $p$, $\phi$, $\Lambda$ and $\Xi^{-}$) from Au+Au collisions at $\sqrt{s_{NN}}$ = 3 GeV and $v_{2}$ for ($\pi^{\pm}$, $K^{\pm}$, $p$ and $\overline{p}$) at 27 and 54.4 GeV with the STAR experiment. While at the two higher energy midcentral collisions the number-of-constituent-quark (NCQ) scaling holds, at 3 GeV the $v_{2}$ at midrapidity is negative for all hadrons and the NCQ scaling is absent. In addition, the $v_1$ slopes at midrapidity for almost all observed hadrons are found to be positive, implying dominant repulsive baryonic interactions. The features of negative $v_2$ and positive $v_1$ slope at 3 GeV can be reproduced with a baryonic mean-field in transport model calculations. These results imply that the medium in such collisions is likely characterized by baryonic interactions.

In high-energy heavy-ion collisions, partonic collectivity is evidenced by the constituent quark number scaling of elliptic flow anisotropy for identified hadrons. A breaking of this scaling and dominance of baryonic interactions is found for identified hadron collective flow measurements in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. In this paper, we report measurements of the first- and second-order azimuthal anisotropic parameters, $v_1$ and $v_2$, of light nuclei ($d$, $t$, $^{3}$He, $^{4}$He) produced in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions at the STAR experiment. An atomic mass number scaling is found in the measured $v_1$ slopes of light nuclei at mid-rapidity. For the measured $v_2$ magnitude, a strong rapidity dependence is observed. Unlike $v_2$ at higher collision energies, the $v_2$ values at mid-rapidity for all light nuclei are negative and no scaling is observed with the atomic mass number. Calculations by the Jet AA Microscopic Transport Model (JAM), with baryonic mean-field plus nucleon coalescence, are in good agreement with our observations, implying baryonic interactions dominate the collective dynamics in 3 GeV Au+Au collisions at RHIC.

Global hyperon polarization, $\overline{P}_\mathrm{H}$, in Au+Au collisions over a large range of collision energy, $\sqrt{s_\mathrm{NN}}$, was recently measured and successfully reproduced by hydrodynamic and transport models with intense fluid vorticity of the quark-gluon plasma. While naïve extrapolation of data trends suggests a large $\overline{P}_\mathrm{H}$ as the collision energy is reduced, the behavior of $\overline{P}_\mathrm{H}$ at small $\sqrt{s_\mathrm{NN}}<7.7$ GeV is unknown. Operating the STAR experiment in fixed-target mode, we measured the polarization of $\Lambda$ hyperons along the direction of global angular momentum in Au+Au collisions at $\sqrt{s_\mathrm{NN}}=3$ GeV. The observation of substantial polarization of $4.91\pm0.81(\rm stat.)\pm0.15(\rm syst.)$% in these collisions may require a reexamination of the viscosity of any fluid created in the collision, of the thermalization timescale of rotational modes, and of hadronic mechanisms to produce global polarization.

We present the centrality dependent measurement of multiplicity and pseudorapidity distributions of charged particles and photons in Au + Au collisions at sqrt{s_NN} = 62.4 GeV. The charged particles and photons are measured in the pseudorapidity region 2.9 < eta < 3.9 and 2.3 < eta < 3.7, respectively. We have studied the scaling of particle production with the number of participating nucleons and the number of binary collisions. The photon and charged particle production in the measured pseudorapidity range has been shown to be consistent with energy independent limiting fragmentation behavior. The photons are observed to follow a centrality independent limiting fragmentation behavior while for the charged particles it is centrality dependent. We have carried out a comparative study of the pseudorapidity distributions of positively charged hadrons, negatively charged hadrons, photons, pions, net protons in nucleus--nucleus collisions and pseudorapidity distributions from p+p collisions. From these comparisons we conclude that baryons in the inclusive charged particle distribution are responsible for the observed centrality dependence of limiting fragmentation. The mesons are found to follow an energy independent behavior of limiting fragmentation while the behavior of baryons seems to be energy dependent.

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in PbPb collisions at nucleon-nucleon center-of-mass energies $\sqrt{s_\mathrm{NN}} =$ 2.76 and 5.02 TeV with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

We present the first measurement of the two-particle transverse momentum differential correlation function, $P_2\equiv\langle \Delta p_{\rm T} \Delta p_{\rm T} \rangle /\langle p_{\rm T} \rangle^2$, in Pb-Pb collisions at $\sqrt{s_{_{\rm NN}}} =$ 2.76 TeV. Results for $P_2$ are reported as a function of relative pseudorapidity ($\Delta \eta$) and azimuthal angle ($\Delta \varphi$) between two particles for different collision centralities. The $\Delta \phi$ dependence is found to be largely independent of $\Delta \eta$ for $|\Delta \eta| \geq$ 0.9. In 5% most central Pb-Pb collisions, the two-particle transverse momentum correlation function exhibits a clear double-hump structure around $\Delta \varphi = \pi$ (i.e., on the away side), which is not observed in number correlations in the same centrality range, and thus provides an indication of the dominance of triangular flow in this collision centrality. Fourier decompositions of $P_2$, studied as a function of collision centrality, show that correlations at $|\Delta \eta| \geq$ 0.9 can be well reproduced by a flow ansatz based on the notion that measured momentum correlations are strictly determined by the collective motion of the system.

We present results on transverse momentum ($p_{\rm T}$) and rapidity ($y$) differential production cross sections, mean transverse momentum and mean transverse momentum square of inclusive J/$\psi$ and $\psi(2S)$ at forward rapidity ($2.5<y<4$) as well as $\psi(2S)$-to-J/$\psi$ cross section ratios. These quantities are measured in pp collisions at center of mass energies $\sqrt{s}=5.02$ and 13 TeV with the ALICE detector. Both charmonium states are reconstructed in the dimuon decay channel, using the muon spectrometer. A comprehensive comparison to inclusive charmonium cross sections measured at $\sqrt{s}=2.76$, 7 and 8 TeV is performed. A comparison to non-relativistic quantum chromodynamics and fixed-order next-to-leading logarithm calculations, which describe prompt and non-prompt charmonium production respectively, is also presented. A good description of the data is obtained over the full $p_{\rm T}$ range, provided that both contributions are summed. In particular, it is found that for $p_{\rm T}>15$ GeV/$c$ the non-prompt contribution reaches up to 50% of the total charmonium yield.

The ALICE Collaboration has measured inclusive J/psi production in pp collisions at a center of mass energy sqrt(s)=2.76 TeV at the LHC. The results presented in this Letter refer to the rapidity ranges |y|<0.9 and 2.5<y<4 and have been obtained by measuring the electron and muon pair decay channels, respectively. The integrated luminosities for the two channels are L^e_int=1.1 nb^-1 and L^mu_int=19.9 nb^-1, and the corresponding signal statistics are N_J/psi^e+e-=59 +/- 14 and N_J/psi^mu+mu-=1364 +/- 53. We present dsigma_J/psi/dy for the two rapidity regions under study and, for the forward-y range, d^2sigma_J/psi/dydp_t in the transverse momentum domain 0<p_t<8 GeV/c. The results are compared with previously published results at sqrt(s)=7 TeV and with theoretical calculations.

Transverse momentum ($p_{\rm{T}}$) spectra of pions, kaons, and protons up to $p_{\rm{T}} = 20$ GeV/$c$ have been measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV using the ALICE detector for six different centrality classes covering 0-80%. The proton-to-pion and the kaon-to-pion ratios both show a distinct peak at $p_{\rm{T}} \approx 3$ GeV/$c$ in central Pb-Pb collisions that decreases towards more peripheral collisions. For $p_{\rm{T}} > 10$ GeV/$c$, the nuclear modification factor is found to be the same for all three particle species in each centrality interval within systematic uncertainties of 10-20%. This suggests there is no direct interplay between the energy loss in the medium and the particle species composition in the hard core of the quenched jet. For $p_{\rm{T}} < 10$ GeV/$c$, the data provide important constraints for models aimed at describing the transition from soft to hard physics.

The second-order azimuthal anisotropy Fourier harmonics, v2, are obtained in pPb and PbPb collisions over a wide pseudorapidity (eta) range based on correlations among six or more charged particles. The pPb data, corresponding to an integrated luminosity of 35 inverse nanobarns, were collected during the 2013 LHC pPb run at a nucleon-nucleon center-of-mass energy of 5.02 TeV by the CMS experiment. A sample of semi-peripheral PbPb collision data at sqrt(s[NN])= 2.76 TeV, corresponding to an integrated luminosity of 2.5 inverse microbarns and covering a similar range of particle multiplicities as the pPb data, is also analyzed for comparison. The six- and eight-particle cumulant and the Lee-Yang zeros methods are used to extract the v2 coefficients, extending previous studies of two- and four-particle correlations. For both the pPb and PbPb systems, the v2 values obtained with correlations among more than four particles are consistent with previously published four-particle results. These data support the interpretation of a collective origin for the previously observed long-range (large Delta[eta]) correlations in both systems. The ratios of v2 values corresponding to correlations including different numbers of particles are compared to theoretical predictions that assume a hydrodynamic behavior of a pPb system dominated by fluctuations in the positions of participant nucleons. These results provide new insights into the multi-particle dynamics of collision systems with a very small overlapping region.

The first measurement of two-pion Bose-Einstein correlations in central Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV at the Large Hadron Collider is presented. We observe a growing trend with energy now not only for the longitudinal and the outward but also for the sideward pion source radius. The pion homogeneity volume and the decoupling time are significantly larger than those measured at RHIC.

The differential cross section for the process $\gamma \gamma \to \eta \pi^0$ has been measured in the kinematic range $0.84 \GeV < W < 4.0 \GeV$, $|\cos \theta^*|<0.8$, where $W$ and $\theta^*$ are the energy and $\pi^0$ (or $\eta$) scattering angle, respectively, in the $\gamma\gamma$ center-of-mass system. The results are based on a 223 fb$^{-1}$ data sample collected with the Belle detector at the KEKB $e^+ e^-$ collider. Clear peaks due to the $a_0(980)$ and $a_2(1320)$ are visible. The differential cross sections are fitted in the energy region $0.9 \GeV < W < 1.46 \GeV$ to obtain the parameters of the $a_0(980)$. Its mass, width and $\Gamma_{\gamma \gamma} \B (\eta \pi^0)$ are measured to be $982.3 ^{+0.6}_{-0.7} ^{+3.1}_{-4.7} \MeV/c^2$, $75.6 \pm 1.6 ^{+17.4}_{-10.0} \MeV$ and $128 ^{+3}_{-2} ^{+502}_{-43} \eV$, respectively. The energy and angular dependences above 3.1 GeV are compared with those measured in the $\pi^0 \pi^0$ channel. The integrated cross section over $|\cos \theta^*|<0.8$ has a $W^{-n}$ dependence with $n = 10.5 \pm 1.2 \pm 0.5$, which is slightly larger than that for $\pi^0 \pi^0$. The differential cross sections show a $\sin^{-4} \theta^*$ dependence similar to $\gamma \gamma \to \pi^0 \pi^0$. The measured cross section ratio, $\sigma(\eta \pi^0)/\sigma(\pi^0 \pi^0) = 0.48 \pm 0.05 \pm 0.04$, is consistent with a QCD-based prediction.

We report a high-statistics measurement of differential cross sections for the process gamma gamma -> pi^0 pi^0 in the kinematic range 0.6 GeV <= W <= 4.0 GeV and |cos theta*| <= 0.8, where W and theta* are the energy and pion scattering angle, respectively, in the gamma gamma center-of-mass system. Differential cross sections are fitted to obtain information on S, D_0, D_2, G_0 and G_2 waves. The G waves are important above W ~= 1.6 GeV. For W <= 1.6 GeV the D_2 wave is dominated by the f_2(1270) resonance while the S wave requires at least one additional resonance besides the f_0(980), which may be the f_0(1370) or f_0(1500). The differential cross sections are fitted with a simple parameterization to determine the parameters (the mass, total width and Gamma_{gamma gamma}B(f_0 -> pi^0 pi^0)) of this scalar meson as well as the f_0(980). The helicity 0 fraction of the f_2(1270) meson, taking into account interference for the first time, is also obtained.

The reactions gamma gamma -> pi^+pi^-pi^+pi^- and gamma gamma -> pi^+pi^0pi^-pi^0 are studied with the L3 detector at LEP in a data sample collected at centre-of-mass energies from 161GeV to 209GeV with a total integrated luminosity of 698/pb. A spin-parity-helicity analysis of the rho^0 rho^0 and rho^+ rho^- systems for two-photon centre-of-mass energies between 1GeV and 3GeV shows the dominance of the spin-parity state 2+ with helicity 2. The contribution of 0+ and 0- spin-parity states is also observed, whereas contributions of 2- states and of a state with spin-parity 2+ and zero helicity are found to be negligible.

Exclusive rho^+ rho^- production in two-photon collisions involving a single highly-virtual photon is studied with data collected at LEP at centre-of-mass energies 89 GeV &lt; \sqrt{s} &lt; 209 GeV with a total integrated luminosity of 854.7 pb^-1. The cross section of the process gamma gamma^* -> rho^+ rho^- is determined as a function of the photon virtuality, Q^2, and the two-photon centre-of-mass energy, W_gg, in the kinematic region: 1.2 GeV^2 &lt; Q^2 &lt; 30 GeV^2 and 1.1 GeV &lt; W_gg &lt; 3 GeV. The \rho^+\rho^- production cross section is found to be of the same magnitude as the cross section of the process gamma gamma^* -> rho^0 rho^0, measured in the same kinematic region by L3, and to have similar W_gg and Q^2 dependences.