We study correlation functions of spatially separated static quark-antiquark pairs in 2+1 flavor QCD in order to investigate onset and nature of color screening at high temperatures. We perform lattice calculations in a wide temperature range, $140\,{\rm MeV}\le T \le 5814\,{\rm MeV}$, using the highly improved staggered quark (HISQ) action and several lattice spacings to control discretization effects. By comparing at high temperatures our lattice results to weak-coupling calculations as well as to the zero temperature result for the energy of a static quark-antiquark pair, we observe that color screening sets in at $rT \approx 0.3$. Furthermore, we also observe that in the range $0.3 \lesssim r T \lesssim 0.6$ weak-coupling calculations in the framework of suitable effective field theories provide an adequate picture of color screening.

We report on the measurement of the all-particle cosmic ray energy spectrum with the High Altitude Water Cherenkov (HAWC) Observatory in the energy range 10 to 500 TeV. HAWC is a ground based air-shower array deployed on the slopes of Volcan Sierra Negra in the state of Puebla, Mexico, and is sensitive to gamma rays and cosmic rays at TeV energies. The data used in this work were taken from 234 days between June 2016 to February 2017. The primary cosmic-ray energy is determined with a maximum likelihood approach using the particle density as a function of distance to the shower core. Introducing quality cuts to isolate events with shower cores landing on the array, the reconstructed energy distribution is unfolded iteratively. The measured all-particle spectrum is consistent with a broken power law with an index of $-2.49\pm0.01$ prior to a break at $(45.7\pm0.1$) TeV, followed by an index of $-2.71\pm0.01$. The spectrum also respresents a single measurement that spans the energy range between direct detection and ground based experiments. As a verification of the detector response, the energy scale and angular resolution are validated by observation of the cosmic ray Moon shadow's dependence on energy.

A simple but novel driver system has been developed to operate the wire gating grid of a Time Projection Chamber (TPC). This system connects the wires of the gating grid to its driver via low impedance transmission lines. When the gating grid is open, all wires have the same voltage allowing drift electrons, produced by the ionization of the detector gas molecules, to pass through to the anode wires. When the grid is closed, the wires have alternating higher and lower voltages causing the drift electrons to terminate at the more positive wires. Rapid opening of the gating grid with low pickup noise is achieved by quickly shorting the positive and negative wires to attain the average bias potential with N-type and P-type MOSFET switches. The circuit analysis and simulation software SPICE shows that the driver restores the gating grid voltage to 90% of the opening voltage in less than 0.20 $\mu$s. When tested in the experimental environment of a time projection chamber larger termination resistors were chosen so that the driver opens the gating grid in 0.35 $\mu$s. In each case, opening time is basically characterized by the RC constant given by the resistance of the switches and terminating resistors and the capacitance of the gating grid and its transmission line. By adding a second pair of N-type and P-type MOSFET switches, the gating grid is closed by restoring 99% of the original charges to the wires within 3 $\mu$s.

A simple but novel driver system has been developed to operate the wire gating grid of a Time Projection Chamber (TPC). This system connects the wires of the gating grid to its driver via low impedance transmission lines. When the gating grid is open, all wires have the same voltage allowing drift electrons, produced by the ionization of the detector gas molecules, to pass through to the anode wires. When the grid is closed, the wires have alternating higher and lower voltages causing the drift electrons to terminate at the more positive wires. Rapid opening of the gating grid with low pickup noise is achieved by quickly shorting the positive and negative wires to attain the average bias potential with N-type and P-type MOSFET switches. The circuit analysis and simulation software SPICE shows that the driver restores the gating grid voltage to 90% of the opening voltage in less than 0.20 $\mu$s. When tested in the experimental environment of a time projection chamber larger termination resistors were chosen so that the driver opens the gating grid in 0.35 $\mu$s. In each case, opening time is basically characterized by the RC constant given by the resistance of the switches and terminating resistors and the capacitance of the gating grid and its transmission line. By adding a second pair of N-type and P-type MOSFET switches, the gating grid is closed by restoring 99% of the original charges to the wires within 3 $\mu$s.

We present the widest-field resolved stellar map to date of the closest ($D\sim3.8$ Mpc) massive elliptical galaxy NGC 5128 (Centaurus A; Cen A), extending out to a projected galactocentric radius of $\sim150$ kpc. The dataset is part of our ongoing Panoramic Imaging Survey of Centaurus and Sculptor (PISCeS) utilizing the Magellan/Megacam imager. We resolve a population of old red giant branch stars down to $\sim1.5$ mag below the tip of the red giant branch, reaching surface brightness limits as low as $\mu_{V,0}\sim32$ mag arcsec$^{-2}$. The resulting spatial stellar density map highlights a plethora of previously unknown streams, shells, and satellites, including the first tidally disrupting dwarf around Cen A (CenA-MM-Dw3), which underline its active accretion history. We report 13 previously unknown dwarf satellite candidates, of which 9 are confirmed to be at the distance of Cen A (the remaining 4 are not resolved into stars), with magnitudes in the range $M_V=-7.2$ to $-13.0$, central surface brightness values of $\mu_{V,0}=25.4-26.9$ mag arcsec$^{-2}$, and half-light radii of $r_h=0.22-2.92$ kpc. These values are in line with Local Group dwarfs but also lie at the faint/diffuse end of their distribution; interestingly, CenA-MM-Dw3 has similar properties to the recently discovered ultra-diffuse galaxies in Virgo and Coma. Most of the new dwarfs are fainter than the previously known Cen A satellites. The newly discovered dwarfs and halo substructures are discussed in light of their stellar populations, and they are compared to those discovered by the PAndAS survey of M31.

In the wide-field Panoramic Imaging Survey of Centaurus and Sculptor (PISCeS), we investigate the resolved stellar halos of two nearby galaxies (the elliptical Centaurus A and the spiral Sculptor, D $\sim3.7$ Mpc) out to a projected galactocentric radius of 150 kpc with Magellan/Megacam. The survey has led to the discovery of $\sim$20 faint satellites to date, plus prominent streams and substructures in two environments that are substantially different from the Local Group, i.e. the Centaurus A group dominated by an elliptical and the loose Sculptor group of galaxies. These discoveries clearly attest to the importance of past and ongoing accretion processes in shaping the halos of these nearby galaxies, and provide the first census of their satellite systems down to an unprecedented $M_V<-8$. The detailed characterization of the stellar content, shape and gradients in the extended halos of Sculptor, Centaurus A, and their dwarf satellites provides key constraints on theoretical models of galaxy formation and evolution.

As part of the Panoramic Imaging Survey of Centaurus and Sculptor (PISCeS) we report the discovery of a pair of faint dwarf galaxies (CenA-MM-Dw1 and CenA-MM-Dw2) at a projected distance of $\sim$90 kpc from the nearby elliptical galaxy NGC5128 (CenA). We measure a tip of the red giant branch distance to each dwarf, finding $D=3.63 \pm 0.41$ Mpc for CenA-MM-Dw1 and $D=3.60 \pm 0.41$ Mpc for CenA-MM-Dw2, both of which are consistent with the distance to NGC5128. A qualitative analysis of the color magnitude diagrams indicates stellar populations consisting of an old, metal-poor red giant branch ($\gtrsim 12$ Gyr, [Fe/H]$\sim-1.7$ to -1.9). In addition, CenA-MM-Dw1 seems to host an intermediate-age population as indicated by its candidate asymptotic giant branch stars. The derived luminosities ($M_V=-10.9\pm0.3$ for CenA-MM-Dw1 and $-8.4\pm0.6$ for CenA-MM-Dw2) and half-light radii ($r_{h}=1.4\pm0.04$ kpc for CenA-MM-Dw1 and $0.36\pm0.08$ kpc for CenA-MM-Dw2) are consistent with those of Local Group dwarfs. CenA-MM-Dw1's low central surface brightness ($\mu_{V,0}=27.3\pm0.1$ mag/arcsec$^2$) places it among the faintest and most extended M31 satellites. Most intriguingly, CenA-MM-Dw1 and CenA-MM-Dw2 have a projected separation of only 3 arcmin ($\sim3$ kpc): we are possibly observing the first, faint satellite of a satellite in an external group of galaxies.

We introduce the notion of a Mahonian pair. Consider the set, P^*, of all words having the positive integers as alphabet. Given finite subsets S,T of P^*, we say that (S,T) is a Mahonian pair if the distribution of the major index, maj, over S is the same as the distribution of the inversion number, inv, over T. So the well-known fact that maj and inv are equidistributed over the symmetric group, S_n, can be expressed by saying that (S_n,S_n) is a Mahonian pair. We investigate various Mahonian pairs (S,T) with S different from T. Our principal tool is Foata's fundamental bijection f: P^* -> P^* since it has the property that maj w = inv f(w) for any word w. We consider various families of words associated with Catalan and Fibonacci numbers. We show that, when restricted to words in {1,2}^*, f transforms familiar statistics on words into natural statistics on integer partitions such as the size of the Durfee square. The Rogers-Ramanujan identities, the Catalan triangle, and various q-analogues also make an appearance. We generalize the definition of Mahonian pairs to infinite sets and use this as a tool to connect a partition bijection of Corteel-Savage-Venkatraman with the Greene-Kleitman decomposition of a Boolean algebra into symmetric chains. We close with comments about future work and open problems.

A rich harvest of RR Lyrae stars has been identified for the first time in B514, a metal-poor ([Fe/H] = 1.95 +/- 0.10 dex) globular cluster of the Andromeda galaxy (M31), based on Hubble Space Telescope Wide Field Planetary Camera 2 and Advanced Camera for Surveys time-series observations. We have detected and derived periods for 89 RR Lyrae stars (82 fundamental-mode -RRab- and 7 first-overtone -RRc- pulsators, respectively) among 161 candidate variables identified in the cluster. The average period of the RR Lyrae variables (<Pab> = 0.58 days and <Pc> = 0.35 days, for RRab and RRc pulsators, respectively) and the position in the period-amplitude diagram both suggest that B514 is likely an Oosterhoff type I cluster. This appears to be in disagreement with the general behaviour of the metal-poor globular clusters in the Milky Way, which show instead Oosterhoff type II pulsation properties. The average apparent magnitude of the RR Lyrae stars sets the mean level of the cluster horizontal branch at <V(RR)> = 25.18 +/- 0.02 (sigma=0.16 mag, on 81 stars). By adopting a reddening E(B-V) = 0.07 +/- 0.02 mag, the above metallicity and M_V=0.44 +/- 0.05 mag for the RR Lyrae variables of this metallicity, we derive a distance modulus of mu_0=24.52 +/- 0.08 mag, corresponding to a distance of about 800 +/- 30 kpc, based on a value of M_V that sets mu_0(LMC)=18.52.

Features of anisotropic collective flow and spectral temperatures have been determined for lambda hyperons emitted from C + C collisions, at incident momentum of 4.2 AGeV/c, measured using the Propane Bubble Chamber of JINR at Dubna. Moreover, characteristics of protons and of negative pions, emitted from those collisions, have been determined and provided for comparison. The directed and elliptic flows of lambdas both agree in sign with the corresponding flows of protons. Parameters of the directed and elliptic flows for lambdas agree further, within errors, with the corresponding parameters for the co-produced protons. This contrasts an earlier finding by the E895 Collaboration of the directed flow being significantly weaker for lambdas than protons, in the much heavier Au + Au system, at comparable incident momentum. Particle spectral temperatures in the C + C collisions have been determined focusing independently on either center-of-mass energy, transverse energy or transverse momentum distributions. For either protons or negative pions, the temperatures were found to be approximately the same, no matter whether the emission of those particles was associated with lambda production or not. Results of the measurements have been compared to the results of simulations within the Quark-Gluon String Model.

The minimal cooling paradigm for neutron star cooling assumes that enhanced cooling due to neutrino emission from any direct Urca process, due either to nucleons or to exotica such as hyperons, Bose condensates, or deconfined quarks, does not occur. This scenario was developed to replace and extend the so-called standard cooling scenario to include neutrino emission from the Cooper pair breaking and formation processes that occur near the critical temperature for superfluid/superconductor pairing. Recently, it has been found that Cooper-pair neutrino emission from the vector channel is suppressed by a large factor compared to the original estimates that violated vector current conservation. We show that Cooper-pair neutrino emission remains, nevertheless, an efficient cooling mechanism through the axial channel. As a result, the elimination of neutrino emission from Cooper-paired nucleons through the vector channel has only minor effects on the long-term cooling of neutron stars within the minimal cooling paradigm. We further quantify precisely the effect of the size of the neutron 3P2 gap and demonstrate that consistency between observations and the minimal cooling paradigm requires that the critical temperature T_c for this gap covers a range of values between T_c^min < 0.2 x 10^9 K up to T_c^max > 0.5 \times 10^9 K in the core of the star. In addition, it is required that young neutron stars have heterogenous envelope compositions: some must have light-element compositions and others must have heavy-element compositions. Unless these two conditions are fulfilled, about half of the observed young cooling neutron stars are inconsistent with the minimal cooling paradigm and provide evidence for the existence of enhanced cooling.

It is proposed that a study of the ratio of shear viscosity to entropy density $\frac{\eta}{s}$ as a function of the baryon chemical potential $\mu_B$, and temperature T, provides a dynamic probe for the critical end point (CEP) in hot and dense QCD matter. An initial estimate from an elliptic flow excitation function gives $\mu^{\text{cep}}_B \sim 150-180$ MeV and $T_{\text{cep}} \sim 165 - 170$ MeV for the location of the the CEP. These values place the CEP in the range for "immediate" validation at RHIC.

We report on the detection of variable stars in the Canes Venatici II (CVn II) dwarf spheroidal galaxy, a new satellite of the Milky Way recently discovered by the Sloan Digital Sky Survey. We also present a V, B-V color-magnitude diagram that reaches V = 25.5 mag, showing the galaxy's main sequence turn off at V = 24.5 mag and revealing several candidate blue straggler stars. Two RR Lyrae stars have been identified within the half-light radius of CVn II,a fundamental-mode variable (RRab) with period P_ab = 0.743 days, and a first-overtone (RRc) RR Lyrae star with P_c = 0.358 days. The rather long periods of these variables along with their position on the period-amplitude diagram support an Oosterhoff type II classification for CVn II. The average apparent magnitude of the RR Lyrae stars, <V> = 21.48 +/- 0.02 mag, is used to obtain a precision distance modulus of mu_0 = 21.02 +/- 0.06 mag and a corresponding distance of 160(+4,-5} kpc, for an adopted reddening E(B-V) = 0.015 mag.

We report on the detection of variable stars in the Canes Venatici II (CVn II) dwarf spheroidal galaxy, a new satellite of the Milky Way recently discovered by the Sloan Digital Sky Survey. We also present a V, B-V color-magnitude diagram that reaches V = 25.5 mag, showing the galaxy's main sequence turn off at V = 24.5 mag and revealing several candidate blue straggler stars. Two RR Lyrae stars have been identified within the half-light radius of CVn II,a fundamental-mode variable (RRab) with period P_ab = 0.743 days, and a first-overtone (RRc) RR Lyrae star with P_c = 0.358 days. The rather long periods of these variables along with their position on the period-amplitude diagram support an Oosterhoff type II classification for CVn II. The average apparent magnitude of the RR Lyrae stars, <V> = 21.48 +/- 0.02 mag, is used to obtain a precision distance modulus of mu_0 = 21.02 +/- 0.06 mag and a corresponding distance of 160(+4,-5} kpc, for an adopted reddening E(B-V) = 0.015 mag.

Collective flow of protons and negative pions has been studied within the momentum region of $4.2 \div 4.5$ AGeV/c ($E =3.4 \div 3.7$ AGeV) for different projectile-target combinations involving carbon and, specifically, He-C, C-C, C-Ne, C-Cu and C-Ta. The data stem from the SKM-200-GIBS streamer chamber and from Propane Bubble Chamber systems utilized at JINR. The directed flow of protons grows dramatically in the carbon region when the counterpart nucleus grows in mass between He and Ta. The elliptic proton flow points out of the reaction plane and also strengthens as system mass increases. Within the reaction plane, the negative pions flow in the same direction as protons for the lighter of the investigated systems, He-C, C-C and C-Ne, and in the opposite direction for the heavier, C-Cu and C-Ta. The Quark-Gluon String Model reproduces observed changes in the flow with system mass.

The NA49 experiment at CERN SPS has acquired a huge data set of Pb+Pb events over a broad range of energy and centrality during the last several years. This high statistics data set, coupled with a state-of-the-art analysis technique, allows for the first model-independent extraction and energy scan of 3D emission sources for pion pairs at SPS energies. These 3D pion emission sources provide new insights into the nature of a long-range source previously reported by PHENIX at RHIC. The new results indicate that the pion source displays significant non-Gaussian tails in the longitudinal direction at 40 and 158 AGeV and in the outward direction at 158 AGeV.

We apply time-distance helioseismology, local correlation tracking and Fourier spatial-temporal filtering methods to realistic supergranule scale simulations of solar convection and compare the results with high-resolution observations from the SOHO Michelson Doppler Imager (MDI). Our objective is to investigate the surface and sub-surface convective structures and test helioseismic measurements. The size and grid of the computational domain are sufficient to resolve various convective scales from granulation to supergranulation. The spatial velocity spectrum is approximately a power law for scales larger than granules, with a continuous decrease in velocity amplitude with increasing size. Aside from granulation no special scales exist, although a small enhancement in power at supergranulation scales can be seen. We calculate the time-distance diagram for f- and p-modes and show that it is consistent with the SOHO/MDI observations. From the simulation data we calculate travel time maps for surface gravity waves (f-mode). We also apply correlation tracking to the simulated vertical velocity in the photosphere to calculate the corresponding horizontal flows. We compare both of these to the actual large-scale (filtered) simulation velocities. All three methods reveal similar large scale convective patterns and provide an initial test of time-distance methods.

We report on a programme to monitor the radial velocities of a sample of candidate and confirmed carbon-enhanced, metal-poor (CEMP) stars. We observed 45 targets using the Echelle Spectrographs of three 4m class telescopes. Radial velocities for these objects were calculated by cross-correlation of their spectra with the spectrum of HD 140283, and have errors < 1 km/s. Sixteen of our programme's targets have reported carbon excess, and nine of these objects also exhibit s-process enhancements (CEMP-s). We combine these stars' radial velocities with other literature studies in search of binarity. The search reveals that four of our CEMP-s stars (44%) are in binary systems. Using the analysis of Lucatello et al. (2004), we find that all the CEMP-s stars in our sample are binaries. This conclusion implies that CEMP-s stars may be the very metal-poor relatives of CH and Ba II stars, which are believed to have acquired their peculiar abundance patterns by mass transfer from a thermally-pulsing AGB companion.

In this paper we show how effective parameters such as effective binding energies can be defined for a proton in the combined nuclear-Coulomb potential, including also the target potential, in the case in which the proton is bound in a nucleus which is partner of a nuclear reaction. Using such effective parameters the proton behaves similarly to a neutron. In this way some unexpected results obtained from dynamical calculations for reactions initiated by very weakly bound proton halo nuclei can be interpreted. Namely the fact that stripping dominates the nuclear breakup cross section which in turn dominates over the Coulomb breakup even when the target is heavy at medium to high incident energies. Our interpretation helps also clarifying why the existence and characteristics of a proton halo extracted from different types of data have sometimes appeared contradictory.

A 3D simulation of the upper part of the solar convective zone is used to obtain information on the frequency component, chi_k, of the correlation product of the turbulent velocity field. This component plays an important role in the stochastic excitation of acoustic oscillations. A time analysis of the solar simulation shows that a gaussian function does not correctly reproduce the nu-dependency of chi_k inferred from the 3D simuation in the frequency range where the acoustic energy injected into the solar p modes is important (nu ~ 2 - 4 mHz). The nu-dependency of chi_k is fitted with different analytical functions which can then conveniently be used to compute the acoustic energy supply rate P injected into the solar radial oscillations. With constraints from a 3D simulation, adjustment of free parameters to solar data is no longer necessary and is not performed here. The result is compared with solar seismic data. Computed values of P obtained with the analytical function which fits best chi_k are found ~ 2.7 times larger than those obtained with the gaussian model and reproduce better the solar seismic observations. This non-gaussian description also leads to a Reynolds stress contribution of the same order as the one arising from the advection of the turbulent fluctuations of entropy by the turbulent motions. Some discrepancy between observed and computed P values still exist at high frequency and possible causes for this discrepancy are discussed.

Analyses of a 3D simulation of the upper layers of a solar convective envelope provide constraints on the physical quantities which enter the theoretical formulation of a stochastic excitation model of solar p modes, for instance the convective velocities and the turbulent kinetic energy spectrum. These constraints are then used to compute the acoustic excitation rate for solar p modes, P. The resulting values are found ~5 times larger than the values resulting from a computation in which convective velocities and entropy fluctuations are obtained with a 1D solar envelope model built with the time-dependent, nonlocal Gough (1977) extension of the mixing length formulation for convection (GMLT). This difference is mainly due to the assumed mean anisotropy properties of the velocity field in the excitation region. The 3D simulation suggests much larger horizontal velocities compared to vertical ones than in the 1D GMLT solar model. The values of P obtained with the 3D simulation constraints however are still too small compared with the values inferred from solar observations. Improvements in the description of the turbulent kinetic energy spectrum and its depth dependence yield further increased theoretical values of P which bring them closer to the observations. It is also found that the source of excitation arising from the advection of the turbulent fluctuations of entropy by the turbulent movements contributes ~ 65-75 % to the excitation and therefore remains dominant over the Reynolds stress contribution. The derived theoretical values of P obtained with the 3D simulation constraints remain smaller by a factor ~3 compared with the solar observations. This shows that the stochastic excitation model still needs to be improved.

The defining characteristics of fragment emission resulting from the non-central collision of 114Cd ions with 92Mo target nuclei at E/A = 50 MeV are presented. Charge correlations and average relative velocities for mid-velocity fragment emission exhibit significant differences when compared to standard statistical decay. These differences associated with similar velocity dissipation are indicative of the influence of the entrance channel dynamics on the fragment production process.

The structure of LiMoS2 has been experimentally determined for the first time. The approach of atomic pair distribution function (PDF) analysis was used because of the lack of well defined Bragg peaks due to the short structural coherence (~ 50 A) in this intercalation compound. This represents one of the first instances where the structure of such poorly diffracting material was determined from the total scattering including diffuse as well as Bragg scattering. The reduction of Mo by Li results in Mo-Mo bonding with the formation of chains of distorted Mo-S6 octahedra. Using the refined structural parameters the electronic band-structure for this material has been calculated and is in good agreement with observed material properties.

The structure of LiMoS2 has been experimentally determined for the first time. The approach of atomic pair distribution function (PDF) analysis was used because of the lack of well defined Bragg peaks due to the short structural coherence (~ 50 A) in this intercalation compound. This represents one of the first instances where the structure of such poorly diffracting material was determined from the total scattering including diffuse as well as Bragg scattering. The reduction of Mo by Li results in Mo-Mo bonding with the formation of chains of distorted Mo-S6 octahedra. Using the refined structural parameters the electronic band-structure for this material has been calculated and is in good agreement with observed material properties.

Nearest and higher neighbor distances as well as bond length distributions (static and thermal) of the In_xGa_(1-x)As (0<x<1) semiconductor alloys have been obtained from high real-space resolution atomic pair distribution functions (PDFs). Using this structural information, we modeled the local atomic displacements in In_xGa_(1-x)As alloys. From a supercell model based on the Kirkwood potential, we obtained 3-D As and (In,Ga) ensemble averaged probability distributions. This clearly shows that As atom displacements are highly directional and can be represented as a combination of <100> and <111> displacements. Examination of the Kirkwood model indicates that the standard deviation (sigma) of the static disorder on the (In,Ga) sublattice is around 60% of the value on the As sublattice and the (In,Ga) atomic displacements are much more isotropic than those on the As sublattice. The single crystal diffuse scattering calculated from the Kirkwood model shows that atomic displacements are most strongly correlated along <110> directions.