Source author record

A. Li

A. Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

76works
13topics
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

76 published item(s)

preprint2026arXiv

Evidence for the rare decay $B^{+} \to \bar Λp μ^{+} μ^{-}$

A search for the rare decay $B^{+} \to \bar Λp μ^{+} μ^{-}$ is performed using proton-proton collision data recorded by the LHCb experiment at a center-of-mass energy of $\sqrt{s}= 13$ TeV, corresponding to an integrated luminosity of 5.4 fb$^-1$. An excess of events is found with respect to the background-only expectation, with a signal significance of 3.5 standard deviations, in the low invariant-mass region of $m(\bar Λp)<2.8$ GeV/$c^2$. The branching fraction is measured to be ${\cal B}_{\rm low}({B^{+}}\to {\bar Λp μ^{+} μ^{-}})=\left(1.70 ^{+0.65}_{-0.56}(\rm stat) \pm 0.17(\rm syst) \pm 0.14(\rm ext)\right) \times 10^{-8}$, where the last uncertainty is due to external inputs on ${\cal B}(B^+\to{J/ψ\bar Λp})\times {\cal B}({J/ψ}\to{μ^+μ^-})$. With no significant signal observed in the high $m(\bar Λp)$ region above 2.8 GeV/$c^2$, an upper limit is set to be ${\cal B}_{\rm high}(B^{+}\to {\bar Λp μ^{+} μ^{-}})<2.8\,(3.7) \times 10^{-9}$ at the $90\%$ ($95\%$) confidence level.

preprint2026arXiv

Light Dark Matter Search with 7.8 Tonne-Year of Ionization-Only Data in XENONnT

We report on a blinded search for dark matter (DM) using ionization-only (S2-only) signals in XENONnT with a total exposure of $7.83\mathrm{tonne}\times\mathrm{year}$ over 579 days in three science runs. Dedicated background suppression techniques and the first complete S2-only background model in XENONnT provide sensitivity to nuclear recoils of [0.5, 5.0] $\mathrm{keV_\mathrm{nr}}$ and electronic recoils of [0.04, 0.7] $\mathrm{keV_\mathrm{ee}}$. No significant excess over the expected background is observed, and we set 90\% confidence level upper limits on spin-independent DM--nucleon and spin-dependent DM--neutron scattering for DM masses between 3 and 8 $\mathrm{GeV}/c^2$, as well as on DM--electron scattering, axion-like particles, and dark photons, improving on previous constraints. For spin-independent DM--nucleon scattering, we exclude cross sections above $6.0\times10^{-45} $cm$^2$ at a DM mass of 5 $\mathrm{GeV}/c^2$, pushing the XENONnT sensitivity closer to the region where coherent elastic neutrino-nucleus scattering ($\text{CE}ν\text{NS}$) becomes an irreducible background.

preprint2026arXiv

Measurement of transverse $Λ$ and $\barΛ$ hyperon polarization in $p$Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV

The transverse polarization of $Λ$ and $\barΛ$ hyperons is measured in $p$Pb collisions collected by the LHCb experiment at a nucleon-nucleon center-of-mass energy of $5.02$ TeV. The polarization is averaged over hyperon transverse momentum in the ranges $0.15 < p_{T} < 6.00 $ GeV/$c$, and Feynman-$x$ in the ranges $0.005 < x_{F} < 0.040$ (forward region) and $-0.10 < x_{F} < -0.01$ (backward region) defined relative to the proton beam direction. The transverse polarization is found to be compatible with zero for both $Λ$ and $\barΛ$ hyperons. The results are also measured as a function of $p_{T}$ and $x_{F}$ with no significant dependence on these variables observed. The results are compared with previous experimental measurements at different center-of-mass energies and collision environments.

preprint2026arXiv

Search for lepton-number-violating $B^-\to D^{(*)+}μ^-μ^-$ decays

A search is performed for lepton-number-violating $B^-\to D^{(*)+}μ^-μ^-$ decays, using data collected by the LHCb experiment in proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 5.4 fb$^{-1}$. No significant signal is observed, and upper limits are set on the branching fractions, ${\cal B}(B^-\to D^{+}μ^-μ^-) < 4.6 \times 10^{-8}$ and ${\cal B}(B^-\to D^{*+}μ^-μ^-) < 5.9 \times 10^{-8}$, at the 95% confidence level.

preprint2026arXiv

Study of charm mixing and CP violation with $D^0\to K^\pmπ^\mpπ^\pmπ^\mp$ decays

A study of charm mixing and CP violation in $D^0\to K^\pmπ^\mpπ^\pmπ^\mp$ decays is performed using data collected by the LHCb experiment in proton-proton collisions from 2015 to 2018, corresponding to an integrated luminosity of 6$\text{fb}^{-1}$. The ratio of promptly produced $D^0\to K^+π^- π^+π^-$ to $D^0\to K^-π^+ π^-π^+$ decay rates is measured as a function of $D^0$ decay time, both inclusive over phase space and in bins of phase space. Taking external inputs for the $D^0 -\overline{D}^0$ mixing parameters $x$ and $y$ allows constraints to be obtained on the hadronic parameters of the charm decay. When combined with previous measurements from charm-threshold experiments and at LHCb, improved knowledge is obtained for these parameters, which is valuable for studies of the angle $γ$ of the Unitarity Triangle. An alternative analysis is also performed, in which external inputs are taken for the hadronic parameters, and the mixing parameters are determined, including $Δx$ and $Δy$, which are nonzero in the presence of CP violation. It is found that $x=\left(0.85^{+0.15}_{-0.24}\right)\%$, $y=\left( 0.21^{+0.29}{-0.27} \right)\%$, $Δx=\left( -0.02\pm {0.04} \right)\% $ and $Δy=\left( 0.02^{+0.04}_{-0.03} \right)\%$. These results are consistent with previous measurements and the hypothesis of CP conservation.

preprint2025arXiv

Coherent photoproduction of $ρ^0, ω$ and excited vector mesons in ultraperipheral PbPb collisions

The invariant-mass distribution for the coherent photoproduction of dipions in ultraperipheral PbPb collisions is measured using data, corresponding to an integrated luminosity of $224.6 \pm 9.6 μ$b$^{-1}$, collected by the LHCb experiment in 2018 at a nucleon-nucleon centre-of-mass energy $\sqrt{s_{\rm NN}}=5.02$ TeV. In the mass range from 400 to 1200 MeV, the results are consistent with previous experiments, with the spectrum dominated by the $ρ^0$ meson, which interferes with a nonresonant component, together with a smaller $ω$ meson contribution. In an extended mass range up to 2300 MeV, models previously used do not fit the data and a consistent description requires the introduction of two resonances at masses of $1350\pm20$ MeV and $1790\pm20$ MeV with widths of about 300 MeV. The cross-section for each meson is measured differentially in twelve bins of rapidity from 2.05 to 4.90. The $ρ^0$ cross-section increases with rapidity from about 400 to 600 mb and is measured with a typical precision of 8\%, while the cross-section times branching fraction for the $ω,ρ^\prime$ and $ρ^{\prime\prime}$, with the statistical precision of the data, do not have a pronounced rapidity dependence and are between 0.5 and 1.5 mb, with uncertainties up to 30\%. A large nuclear suppression is observed for the $ρ^0$ meson compared to expectations based on photoproduction on the proton that use the impulse approximation. Significant suppression is also observed compared to that predicted by elastic scattering described in the Glauber approach, or with the addition of inelastic scattering in a Gribov--Glauber model.

preprint2023arXiv

Modeling Backgrounds for the MAJORANA DEMONSTRATOR

The MAJORANA DEMONSTRATOR is a neutrinoless double-beta decay ($0νββ$) experiment containing $\sim$30 kg of p-type point contact germanium detectors enriched to 88% in 76Ge and $\sim$14 kg of natural germanium detectors. The detectors are housed in two electroformed copper cryostats and surrounded by a graded passive shield with active muon veto. An extensive radioassay campaign was performed prior to installation to insure the use of ultra-clean materials. The DEMONSTRATOR achieved one of the lowest background rates in the region of the $0νββ$ Q-value, 15.7 $\pm$ 1.4 cts/(FWHM t y) from the low-background configuration spanning most of the 64.5 kg-yr active exposure. Nevertheless this background rate is a factor of five higher than the projected background rate. This discrepancy arises from an excess of events from the 232Th decay chain. Background model fits aim to understand this deviation from assay-based projections, potentially determine the source(s) of observed backgrounds, and allow a precision measurement of the two-neutrino double-beta decay half-life. The fits agree with earlier simulation studies, which indicate the origin of the 232Th excess is not from a near-detector component and have informed design decisions for the next-generation LEGEND experiment. Recent findings have narrowed the suspected locations for the excess activity, motivating a final simulation and assay campaign to complete the background model.

preprint2022arXiv

A search for correlated low-energy electron antineutrinos in KamLAND with gamma-ray bursts

We present the results of a time-coincident event search for low-energy electron antineutrinos in the KamLAND detector with gamma-ray bursts from the Gamma-ray Coordinates Network and Fermi Gamma-ray Burst Monitor. Using a variable coincidence time window of $\pm$500s plus the duration of each gamma-ray burst, no statistically significant excess above background is observed. We place the world's most stringent 90% confidence level upper limit on the electron antineutrino fluence below 17.5 MeV. Assuming a Fermi-Dirac neutrino energy spectrum from the gamma-ray burst source, we use the available redshift data to constrain the electron antineutrino luminosity and effective temperature.

preprint2022arXiv

Abundances of uranium and thorium elements in Earth estimated by geoneutrino spectroscopy

The decay of the primordial isotopes $^{238}\mathrm{U}$, $^{235}\mathrm{U}$, $^{232}\mathrm{Th}$, and $^{40}\mathrm{K}$ have contributed to the terrestrial heat budget throughout the Earth's history. Hence the individual abundance of those isotopes are key parameters in reconstructing contemporary Earth model. The geoneutrinos produced by the radioactive decays of uranium and thorium have been observed with the Kamioka Liquid-Scintillator Antineutrino Detector (KamLAND). Those measurements have been improved with more than 18-year observation time, and improvements in detector background levels mainly by an 8-year almost rector-free period now permit spectroscopy with geoneutrinos. Our results yield the first constraint on both uranium and thorium heat contributions. Herein the KamLAND result is consistent with geochemical estimations based on elemental abundances of chondritic meteorites and mantle peridotites. The High-Q model is disfavored at 99.76% C.L. and a fully radiogenic model is excluded at 5.2$σ$ assuming a homogeneous heat producing element distribution in the mantle.

preprint2022arXiv

Constraints on Lightly Ionizing Particles from CDMSlite

The Cryogenic Dark Matter Search low ionization threshold experiment (CDMSlite) achieved efficient detection of very small recoil energies in its germanium target, resulting in sensitivity to Lightly Ionizing Particles (LIPs) in a previously unexplored region of charge, mass, and velocity parameter space. We report first direct-detection limits calculated using the optimum interval method on the vertical intensity of cosmogenically-produced LIPs with an electric charge smaller than $e/(3\times10^5$), as well as the strongest limits for charge $\leq e/160$, with a minimum vertical intensity of $1.36\times10^{-7}$\,cm$^{-2}$s$^{-1}$sr$^{-1}$ at charge $e/160$. These results apply over a wide range of LIP masses (5\,MeV/$c^2$ to 100\,TeV/$c^2$) and cover a wide range of $βγ$ values (0.1 -- $10^6$), thus excluding non-relativistic LIPs with $βγ$ as small as 0.1 for the first time.

preprint2022arXiv

Effective Field Theory Analysis of CDMSlite Run 2 Data

CDMSlite Run 2 was a search for weakly interacting massive particles (WIMPs) with a cryogenic 600 g Ge detector operated in a high-voltage mode to optimize sensitivity to WIMPs of relatively low mass from 2 - 20 GeV/$c^2$. In this article, we present an effective field theory (EFT) analysis of the CDMSlite Run 2 data using an extended energy range and a comprehensive treatment of the expected background. A binned likelihood Bayesian analysis was performed on the recoil energy data, taking into account the parameters of the EFT interactions and optimizing the data selection with respect to the dominant background components. Energy regions within 5$σ$ of known activation peaks were removed from the analysis. The Bayesian evidences resulting from the different operator hypotheses show that the CDMSlite Run 2 data are consistent with the background-only models and do not allow for a signal interpretation assuming any additional EFT interaction. Consequently, upper limits on the WIMP mass and coupling-coefficient amplitudes and phases are presented for each EFT operator. These limits improve previous CDMSlite Run 2 bounds for WIMP masses above 5 GeV/$c^2$.

preprint2022arXiv

Experimental study of 13C(α,n)16O reactions in the Majorana Demonstrator calibration data

Neutron captures and delayed decays of reaction products are common sources of backgrounds in ultra-rare event searches. In this work, we studied $^{13}$C($α,n)^{16}$O reactions induced by $α$-particles emitted within the calibration sources of the \textsc{Majorana Demonstrator}. These sources are thorium-based calibration standards enclosed in carbon-rich materials. The reaction rate was estimated by using the 6129-keV $γ$-rays emitted from the excited $^{16}$O states that are populated when the incoming $α$-particles exceed the reaction Q-value. Thanks to the excellent energy performance of the \textsc{Demonstrator}'s germanium detectors, these characteristic photons can be clearly observed in the calibration data. Facilitated by \textsc{Geant4} simulations, a comparison between the observed 6129-keV photon rates and predictions by a TALYS-based software was performed. The measurements and predictions were found to be consistent, albeit with large statistical uncertainties. This agreement provides support for background projections from ($α,n$)-reactions in future double-beta decay search efforts.

preprint2022arXiv

KamLAND's search for correlated low-energy electron antineutrinos with astrophysical neutrinos from IceCube

We report the results of a search for MeV-scale astrophysical neutrinos in KamLAND presented as an excess in the number of coincident neutrino interactions associated with the publicly available high-energy neutrino datasets from the IceCube Neutrino Observatory. We find no statistically significant excess in the number of observed low-energy electron antineutrinos in KamLAND, given a coincidence time window of $\pm$500s, $\pm$1,000s, $\pm$3,600s, and $\pm$10,000s around each of the IceCube neutrinos. We use this observation to present limits from 1.8 MeV to 100 MeV on the electron antineutrino fluence, assuming a mono-energetic flux. We then compare the results to several astrophysical measurements performed by IceCube and place a limit at the 90% confidence level on the electron antineutrino isotropic thermal luminosity from the TXS 0506+056 blazar.

preprint2022arXiv

Response of a CMS HGCAL silicon-pad electromagnetic calorimeter prototype to 20-300 GeV positrons

The Compact Muon Solenoid Collaboration is designing a new high-granularity endcap calorimeter, HGCAL, to be installed later this decade. As part of this development work, a prototype system was built, with an electromagnetic section consisting of 14 double-sided structures, providing 28 sampling layers. Each sampling layer has an hexagonal module, where a multipad large-area silicon sensor is glued between an electronics circuit board and a metal baseplate. The sensor pads of approximately 1 cm$^2$ are wire-bonded to the circuit board and are readout by custom integrated circuits. The prototype was extensively tested with beams at CERN's Super Proton Synchrotron in 2018. Based on the data collected with beams of positrons, with energies ranging from 20 to 300 GeV, measurements of the energy resolution and linearity, the position and angular resolutions, and the shower shapes are presented and compared to a detailed Geant4 simulation.

preprint2022arXiv

Search for solar axions via axion-photon coupling with the Majorana Demonstrator

Axions were originally proposed to explain the strong-CP problem in QCD. Through the axion-photon coupling, the Sun could be a major source of axions, which could be measured in solid state detection experiments with enhancements due to coherent Primakoff-Bragg scattering. The Majorana Demonstrator experiment has searched for solar axions with a set of $^{76}$Ge-enriched high purity germanium detectors using a 33 kg-yr exposure collected between Jan. 2017 and Nov. 2019. A temporal-energy analysis gives a new limit on the axion-photon coupling as $g_{aγ}<1.45\times 10^{-9}$ GeV$^{-1}$ (95% C.I.) for axions with mass up to 100 eV/$c^2$. This improves laboratory-based limits between about 1 eV/$c^2$ and 100 eV/$c^2$.

preprint2022arXiv

Search for supernova neutrinos and constraint on the galactic star formation rate with the KamLAND data

We present the results of a search for core-collapse supernova neutrinos, using long-term KamLAND data from 2002 March 9 to 2020 April 25. We focus on the electron antineutrinos emitted from supernovae in the energy range of 1.8--111 MeV. Supernovae will make a neutrino event cluster with the duration of $\sim$10 s in the KamLAND data. We find no neutrino clusters and give the upper limit on the supernova rate as to be 0.15 yr$^{-1}$ with a 90% confidence level. The detectable range, which corresponds to a >95% detection probability, is 40--59 kpc and 65--81 kpc for core-collapse supernovae and failed core-collapse supernovae, respectively. This paper proposes to convert the supernova rate obtained by the neutrino observation to the Galactic star formation rate. Assuming a modified Salpeter-type initial mass function, the upper limit on the Galactic star formation rate is <(17.5--22.7) $M_{\odot} \mathrm{yr}^{-1}$ with a 90% confidence level.

preprint2022arXiv

The MAJORANA DEMONSTRATOR Readout Electronics System

The MAJORANA DEMONSTRATOR comprises two arrays of high-purity germanium detectors constructed to search for neutrinoless double-beta decay in 76-Ge and other physics beyond the Standard Model. Its readout electronics were designed to have low electronic noise, and radioactive backgrounds were minimized by using low-mass components and low-radioactivity materials near the detectors. This paper provides a description of all components of the MAJORANA DEMONSTRATOR readout electronics, spanning the front-end electronics and internal cabling, back-end electronics, digitizer, and power supplies, along with the grounding scheme. The spectroscopic performance achieved with these readout electronics is also demonstrated.

preprint2021arXiv

A Search for Charged Excitation of Dark Matter with the KamLAND-Zen Detector

There are many theories where a dark matter particle is part of a multiplet with an electrically charged state. If WIMP dark matter ($χ^{0}$) is accompanied by a charged excited state ($χ^{-}$) separated by a small mass difference, it can form a stable bound state with a nucleus. In supersymmetric models, the $χ^{0}$ and the $χ^{-}$ could be the neutralino and a charged slepton, such as the neutralino-stau degenerate model. The formation binding process is expected to result in an energy deposition of {\it O}(1--10 MeV), making it suitable for detection in large liquid scintillator detectors. We describe new constraints on the bound state formation with a xenon nucleus using the KamLAND-Zen 400 Phase-II dataset. In order to enlarge the searchable parameter space, all xenon isotopes in the detector were used. For a benchmark parameter set of $m_{χ^{0}} = 100$ GeV and $Δm = 10$ MeV, this study sets the most stringent upper limits on the recombination cross section $\langleσv\rangle$ and the decay-width of $χ^{-}$ of $2.0 \times 10^{-31}$ ${\rm cm^3/s}$ and $1.1 \times 10^{-18}$ GeV, respectively (90\% confidence level).

preprint2021arXiv

Constraints on dark photons and axion-like particles from SuperCDMS Soudan

We present an analysis of electron recoils in cryogenic germanium detectors operated during the SuperCDMS Soudan experiment. The data are used to set new constraints on the axioelectric coupling of axion-like particles and the kinetic mixing parameter of dark photons, assuming the respective species constitutes all of the galactic dark matter. This study covers the mass range from 40 eV/$c^2$ to 500 eV/$c^2$ for both candidates, excluding previously untested parameter space for masses below ~1 keV/$c^2$. For the kinetic mixing of dark photons, values below $10^{-15}$ are reached for particle masses around 100 eV/$c^2$; for the axioelectric coupling of axion-like particles, values below $10^{-12}$ are reached for particles with masses in the range of a few-hundred eV/$c^2$.

preprint2021arXiv

Constraints on low-mass, relic dark matter candidates from a surface-operated SuperCDMS single-charge sensitive detector

This article presents an analysis and the resulting limits on light dark matter inelastically scattering off of electrons, and on dark photon and axion-like particle absorption, using a second-generation SuperCDMS high-voltage eV-resolution detector. The 0.93 gram Si detector achieved a 3 eV phonon energy resolution; for a detector bias of 100 V, this corresponds to a charge resolution of 3% of a single electron-hole pair. The energy spectrum is reported from a blind analysis with 1.2 gram-days of exposure acquired in an above-ground laboratory. With charge carrier trapping and impact ionization effects incorporated into the dark matter signal models, the dark matter-electron cross section $\barσ_{e}$ is constrained for dark matter masses from 0.5--$10^{4} $MeV$/c^{2}$; in the mass range from 1.2--50 eV$/c^{2}$ the dark photon kinetic mixing parameter $\varepsilon$ and the axioelectric coupling constant $g_{ae}$ are constrained. The minimum 90% confidence-level upper limits within the above mentioned mass ranges are $\barσ_{e}\,=\,8.7\times10^{-34}$ cm$^{2}$, $\varepsilon\,=\,3.3\times10^{-14}$, and $g_{ae}\,=\,1.0\times10^{-9}$.

preprint2021arXiv

Limits on astrophysical antineutrinos with the KamLAND experiment

We report on a search for electron antineutrinos ($\barν_e$) from astrophysical sources in the neutrino energy range 8.3 to 30.8 MeV with the KamLAND detector. In an exposure of 6.72 kton-year of the liquid scintillator, we observe 18 candidate events via the inverse beta decay reaction. Although there is a large background uncertainty from neutral current atmospheric neutrino interactions, we find no significant excess over background model predictions. Assuming several supernova relic neutrino spectra, we give upper flux limits of 60--110 cm$^{-2}$ s$^{-1}$ (90% CL) in the analysis range and present a model-independent flux. We also set limits on the annihilation rates for light dark matter pairs to neutrino pairs. These data improves on the upper probability limit of $^{8}$B solar neutrinos converting into $\barν_e$'s, $P_{ν_e \rightarrow \barν_e} < 3.5\times10^{-5}$ (90% CL) assuming an undistorted $\barν_e$ shape. This corresponds to a solar $\barν_e$ flux of 60 cm$^{-2}$ s$^{-1}$ (90% CL) in the analysis energy range.

preprint2021arXiv

Search for Solar Flare Neutrinos with the KamLAND detector

We report the result of a search for neutrinos in coincidence with solar flares from the GOES flare database. The search was performed on a 10.8 kton-year exposure of KamLAND collected from 2002 to 2019. This large exposure allows us to explore previously unconstrained parameter space for solar flare neutrinos. We found no statistical excess of neutrinos and established 90% confidence level upper limits of $8.4 \times 10^7$ cm$^{-2}$ ($3.0 \times 10^{9}$ cm$^{-2}$) on electron anti-neutrino (electron neutrino) fluence at 20 MeV normalized to the X12 flare, assuming that the neutrino fluence is proportional to the X-ray intensity.

preprint2021arXiv

Ultrafast Hole Spin Qubit with Gate-Tunable Spin-Orbit Switch

A key challenge in quantum computation is the implementation of fast and local qubit control while simultaneously maintaining coherence. Qubits based on hole spins offer, through their strong spin-orbit interaction, a way to implement fast quantum gates. Strikingly, for hole spins in one-dimensional germanium and silicon devices, the spin-orbit interaction has been predicted to be exceptionally strong yet highly tunable with gate voltages. Such electrical control would make it possible to switch on demand between qubit idling and manipulation modes. Here, we demonstrate ultrafast and universal quantum control of a hole spin qubit in a germanium/silicon core/shell nanowire, with Rabi frequencies of several hundreds of megahertz, corresponding to spin-flipping times as short as ~1 ns - a new record for a single-spin qubit. Next, we show a large degree of electrical control over the Rabi frequency, Zeeman energy, and coherence time - thus implementing a switch toggling from a rapid qubit manipulation mode to a more coherent idling mode. We identify an exceptionally strong but gate-tunable spin-orbit interaction as the underlying mechanism, with a short associated spin-orbit length that can be tuned over a large range down to 3 nm for holes of heavy-hole mass. Our work demonstrates a spin-orbit qubit switch and establishes hole spin qubits defined in one-dimensional germanium/silicon nanostructures as a fast and highly tunable platform for quantum computation.

preprint2020arXiv

GW190425: Observation of a Compact Binary Coalescence with Total Mass $\sim 3.4 M_{\odot}$

On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from 1.12 to 2.52 $M_{\odot}$ (1.45 to 1.88 $M_{\odot}$ if we restrict the dimensionless component spin magnitudes to be smaller than 0.05). These mass parameters are consistent with the individual binary components being neutron stars. However, both the source-frame chirp mass $1.44^{+0.02}_{-0.02} M_{\odot}$ and the total mass $3.4^{+0.3}_{-0.1}\,M_{\odot}$ of this system are significantly larger than those of any other known binary neutron star system. The possibility that one or both binary components of the system are black holes cannot be ruled out from gravitational-wave data. We discuss possible origins of the system based on its inconsistency with the known Galactic binary neutron star population. Under the assumption that the signal was produced by a binary neutron star coalescence, the local rate of neutron star mergers is updated to $250-2810 \text{Gpc}^{-3}\text{yr}^{-1}$.

preprint2020arXiv

Measurement of neutron-proton capture in the SNO+ water phase

The SNO+ experiment collected data as a low-threshold water Cherenkov detector from September 2017 to July 2019. Measurements of the 2.2-MeV $γ$ produced by neutron capture on hydrogen have been made using an Am-Be calibration source, for which a large fraction of emitted neutrons are produced simultaneously with a 4.4-MeV $γ$. Analysis of the delayed coincidence between the 4.4-MeV $γ$ and the 2.2-MeV capture $γ$ revealed a neutron detection efficiency that is centered around 50% and varies at the level of 1% across the inner region of the detector, which to our knowledge is the highest efficiency achieved among pure water Cherenkov detectors. In addition, the neutron capture time constant was measured and converted to a thermal neutron-proton capture cross section of $336.3^{+1.2}_{-1.5}$ mb.

preprint2020arXiv

Model comparison from LIGO-Virgo data on GW170817's binary components and consequences for the merger remnant

GW170817 is the very first observation of gravitational waves originating from the coalescence of two compact objects in the mass range of neutron stars, accompanied by electromagnetic counterparts, and offers an opportunity to directly probe the internal structure of neutron stars. We perform Bayesian model selection on a wide range of theoretical predictions for the neutron star equation of state. For the binary neutron star hypothesis, we find that we cannot rule out the majority of theoretical models considered. In addition, the gravitational-wave data alone does not rule out the possibility that one or both objects were low-mass black holes. We discuss the possible outcomes in the case of a binary neutron star merger, finding that all scenarios from prompt collapse to long-lived or even stable remnants are possible. For long-lived remnants, we place an upper limit of 1.9 kHz on the rotation rate. If a black hole was formed any time after merger and the coalescing stars were slowly rotating, then the maximum baryonic mass of non-rotating neutron stars is at most 3.05 $M_\odot$, and three equations of state considered here can be ruled out. We obtain a tighter limit of 2.67 $M_\odot$ for the case that the merger results in a hypermassive neutron star.

preprint2020arXiv

Nucleon effective mass in hot dense matter

Nucleon effective masses are studied in the framework of the Brueckner-Hartree-Fock many-body approach at finite temperature. Self-consistent calculations using the Argonne $V_{18}$ interaction including microscopic three-body forces are reported for varying temperature and proton fraction up to several times the nuclear saturation density. Our calculations are based on the exact treatment of the center-of-mass momentum instead of the average-momentum approximation employed in previous works. We discuss in detail the effects of the temperature together with those of the three-body forces, the density, and the isospin asymmetry. We also provide an analytical fit of the effective mass taking these dependencies into account. The temperature effects on the cooling of neutron stars are briefly discussed based on the results for betastable matter.

preprint2020arXiv

Search for Eccentric Binary Black Hole Mergers with Advanced LIGO and Advanced Virgo during their First and Second Observing Runs

When formed through dynamical interactions, stellar-mass binary black holes may retain eccentric orbits ($e>0.1$ at 10 Hz) detectable by ground-based gravitational-wave detectors. Eccentricity can therefore be used to differentiate dynamically-formed binaries from isolated binary black hole mergers. Current template-based gravitational-wave searches do not use waveform models associated to eccentric orbits, rendering the search less efficient to eccentric binary systems. Here we present results of a search for binary black hole mergers that inspiral in eccentric orbits using data from the first and second observing runs (O1 and O2) of Advanced LIGO and Advanced Virgo. The search uses minimal assumptions on the morphology of the transient gravitational waveform. We show that it is sensitive to binary mergers with a detection range that is weakly dependent on eccentricity for all bound systems. Our search did not identify any new binary merger candidates. We interpret these results in light of eccentric binary formation models.

preprint2020arXiv

Search for intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network

Gravitational wave astronomy has been firmly established with the detection of gravitational waves from the merger of ten stellar mass binary black holes and a neutron star binary. This paper reports on the all-sky search for gravitational waves from intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network. The search uses three independent algorithms: two based on matched filtering of the data with waveform templates of gravitational wave signals from compact binaries, and a third, model-independent algorithm that employs no signal model for the incoming signal. No intermediate mass black hole binary event was detected in this search. Consequently, we place upper limits on the merger rate density for a family of intermediate mass black hole binaries. In particular, we choose sources with total masses $M=m_1+m_2\in[120,800]$M$_\odot$ and mass ratios $q = m_2/m_1 \in[0.1,1.0]$. For the first time, this calculation is done using numerical relativity waveforms (which include higher modes) as models of the real emitted signal. We place a most stringent upper limit of $0.20$~Gpc$^{-3}$yr$^{-1}$ (in co-moving units at the 90% confidence level) for equal-mass binaries with individual masses $m_{1,2}=100$M$_\odot$ and dimensionless spins $χ_{1,2}= 0.8$ aligned with the orbital angular momentum of the binary. This improves by a factor of $\sim 5$ that reported after Advanced LIGO's first observing run.

preprint2020arXiv

Strong spin-orbit interaction and $g$-factor renormalization of hole spins in Ge/Si nanowire quantum dots

The spin-orbit interaction lies at the heart of quantum computation with spin qubits, research on topologically non-trivial states, and various applications in spintronics. Hole spins in Ge/Si core/shell nanowires experience a spin-orbit interaction that has been predicted to be both strong and electrically tunable, making them a particularly promising platform for research in these fields. We experimentally determine the strength of spin-orbit interaction of hole spins confined to a double quantum dot in a Ge/Si nanowire by measuring spin-mixing transitions inside a regime of spin-blockaded transport. We find a remarkably short spin-orbit length of $\sim$65 nm, comparable to the quantum dot length and the interdot distance. We additionally observe a large orbital effect of the applied magnetic field on the hole states, resulting in a large magnetic field dependence of the spin-mixing transition energies. Strikingly, together with these orbital effects, the strong spin-orbit interaction causes a significant enhancement of the $g$-factor with magnetic field.The large spin-orbit interaction strength demonstrated is consistent with the predicted direct Rashba spin-orbit interaction in this material system and is expected to enable ultrafast Rabi oscillations of spin qubits and efficient qubit-qubit interactions, as well as provide a platform suitable for studying Majorana zero modes.

preprint2019arXiv

An Optically Targeted Search for Gravitational Waves emitted by Core-Collapse Supernovae during the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed within a source distance of approximately 20 Mpc during the first and second observing runs of Advanced LIGO and Advanced Virgo. No significant gravitational-wave candidate was detected. We report the detection efficiencies as a function of the distance for waveforms derived from multidimensional numerical simulations and phenomenological extreme emission models. For neutrino-driven explosions the distance at which we reach 50% detection efficiency is approaching 5 kpc, and for magnetorotationally-driven explosions is up to 54 kpc. However, waveforms for extreme emission models are detectable up to 28 Mpc. For the first time, the gravitational-wave data enabled us to exclude part of the parameter spaces of two extreme emission models with confidence up to 83%, limited by coincident data coverage. Besides, using ad hoc harmonic signals windowed with Gaussian envelopes we constrained the gravitational-wave energy emitted during core-collapse at the levels of $4.27\times 10^{-4}\,M_\odot c^2$ and $1.28\times 10^{-1}\,M_\odot c^2$ for emissions at 235 Hz and 1304 Hz respectively. These constraints are two orders of magnitude more stringent than previously derived in the corresponding analysis using initial LIGO, initial Virgo and GEO 600 data.

preprint2019arXiv

Search for sub-solar mass ultracompact binaries in Advanced LIGO's second observing run

We present an Advanced LIGO and Advanced Virgo search for sub-solar mass ultracompact objects in data obtained during Advanced LIGO's second observing run. In contrast to a previous search of Advanced LIGO data from the first observing run, this search includes the effects of component spin on the gravitational waveform. We identify no viable gravitational wave candidates consistent with sub-solar mass ultracompact binaries with at least one component between 0.2 - 1.0 solar masses. We use the null result to constrain the binary merger rate of (0.2 solar mass, 0.2 solar mass) binaries to be less than 3.7 x 10^5 Gpc^-3 yr^-1 and the binary merger rate of (1.0 solar mass, 1.0 solar mass) binaries to be less than 5.2 x 10^3 Gpc^-3 yr^-1. Sub-solar mass ultracompact objects are not expected to form via known stellar evolution channels, though it has been suggested that primordial density fluctuations or particle dark matter with cooling mechanisms and/or nuclear interactions could form black holes with sub-solar masses. Assuming a particular primordial black hole formation model, we constrain a population of merging 0.2 solar mass black holes to account for less than 16% of the dark matter density and a population of merging 1.0 solar mass black holes to account for less than 2% of the dark matter density. We discuss how constraints on the merger rate and dark matter fraction may be extended to arbitrary black hole population models that predict sub-solar mass binaries.

preprint2016arXiv

High voltage testing for the MAJORANA Demonstrator

The MAJORANA Collaboration is constructing the MAJORANA Demonstrator, an ultra-low background, 44-kg modular high-purity Ge (HPGe) detector array to search for neutrinoless double-beta decay in Ge-76. The phenomenon of surface micro-discharge induced by high-voltage has been studied in the context of the MAJORANA Demonstrator. This effect can damage the front-end electronics or mimic detector signals. To ensure the correct performance, every high-voltage cable and feedthrough must be capable of supplying HPGe detector operating voltages as high as 5 kV without exhibiting discharge. R&D measurements were carried out to understand the testing system and determine the optimum design configuration of the high-voltage path, including different improvements of the cable layout and feedthrough flange model selection. Every cable and feedthrough to be used at the MAJORANA Demonstrator was characterized and the micro-discharge effects during the MAJORANA Demonstrator commissioning phase were studied. A stable configuration has been achieved, and the cables and connectors can supply HPGe detector operating voltages without exhibiting discharge.

preprint2016arXiv

Nonrelativistic nucleon effective masses in nuclear matter: BHF versus RHF

The density and isospin dependences of the nonrelativistic nucleon effective mass ($m^*$) are studied, which is a measure of the nonlocality of the single particle (s.p.) potential. We decouple it further into the so called k-mass ($m^*_k$, i.e., the nonlocality in space) and E-mass ($m^*_E$, i.e., the nonlocality in time). Both masses are determined and compared from the latest versions of the nonrelativistic Brueckner-Hartree Fock (BHF) model and the relativistic Hartree-Fock (RHF) model. The latter are achieved based on the corresponding Schrödinger equivalent s.p. potential in a relativistic framework. We demonstrate the origins of different effective masses and discuss also their neutron-proton splitting in the asymmetric matter in different models. We find that the neutron-proton splittings of both the k-mass and the E-mass have the same asymmetry dependences at considered densities, namely $m^*_{k,n} > m^*_{k,p}$ and $m^*_{E,p} > m^*_{E,n}$. However, the resulting splittings of nucleon effective masses could have different asymmetry dependences in the two models, because they could be dominated either by that of the k-mass (then we have $m^*_n > m^*_p$ in the BHF model) or by that of the E-mass (then we have $m^*_p > m^*_n$ in the RHF model).

preprint2016arXiv

Tip Pressure Induced Incoherent Energy Gap in CaFe2As2

In CaFe2As2, superconductivity can be achieved by applying a modest c-axis pressure of several kbar. Here we use scanning tunneling microscopy/spectroscopy (STM/S) to explore the STM tip pressure effect on single crystals of CaFe2As2. When performing STM/S measurements, the tip-sample interaction can be controlled to act repulsive with reduction of the junction resistance, thus to apply a tip pressure on the sample. We find that an incoherent energy gap emerges at the Fermi level in the differential conductance spectrum when the tip pressure is increased. This energy gap is of the similar order of magnitude as the superconducting gap in the chemical doped compound Ca0.4Na0.6Fe2As2 and disappears at the temperature well below that of the bulk magnetic ordering. Moreover, we also observe the rhombic distortion of the As lattice, which agrees with the orthorhombic distortion of the underlying Fe lattice. These findings suggest that the STM tip pressure can induce the local Cooper pairing in the orthorhombic phase of CaFe2As2.

preprint2015arXiv

Early Recognition of Emerging Flu Strain Clusters

Minimizing time delays in manufacturing vaccines appropriate to rapidly mutating viruses is the key step for improving vaccine effectiveness. The vaccine for the H3N2 flu type has failed for the last two years (~ 15% effective). Here we summarize the state of the predictive art and report the most current results for H3N2 flu vaccine design. Using a 2006 model of dimensional reduction of viral mutational complexity, we show that this model can reduce vaccine time delays by a year or more in some cases.

preprint2015arXiv

Massive hybrid stars with a first order phase transition

We develop our previous study of the transition to deconfined quark phase in neutron stars, including the interaction in the quark equation of state to the leading order in the perturbative expansion within the confinement density-dependent mass model. Using the Gibbs conditions the hadron-quark mixed phase is constructed matching the latter with the hadron equation of state derived from the microscopic Brueckner-Hartree-Fock approximation. The influence of quark interaction parameters on threshold properties and phase diagram of dense neutron star matter are discussed in detail. We find that the leading-order quark interaction expands the density range of the mixed phase, pushing forward the disappearance of the hadron phase. Moreover, since the equation of state could turn out to be stiffer, a high-mass hybrid star is possible with mixed-phase core with typical parameter sets.

preprint2015arXiv

Structures of the Vela pulsar and the glitch crisis from the Brueckner theory

Detailed structures of the Vela pulsar (PSR B0833-45, with a period of $89.33$ milliseconds) are predicted by adopting a recently-constructed unified treatment of all parts of neutron stars: the outer crust, the inner crust and the core based on modern microscopic Brueckner-Hartree-Fock calculations. To take the pulsar mass ranging from $1.0M_{\odot}$ to $2.0M_{\odot}$, we calculate the central density, the core/crust radii, the core/crust mass, the core/crustal thickness, the moment of inertia, and the crustal moment of inertia. Among them, the crustal moment of inertia could be effectively constrained from the accumulated glitch observations, which has been a great debate recently, known as "glitch crisis". Namely, superfluid neutrons contained in the inner crust, which are regarded as the origin of the glitch in the standard two-component model, could be largely entrained in the nuclei lattices, then there may not be enough superfluid neutrons ($\sim 4/5$ less than the previous value) to trigger the large glitches ($Δν/ν_0 \sim 10^{-6}$) in the Vela pulsar. We then provide the first analysis of the crisis based on the microscopic basis, by confronting the glitch observations with the theoretical calculations for the crustal moment of inertia. We find that despite some recent opposition to the crisis argument, the glitch crisis is still present, which means that besides the crust superfluid neutrons, core neutrons might be necessary for explaining the large glitches of the Vela pulsar.

preprint2014arXiv

Charmonium, $D_s$ and $D_s^*$ from overlap fermion on DWF configurations

We take a new approach to determine the scale parameter $r_0$, the physical masses of strange and charm quarks through a global fit which incorporates continuum extrapolation, chiral extrapolation and quark mass interpolation to the lattice data. The charmonium and charm-strange meson spectrum are calculated with overlap valence quarks on $2+1$-flavor domain-wall fermion gauge configurations generated by the RBC and UKQCD Collaboration. We use the masses of $D_s$, $D_s^*$ and $J/ψ$ as inputs and obtain $m_c^{\overline{\rm MS}}(2\,{\rm GeV})=1.110(24)\,{\rm GeV}$, $m_s^{\overline{\rm MS}}(2\,{\rm GeV})=0.104(9)\,{\rm GeV}$ and $r_0=0.458(11)\,{\rm fm}$. Subsequently, the hyperfine-splitting of charmonium and $f_{D_s}$ are predicted to be $112(5)\,{\rm MeV}$ and $254(5)\,{\rm MeV}$, respectively.

preprint2014arXiv

Extended quark mean-field model for neutron stars

We extend the quark mean-field (QMF) model to strangeness freedom to study the properties of hyperons ($Λ,Σ,Ξ$) in infinite baryon matter and neutron star properties. The baryon-scalar meson couplings in the QMF model are determined self-consistently from the quark level, where the quark confinement is taken into account in terms of a scalar-vector harmonic oscillator potential. The strength of such confinement potential for $u,d$ quarks is constrained by the properties of finite nuclei, while the one for $s$ quark is limited by the properties of nuclei with a $Λ$ hyperon. These two strengths are not same, which represents the SU(3) symmetry breaking effectively in the QMF model. Also, we use an enhanced $Σ$ coupling with the vector meson, and both $Σ$ and $Ξ$ hyperon potentials can be properly described in the model. The effects of the SU(3) symmetry breaking on the neutron star structures are then studied. We find that the SU(3) breaking shifts earlier the hyperon onset density and makes hyperons more abundant in the star, in comparisons with the results of the SU(3) symmetry case. However, it does not affect much the star's maximum mass. The maximum masses are found to be $1.62 M_{\odot}$ with hyperons and $1.88 M_{\odot}$ without hyperons. The present neutron star model is shown to have limitations on explaining the recently measured heavy pulsar.

preprint2014arXiv

Searching for Debris Disks around Seven Radio Pulsars

We report on our searches for debris disks around seven relatively nearby radio pulsars, which are isolated sources and were carefully selected as the targets on the basis of our deep $K_s$-band imaging survey. The $K_s$ images obtained with the 6.5\,m Baade Magellan Telescope at Las Campanas Observatory are analyzed together with the \textit{Spitzer}/IRAC images at 4.5 and 8.0~$μ$m and the \textit{WISE} images at 3.4, 4.6, 12 and 22~$μ$m. No infrared (IR) counterparts of these pulsars are found, with flux upper limits of $\sim μ$Jy at near-infrared ($λ<10\ μ$m) and $\sim$10--1000\,$μ$Jy at mid-infrared wavelengths ($λ>10\ μ$m). The results of this search are discussed in terms of the efficiency of converting the pulsar spin-down energy to thermal energy and X-ray heating of debris disks, with comparison made to the two magnetars 4U~0142+61 and 1E~2259+586 which are suggested to harbor a debris disk.

preprint2013arXiv

Dark matter effect on the mass measurement of neutron stars

Newly-determined mass of 1.97 $\pm$ 0.04 $M_{\odot}$ for PSR J1614-2230 has been a challenge for the neutron star with a hyperon core (namely hyperon star), since hyperons usually reduce the theoretical maximum mass of the star. In this article, we consider dark matter as another possible constituent in hyperon stars' interior to loose this mass constrain. We take dark matter as self-interacting Fermi gas with certain repulsive interaction among the dark matter particles and non-interaction between dark matter and ordinary matter as is generally assumed. We find that the star maximum mass is sensitive to the particle mass of dark matter, and a high enough star mass larger than 2 $M_{\odot}$ could be achieved when the particle mass is small enough. In this particular model, a strong upper limit 0.64 GeV for dark matter mass is obtained in strongly-interacting dark matter and 0.16 GeV for dark matter mass in weakly-interacting dark matter. Dark matter accumulated around the star could also contribute to the mass measurement, however, such contribution could be safely ignored when the generally used dark matter density is assumed.

preprint2013arXiv

Experimental investigation of the electronic structure of Ca$_{0.83}$La$_{0.17}$Fe$_2$As$_2$

We performed a combined angle-resolved photoemission spectroscopy and scanning tunneling microscopy study of the electronic structure of electron-doped Ca$_{0.83}$La$_{0.17}$Fe$_2$As$_2$. A surface reconstruction associated with the dimerization of As atoms is observed directly in the real space, as well as the consequent band folding in the momentum space. Besides this band folding effect, the Fermi surface topology of this material is similar to that reported previously for BaFe$_{1.85}$Co$_{0.15}$As$_2$, with $Γ$-centred hole pockets quasi-nested to M-centred electron pockets by the antiferromagnetic wave vector. Although no superconducting gap is observed by ARPES possibly due to low superconducting volume fraction, a gap-like density of states depression of $7.7\pm 2.9$ meV is determined by scanning tunneling microscopy.

preprint2013arXiv

Quark mean-field model for single and double $Λ$ and $Ξ$ hypernuclei

The study of hypernuclei in the quark mean-field model are extended from single-$Λ$ hypernuclei to double-$Λ$ hypernuclei as well as $Ξ$ hypernuclei (double strangeness nuclei), with a broken SU(3) symmetry for the quark confinement potential. The strength of the potential for $u,d$ quarks is constrained by the properties of finite nuclei, while the one for $s$ quark is fixed by the single $Λ$-nucleus potential at the nuclear saturation density, which has a slightly different value. Compared to our previous work, we find that the introduction of such symmetry breaking improves efficiently the description of the single $Λ$ energies for a wide range of mass region, which demonstrates the importance of this effect on the hypernuclei study. Predictions for double-$Λ$ hypernuclei and $Ξ$ hypernuclei are then made, confronting with a few available experiential double-$Λ$ hypernuclei data and other model calculations. The consistency of our results with the available double-$Λ$ hypernucleus is found to be surprisingly good. Also, there is generally a bound state for $Ξ^0$ hypernuclei with both light and heavy core nuclei in our model.

preprint2013arXiv

Revisiting the hot matter in the center of gamma-ray bursts and supernova

Hot matter with nucleons can be produced in the inner region of the neutrino-dominated accretion flow in gamma-ray bursts or during the proto-neutron star birth in successful supernova. The composition and equation of state of the matter depend on the dynamic $β$ equilibrium under various neutrino opacities. The strong interaction between nucleons may also play an important role. We plan to extend the previous studies by incorporating these two aspects in our model. The modification of the $β$-equilibrium condition from neutrino optically thin to thick has been modeled by an equilibrium factor $χ$ ranging between the neutrino-freely-escaping case and the neutrino-trapped case. We employ the microscopic Brueckner-Hartree-Fock approach extended to the finite temperature regime to study the interacting nucleons. We show that the composition and chemical potentials of the hot nuclear matter for different densities and temperatures at each stage of $β$ equilibrium. We also compare our realistic equation of states with those of the free gas model. We find the neutrino opacity and the strong interaction between nucleons are important for the description and should be taken into account in model calculations.

preprint2013arXiv

Strangeness and charmness content of nucleon from overlap fermions on 2+1-flavor domain-wall fermion configurations

We present a calculation of the strangeness and charmness contents <N|\bar{s}s|N> and <N|\bar{c}c|N> of the nucleon from dynamical lattice QCD with 2+1 flavors. The calculation is performed with overlap valence quarks on 2+1-flavor domain-wall fermion gauge configurations. The configurations are generated by the RBC collaboration on a 24^3*64 lattice with sea quark mass am_l=0.005, am_s=0.04, and inverse lattice spacing a^{-1}=1.73GeV. Both actions have chiral symmetry which is essential in avoiding contamination due to the operator mixing with other flavors. Nucleon propagator and the quark loops are both computed with stochastic grid sources, while low-mode substitution and low-mode averaging methods are used respectively which substantially improve the signal to noise ratio. We obtain the strangeness matrix element f_{T_{s}} = m_s <N|\bar{s}s|N> / M_N = 0.0334(62), and the charmness content f_{T_{c}} = m_c <N|\bar{c}c|N> / M_N = 0.094(31) which is resolved from zero by 3σprecision for the first time.

preprint2013arXiv

Tensor correlation, pairing interaction and deformation in Ne isotopes and Ne hypernuclei

We study tensor and pairing effects on the quadruple deformation of neon isotopes based on a deformed Skyrme-Hartree-Fock model with BCS approximation for the pairing channel. We extend the Skyrme-Hartree-Fock formalism for the description of single- and double-lambda hypernuclei adopting two different hyperon-nucleon interactions. It is found that the interplay of pairing and tensor interactions is crucial to derive the deformations in several neon isotopes. Especially, the shapes of $^{26,30}$Ne are studied in details in comparisons with experimentally observed shapes. Furthermore the deformations of the hypernuclei are compared with the corresponding neon isotopic cores in the presence of tensor force. We find the same shapes with somewhat smaller deformations for single $Λ$-hypernuclei compared with their core deformations. It is also pointed out that the latest version of hyperon interaction, the ESC08b model, having a deeper $Λ$ potential makes smaller deformations for hypernuclei than those of another NSC97f model.

preprint2013arXiv

Tensor force and shape evolution of Si isotopes in Skyrme-Hartree-Fock model

Much interest has been devoted to the shape evolution of neutron-rich Si isotopes recently both experimentally and theoretically. We provide our study of $^{28-42}$Si by Skyrme-Hartree-Fock model with BCS approximation for the pairing channel. We use an empirical pairing parameter deduced from the experimental binding energy data for each nuclei and each Skyrme parametrization. The recent highlight of the deformation change from a prolate shape of $^{38}$Si to an oblate shape of $^{42}$Si is confirmed in the present model. We also emphasize the role of tensor force on deformations of neutron-rich $^{30,32}$Si, and discuss its underlying physics.

preprint2013arXiv

The thermal dust emission in the N158-N159-N160 (LMC) star forming complex mapped by Spitzer, Herschel and LABOCA

We present a study of the infrared/submm emission of the LMC star forming complex N158-N159-N160. Combining observations from the Spitzer Space Telescope (3.6-70um), the Herschel Space Observatory (100-500um) and LABOCA (870um) allows us to work at the best angular resolution available now for an extragalactic source. We observe a remarkably good correlation between SPIRE and LABOCA emission and resolve the low surface brightnesses emission. We use the Spitzer and Herschel data to perform a resolved Spectral Energy Distribution (SED) modelling of the complex. Using MBB, we derive a global emissivity index beta_c of 1.47. If beta cold is fixed to 1.5, we find an average temperature of 27K. We also apply the Galliano et al. (2011) modelling technique (and amorphous carbon to model carbon dust) to derive maps of the star formation rate, the mean starlight intensity, the fraction of PAHs or the dust mass surface density of the region. We observe that the PAH fraction strongly decreases in the HII regions. This decrease coincides with peaks in the mean radiation field intensity map. The dust surface densities follow the FIR distribution, with a total dust mass of 2.1x10^4 Msolar (2.8 times less than when using graphite grains) in the resolved elements we model. We find a non-negligible amount of dust in the molecular cloud N159 South (showing no massive SF). We also investigate the drivers of the Herschel/PACS and SPIRE submm colours as well as the variations in the gas-to-dust mass ratio (G/D) and the XCO conversion factor in the region N159. We finally model individual regions to analyse variations in the SED shape across the complex and the 870um emission in more details. No measurable submm excess emission at 870um seems to be detected in these regions.

preprint2012arXiv

The $Δ_{mix}$ parameter in the overlap on domain-wall mixed action

A direct calculation of the mixed-action parameter $Δ_{mix}$ with valence overlap fermions on a domain-wall fermion sea is presented. The calculation is performed on four ensembles of the 2+1-flavor domain-wall gauge configurations: $24^3 \times 64$ ($a m_l= 0.005$, $a=0.114\fm$) and $32^3 \times 64$ ($a m_l = 0.004, 0.006, 0.008$, $a=0.085\fm$). For pion masses close to $300\MeV$ we find \hbox{$Δ_{mix}=0.030(6)\GeV^4$} at $a=0.114\fm$ and $Δ_{mix}=0.033(12)\GeV^4$ at $a=0.085\fm$. The results are quite independent of the lattice spacing and they are significantly smaller than the results for valence domain-wall fermions on Asqtad sea or those of valence overlap fermions on clover sea. Combining the results extracted from these two ensembles, we get $Δ_{mix}=0.030(6)(5)\GeV^4$, where the first error is statistical and the second is the systematic error associated with the fitting method.

preprint2012arXiv

The Strangeness and Charmness of Nucleon from Overlap Fermions

The calculation of the strangeness and charmness of the nucleon is presented with overlap fermion action on 2+1 flavor domain wall fermion configurations. We adopt stochastic grid sources and the low mode substitution technique to improve the signals of nucleon correlation functions and the loops. The calculation is done on a $24^3\times 64$ lattice with $m_l=0.005$, $m_h=0.04$, and $a^{-1}=1.73\,{\rm GeV}$. We find $ f_{T_{s}} = 0.048(15)$ and $f_{T_{c}} = 0.029(43)$.

preprint2011arXiv

Dust Emission from Evolved and Unevolved HII Regions in the Large Magellanic Cloud

We present a study of the dust properties of 12 classical and superbubble HII regions in the Large Magellanic Cloud. We use infrared photometry from Spitzer (8, 24, 70, and 160 \mum bands), obtained as part of the Surveying the Agents of a Galaxy's Evolution (SAGE) program, along with archival spectroscopic classifications of the ionizing stars to examine the role of stellar sources on dust heating and processing. Our infrared observations show surprisingly little correlation between the emission properties of the dust and the effective temperatures or bolometric magnitudes of stars in the HII regions, suggesting that the HII region evolutionary timescale is not on the order of the dust processing timescale. We find that the infrared emission of superbubbles and classical HII regions shows little differentiation between the two classes, despite the significant differences in age and morphology. We do detect a correlation of the 24 \mum emission from hot dust with the ratio of 70 to 160 \mum flux. This correlation can be modeled as a trend in the temperature of a minority hot dust component, while a majority of the dust remains significantly cooler.

preprint2011arXiv

Non-Standard Grain Properties, Dark Gas Reservoir, and Extended Submillimeter Excess, Probed by Herschel in the Large Magellanic Cloud

Aims: In this paper, we perform detailed modelling of the Spitzer and Herschel observations of the LMC, in order to: (i) systematically study the uncertainties and biases affecting dust mass estimates; and to (ii) explore the peculiar ISM properties of the LMC. Methods: To achieve these goals, we have modelled the spatially resolved SEDs with two alternate grain compositions, to study the impact of different submillimetre opacities on the dust mass. We have rigorously propagated the observational errors (noise and calibration) through the entire fitting process, in order to derive consistent parameter uncertainties. Results: First, we show that using the integrated SED leads to underestimating the dust mass by ~50 % compared to the value obtained with sufficient spatial resolution, for the region we studied. This might be the case, in general, for unresolved galaxies. Second, we show that Milky Way type grains produce higher gas-to-dust mass ratios than what seems possible according to the element abundances in the LMC. A spatial analysis shows that this dilemma is the result of an exceptional property: the grains of the LMC have on average a larger intrinsic submm opacity (emissivity index beta~1.7 and opacity kappa_abs(160 microns)=1.6 m2/kg) than those of the Galaxy. By studying the spatial distribution of the gas-to-dust mass ratio, we are able to constrain the fraction of unseen gas mass between ~10, and ~100 % and show that it is not sufficient to explain the gas-to-dust mass ratio obtained with Milky Way type grains. Finally, we confirm the detection of a 500 microns extended emission excess with an average relative amplitude of ~15 %, varying up to 40 %. This excess anticorrelates well with the dust mass surface density. Although we do not know the origin of this excess, we show that it is unlikely the result of very cold dust, or CMB fluctuations.

preprint2011arXiv

Study of the scalar charmed-strange meson Ds0*(2317) with chiral fermions

The recently discovered charmed-strangemeson Ds0*(2317) has been speculated to be a tetraquark mesonium. We study this suggestion with overlap fermions on 2+1 flavor domain wall fermion configurations. We use 4-quark interpolating operatorswith Z4 grid sources on two lattices (16\times16\times16\times32 and 24\times24\times24 times64) to study the volume dependence of the states in an attempt to discern the nature of the states in the four-quark correlator to see if they are all two-meson scattering states or if one is a tetraquark mesonium. We also use the hybrid boundary condition method for this purpose which is designed to lift the two-meson states in energy while leaving the tetraquark mesonium unchanged. We find that the volume method is not effective in the present case due to the fact that the scattering states spectrum is closely packed for such heavy states so that one cannot separate out individual scattering states since the volume dependence is skewed as a result. However, the hybrid boundary condition method works and we found that the four-quark correlators can be fitted with a tower of two-meson scattering states. We conclude that we do not see a tetraquark mesonium in the Ds0*(2317) meson region.

preprint2010arXiv

Cold dust clumps in dynamically hot gas

We present clumps of dust emission from Herschel observations of the Large Magellanic Cloud (LMC) and their physical and statistical properties. We catalog cloud features seen in the dust emission from Herschel observations of the LMC, the Magellanic type irregular galaxy closest to the Milky Way, and compare these features with HI catalogs from the ATCA+Parkes HI survey. Using an automated cloud-finding algorithm, we identify clouds and clumps of dust emission and examine the cumulative mass distribution of the detected dust clouds. The mass of cold dust is determined from physical parameters that we derive by performing spectral energy distribution fits to 250, 350, and 500 micronm emission from SPIRE observations using DUSTY and GRASIL radiative transfer calculation with dust grain size distributions for graphite/silicate in low-metallicity extragalactic environments. The dust cloud mass spectrum follows a power law distribution with an exponent of gamma=-1.8 for clumps larger than 400 solar mass and is similar to the HI mass distribution. This is expected from the theory of ISM structure in the vicinity of star formation.

preprint2010arXiv

Determining dust temperatures and masses in the Herschel era: the importance of observations longward of 200 micron

The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (>100 micron). Only with SPIRE on Herschel has it been possible to get high spatial resolution at 200 to 500 micron that is beyond the peak (~160 micron) of dust emission in most galaxies. We investigate the differences in the fitted dust temperatures and masses determined using only <200 micron data and then also including >200 micron data (new SPIRE observations) to determine how important having >200 micron data is for deriving these dust properties. We fit the 100 to 350 micron observations of the Large Magellanic Cloud (LMC) point-by-point with a model that consists of a single temperature and fixed emissivity law. The data used are existing observations at 100 and 160 micron (from IRAS and Spitzer) and new SPIRE observations of 1/4 of the LMC observed for the HERITAGE Key Project as part of the Herschel Science Demonstration phase. The dust temperatures and masses computed using only 100 and 160 micron data can differ by up to 10% and 36%, respectively, from those that also include the SPIRE 250 & 350 micron data. We find that an emissivity law proportional to lambda^-1.5 minimizes the 100-350 micron fractional residuals. We find that the emission at 500 micron is ~10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the fractional 500 micron excess is weakly anti-correlated with MIPS 24 micron flux and the total gas surface density. This argues against a flux calibration error as the origin of the 500 micron excess. Our results do not allow us to distinguish between a systematic variation in the wavelength dependent emissivity law or a population of very cold dust only detectable at lambda > 500 micron for the origin of the 500 micron excess.

preprint2010arXiv

Herschel Inventory of The Agents of Galaxy Evolution (HERITAGE): the Large Magellanic Cloud dust

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory's PACS and SPIRE photometry cameras, we imaged a 2x8 square degree strip through the LMC, at a position angle of ~22.5 degrees as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 μm and SPIRE 250, 350 and 500 μm. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 μm emission is in excess of the models by 6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_PAH) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J=1-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65(+15,-18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287(+25,-42) that is more consistent with a GDR inferred from extinction.

preprint2010arXiv

Meson spectra from overlap fermion on domain wall gauge configurations

We report meson spectra obtained by using valence overlap fermion propagators generated on a background of 2+1 flavor domain wall fermion gauge configurations on 16^3 X 32, 24^3 X 64 and 32^3 X 64 lattices. We use many-to-all correlators with Z3 grid source and low eigenmode substitution which is efficient in reducing errors for the hadron correlators. The preliminary results on meson spectrum, a0 correlators, and charmonium hyperfine splitting for three sea quark masses are reported here.

preprint2010arXiv

Overlap Valence on 2+1 Flavor Domain Wall Fermion Configurations with Deflation and Low-mode Substitution

The overlap fermion propagator is calculated on 2+1 flavor domain wall fermion gauge configurations on 16^3 x 32, 24^3 x 64 and 32^3 x 64 lattices. With HYP smearing and low eigenmode deflation, it is shown that the inversion of the overlap operator can be expedited by ~ 20 times for the 16^3 x 32 lattice and ~ 80 times for the 32^3 x 64 lattice. Through the study of hyperfine splitting, we found that the O(m^2a^2) error is small and these dynamical fermion lattices can adequately accommodate quark mass up to the charm quark. The low energy constant Δ_{mix} which characterizes the discretization error of the pion made up of a pair of sea and valence quarks in this mixed action approach is calculated via the scalar correlator with periodic and anti-periodic boundary conditions. It is found to be small which shifts a 300 MeV pion mass by ~ 10 to 19 MeV on these sets of lattices. We have studied the signal-to-noise issue of the noise source for the meson and baryon. It is found that the many-to-all meson and baryon correlators with Z_3 grid source and low eigenmode substitution is efficient in reducing errors for the correlators of both mesons and baryons. With 64-point Z_3 grid source and low-mode substitution, it can reduce the statistical errors of the light quark (m_π ~ 200 - 300 MeV) meson and nucleon correlators by a factor of ~ 3-4 as compared to the point source. The Z_3 grid source itself can reduce the errors of the charmonium correlators by a factor of ~ 3.

preprint2010arXiv

Protoneutron stars in the Brueckner-Hartree-Fock approach and finite-temperature kaon condensation

We study the properties of hot neutrino-trapped beta-stable stellar matter using an equation of state of nuclear matter within the Brueckner-Hartree-Fock approach including three-body forces, combined with a standard chiral model for kaon condensation at finite temperature. The properties of (proto)neutron stars are then investigated within this framework.

preprint2010arXiv

Strange stars with different quark mass scaling

We investigate the stability of strange quark matter and the properties of the cor- responding strange stars, within a wide range of quark mass scaling. The calculation shows that the resulting maximum mass always lies between 1.5M_{sun} and 1.8M_{sun} for all the scalings chosen here. Strange star sequences with a linear scaling would sup- port less gravitational mass, and a change (increase or decrease) of the scaling around the linear scaling would lead to a higher maximum mass. Radii invariably decrease with the mass scaling. Then the larger the scaling, the faster the star might spin. In addition, the variation of the scaling would cause an order of magnitude change of the strong electric field on quark surface, which is essential to support possible crusts of strange stars against gravity and may then have some astrophysical implications.

preprint2010arXiv

The Herschel revolution: unveiling the morphology of the high mass star formation sites N44 and N63 in the LMC

We study the structure of the medium surrounding sites of high-mass star formation to determine the interrelation between the HII regions and the environment from which they were formed. The density distribution of the surroundings is key in determining how the radiation of the newly formed stars interacts with the surrounds in a way that allows it to be used as a star formation tracer. We present new Herschel/SPIRE 250, 350 and 500 mum data of LHA 120-N44 and LHA 120-N63 in the LMC. We construct average spectral energy distributions (SEDs) for annuli centered on the IR bright part of the star formation sites. The annuli cover ~10-~100 pc. We use a phenomenological dust model to fit these SEDs to derive the dust column densities, characterise the incident radiation field and the abundance of polycyclic aromatic hydrocarbon molecules. We see a factor 5 decrease in the radiation field energy density as a function of radial distance around N63. N44 does not show a systematic trend. We construct a simple geometrical model to derive the 3-D density profile of the surroundings of these two regions. Herschel/SPIRE data have proven very efficient in deriving the dust mass distribution. We find that the radiation field in the two sources behaves very differently. N63 is more or less spherically symmetric and the average radiation field drops with distance. N44 shows no systematic decrease of the radiation intensity which is probably due to the inhomogeneity of the surrounding molecular material and to the complex distribution of several star forming clusters in the region.

preprint2009arXiv

Mid-Infrared Extinction and its Variation with Galactic Longitude

Based on the data obtained from the Spitzer/GLIPMSE Legacy Program and the 2MASS project, we derive the extinction in the four IRAC bands, [3.6], [4.5], [5.8] and [8.0] micron, relative to the 2MASS Ks band (at 2.16 micron) for 131 GLIPMSE fields along the Galactic plane within |l|<65 deg, using red giants and red clump giants as tracers. As a whole, the mean extinction in the IRAC bands (normalized to the 2MASS Ks band), A_[3.6]/A_Ks=0.63, A_[4.5]/A_Ks=0.57, A_[5.8]/A_Ks=0.49, A_[8.0]/A_Ks=0.55, exhibits little variation with wavelength (i.e. the extinction is somewhat flat or gray). This is consistent with previous studies and agrees with that predicted from the standard interstellar grain model for R_V=5.5 by Weingartner & Draine (2001). As far as individual sightline is concerned, however, the wavelength dependence of the mid-infrared interstellar extinction A_{lambda}/A_Ks varies from one sightline to another, suggesting that there may not exist a "universal" IR extinction law. We, for the first time, demonstrate the existence of systematic variations of extinction with Galactic longitude which appears to correlate with the locations of spiral arms as well as with the variation of the far infrared luminosity of interstellar dust.

preprint2009arXiv

On the Carriers of the 21 Micron Emission Feature in Post-Asymptotic Giant Branch Stars

The mysterious 21micron emission feature seen in sixteen C-rich proto-planetary nebulae (PPNe) remains unidentified since its discovery in 1989. Over a dozen of materials are suggested as the carrier candidates. In this work we quantitatively investigate eight inorganic and one organic carrier candidates in terms of elemental abundance constraints, while previous studies mostly focus on their spectral profiles (which could be largely affected by grain size, shape and clustering effects). It is found that: (1) five candidates (TiC nanoclusters, fullerenes coordinated with Ti atoms, SiS$_2$, doped-SiC, and SiO$_2$-coated SiC dust) violate the abundance constraints (i.e. they require too much Ti, S or Si to account for the emission power of the 21micron band, (2) three candidates (carbon and silicon mixtures, Fe$_2$O$_3$, and Fe$_3$O$_4$),while satisfying the abundance constraints, exhibit secondary features which are not detected in the 21micron sources, and (3) nano FeO, neither exceeding the abundance budget nor producing undetected secondary features, seems to be a viable candidate, supporting the suggestions of Posch et al. 2004.

preprint2009arXiv

The charmed-strange meson spectrum from overlap fermions on domain wall dynamical fermion configurations

The charmed-strange meson spectrum is calculated with the overlap valence fermions on 2+1 flavor domain wall dynamical configurations for $32^3\times 64$ lattices with a spatial size of 2.7 fm. Both charm and strange quark propagators are calculated with the overlap fermion action. The calculated scalar meson at 2304(22) MeV and axial-vector meson at 2546(27) MeV are in good agreement with the experimental masses of $D{s0}^*$(2317) and $D_{s1}$(2536).

preprint2008arXiv

Dust Processing and Grain Growth in Protoplanetary Disks in the Taurus-Auriga Star-Forming Region

Mid-infrared spectra of 65 T Tauri stars (TTS) taken with the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope are modeled using dust at two temperatures to probe the radial variation in dust composition in the uppermost layers of protoplanetary disks. Most spectra indicating crystalline silicates require Mg-rich minerals and silica, but a few suggest otherwise. Spectra indicating abundant enstatite at higher temperatures also require crystalline silicates at temperatures lower than those required for spectra showing high abundance of other crystalline silicates. A few spectra show 10 micron complexes of very small equivalent width. They are fit well using abundant crystalline silicates but very few large grains, inconsistent with the expectation that low peak-to-continuum ratio of the 10 micron complex always indicates grain growth. Most spectra in our sample are fit well without using the opacities of large crystalline silicate grains. If large grains grow by agglomeration of submicron grains of all dust types, the amorphous silicate components of these aggregates must typically be more abundant than the crystalline silicate components. Crystalline silicate abundances correlate positively with other such abundances, suggesting that crystalline silicates are processed directly from amorphous silicates and that neither forsterite, enstatite, nor silica are intermediate steps when producing either of the other two. Disks with more dust settling typically have greater crystalline abundances. Large-grain abundance is somewhat correlated with greater settling of disks. The lack of strong correlation is interpreted to mean that settling of large grains is sensitive to individual disk properties. Lower-mass stars have higher abundances of large grains in their inner regions.

preprint2008arXiv

PAH emission from Herbig AeBe stars

We present spectra of a sample of Herbig Ae and Be (HAeBe) stars obtained with the Infrared Spectrograph on the Spitzer Space Telescope. All but one of the Herbig stars show emission from polycyclic aromatic hydrocarbons (PAHs) and seven of the spectra show PAH emission, but no silicate emission at 10 microns. The central wavelengths of the 6.2, 7.7--8.2, and 11.3 micron emission features decrease with stellar temperature, indicating that the PAHs are less photo-processed in cooler radiation fields. The apparent low level of photo processing in HAeBe stars, relative to other PAH emission sources, implies that the PAHs are newly exposed to the UV-optical radiation fields from their host stars. HAeBe stars show a variety of PAH emission intensities and ionization fractions, but a narrow range of PAH spectral classifications based on positions of major PAH feature centers. This may indicate that, regardless of their locations relative to the stars, the PAH molecules are altered by the same physical processes in the proto-planetary disks of intermediate-mass stars. Analysis of the mid-IR spectral energy distributions indicates that our sample likely includes both radially flared and more flattened/settled disk systems, but we do not see the expected correlation of overall PAH emission with disk geometry. We suggest that the strength of PAH emission from HAeBe stars may depend not only on the degree of radial flaring, but also on the abundance of PAHs in illuminated regions of the disks and possibly on the vertical structure of the inner disk as well.

preprint2008arXiv

The relations among 8, 24, and 160 micron dust emission within nearby spiral galaxies

We investigate the relations among the stellar continuum-subtracted 8 micron polycyclic aromatic hydrocarbon (PAH 8 micron) emission, 24 micron hot dust emission, and 160 micron cold dust emission in fifteen nearby face-on spiral galaxies in the Spitzer Infrared Nearby Galaxies Survey sample. The relation between PAH 8 and 24 micron emission measured in ~2 kpc regions is found to exhibit a significant amount of scatter, and strong spatial variations are observed in the (PAH 8 micron)/24 micron surface brightness ratio. In particular, the (PAH 8 micron)/24 micron surface brightness ratio is observed to be high in the diffuse interstellar medium and low in bright star-forming regions and other locations with high 24 micron surface brightness. PAH 8 micron emission is found to be well-correlated with 160 micron emission on spatial scales of ~2 kpc, and the (PAH 8 micron)/160 micron surface brightness ratio is generally observed to increase as the 160 micron surface brightness increases. These results suggest that the PAHs are associated with the diffuse, cold dust that produces most of the 160 micron emission in these galaxies, and the variations in the (PAH 8 micron)/160 micron ratio may generally be indicative of either the intensity or the spectrum of the interstellar radiation field that is heating both the PAHs and the diffuse interstellar dust.

preprint2007arXiv

Dust and Atomic Gas in Dwarf Irregular Galaxies of the M81 Group: The SINGS and THINGS view

[abridged] We present observations of the dust and atomic gas phase in seven dwarf irregular galaxies of the M81 group from the SINGS and THINGS surveys. The Spitzer observations provide a first glimpse of the nature of the non-atomic ISM in these metal-poor (Z~0.1 Z_sun), quiescent (SFR~0.001-0.1 M_sun/yr) dwarf galaxies. Dust emission is detected in five out of the seven targets. Most detected dust emission is restricted to HI column densities >1x10^21 cm^-2. Spitzer spectroscopy of two regions in the brightest galaxies (IC 2574 and Holmberg II) show distinctly different spectral shapes. The spectrum of IC 2574 shows aromatic features that are less luminous (relative to the FIR luminosity) compared to an average SINGS spiral galaxy by a factor of \~7 . The aromatic features in Holmberg~II (which has only a slightly lower gas-phase metallicity) are fainter than in IC 2574 by an order of magnitude. This result emphazises that the strength of the aromatic features is not a simple linear function of metallicity. We estimate dust masses of ~10^4-10^6 M_sun for the M81 dwarf galaxies, resulting in an average dust--to--gas ratio (M_dust/M_HI) of ~3x10^-4 (1.5x10^-3 if only the HI that is associated with dust emission is considered); this is an order of magnitude lower than the typical value derived for the SINGS spirals. The dwarf galaxies are underluminous per unit star formation rate at 70um as compared to the more massive galaxies in SINGS by a factor of ~2. However, the average 70um/160um ratio in the sample dwarf galaxies is higher than what is found in the other galaxies of the SINGS sample. This can be explained by a combination of a lower dust content in conjunction with a higher dust temperature in the dwarfs.

preprint2007arXiv

The Calibration of Mid-Infrared Star Formation Rate Indicators

With the goal of investigating the degree to which the mid-infrared emission traces the star formation rate (SFR), we analyze Spitzer 8 um and 24 um data of star-forming regions in a sample of 33 nearby galaxies with available HST/NICMOS images in the Paschen-alpha (1.8756 um) emission line. The galaxies are drawn from the Spitzer Infrared Nearby Galaxies Survey (SINGS) sample, and cover a range of morphologies and a factor ~10 in oxygen abundance. Published data on local low-metallicity starburst galaxies and Luminous Infrared Galaxies are also included in the analysis. Both the stellar-continuum-subtracted 8 um emission and the 24 um emission correlate with the extinction-corrected Pa-alpha line emission, although neither relationship is linear. Simple models of stellar populations and dust extinction and emission are able to reproduce the observed non-linear trend of the 24 um emission versus number of ionizing photons, including the modest deficiency of 24 um emission in the low metallicity regions, which results from a combination of decreasing dust opacity and dust temperature at low luminosities. Conversely, the trend of the 8 um emission as a function of the number of ionizing photons is not well reproduced by the same models. The 8 um emission is contributed, in larger measure than the 24 um emission, by dust heated by non-ionizing stellar populations, in agreement with previous findings. Two SFR calibrations, one using the 24 um emission and the other using a combination of the 24 um and H-alpha luminosities (Kennicutt et al. 2007), are presented. No calibration is presented for the 8 um emission, because of its significant dependence on both metallicity and environment. The calibrations presented here should be directly applicable to systems dominated by on-going star formation.

preprint2007arXiv

Variations in 24 micron morphologies among galaxies in the Spitzer Infrared Nearby Galaxies Survey: New insights into the Hubble sequence

To study the distribution of star formation and dust emission within nearby galaxies, we measured five morphological parameters in the 3.6 and 24 micron wave bands for 65 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS) and 8 galaxies that were serendipitously observed by SINGS. The morphological parameters demonstrate strong variations along the Hubble sequence, including statistically significant differences between S0/a-Sab and Sc-Sd galaxies. Early-type galaxies are generally found to be compact, centralized, symmetric sources in the 24 micron band, while late-type galaxies are generally found to be extended, asymmetric sources. These results suggest that the processes that increase the real or apparent sizes of galaxies' bulges also lead to more centralized 24 micron dust emission. Several phenomena, such as strong nuclear star formation, Seyfert activity, or outer ring structures, may cause galaxies to deviate from the general morphological trends observed at 24 microns. We also note that the 24 micron morphologies of Sdm-Im galaxies are quite varied, with some objects appearing very compact and symmetric while others appear diffuse and asymmetric. These variations reflect the wide variation in star formation in irregular galaxies as observed at other wavelengths. The variations in the 24 micron morphological parameters across the Hubble sequence mirror many of the morphological trends seen in other tracers of the ISM and in stellar emission. However, the 24 micron morphological parameters for the galaxies in this sample do not match the morphological parameters measured in the stellar wave bands. This implies that the distribution of dust emission is related to but not equivalent to the distribution of stellar emission.

preprint2006arXiv

The Mid-Infrared Spectrum of Star-Forming Galaxies: Global Properties of PAH Emission

We present a sample of low-resolution 5-38um Spitzer IRS spectra of the inner few square kiloparsecs of 59 nearby galaxies spanning a large range of star formation properties. A robust method for decomposing mid-infrared galaxy spectra is described, and used to explore the behavior of PAH emission and the prevalence of silicate dust extinction. Evidence for silicate extinction is found in ~1/8 of the sample, at strengths which indicate most normal galaxies undergo A_V < ~3 magnitudes averaged over their centers. The contribution of PAH emission to the total infrared power is found to peak near 10% and extend up to ~20%, and is suppressed at metallicities Z < ~Z_sun/4, as well as in low-luminosity AGN environments. Strong inter-band PAH feature strength variations (2-5x) are observed, with the presence of a weak AGN and, to a lesser degree, increasing metallicity shifting power to the longer wavelength bands. A peculiar PAH emission spectrum with markedly diminished 5-8um features arises among the sample solely in systems with relatively hard radiation fields harboring low-luminosity AGN. The AGN may modify the emitting grain distribution and provide the direct excitation source of the unusual PAH emission, which cautions against using absolute PAH strength to estimate star formation rates in systems harboring active nuclei. Alternatively, the low star formation intensity often associated with weak AGN may affect the spectrum. The effect of variations in the mid-infrared spectrum on broadband infrared surveys is modeled, and points to more than a factor of two uncertainty in results which assume a fixed PAH emission spectrum, for redshifts z=0-2.5.

preprint2004arXiv

Hot Nuclear Matter Equation of State with a Three-body Force

The finite temperature Brueckner-Hartree-Fock approach is extended by introducing a microscopic three-body force. In the framework of the extended model, the equation of state of hot asymmetric nuclear matter and its isospin dependence have been investigated. The critical temperature of liquid-gas phase transition for symmetric nuclear matter has been calculated and compared with other predictions. It turns out that the three-body force gives a repulsive contribution to the equation of state which is stronger at higher density and as a consequence reduces the critical temperature of liquid-gas phase transition. The calculated energy per nucleon of hot asymmetric nuclear matter is shown to satisfy a simple quadratic dependence on asymmetric parameter $β$ as in the zero-temperature case. The symmetry energy and its density dependence have been obtained and discussed. Our results show that the three-body force affects strongly the high-density behavior of the symmetry energy and makes the symmetry energy more sensitive to the variation of temperature. The temperature dependence and the isospin dependence of other physical quantities, such as the proton and neutron single particle potentials and effective masses are also studied. Due to the additional repulsion produced by the three-body force contribution, the proton and neutron single particle potentials are correspondingly enhanced as similar to the zero-temperature case.

preprint2004arXiv

Mid-Infrared IRS Spectroscopy of NGC 7331: A First Look at the SINGS Legacy

The nearby spiral galaxy NGC 7331 was spectrally mapped from 5-38um using all modules of Spitzer's IRS spectrograph. A strong new dust emission feature, presumed due to PAHs, was discovered at 17.1um. The feature's intensity is nearly half that of the ubiquitous 11.3um band. The 7-14um spectral maps revealed significant variation in the 7.7 and 11.3um PAH features between the stellar ring and nucleus. Weak [OIV] 25.9um line emission was found to be centrally concentrated in the nucleus, with an observed strength over 10% of the combined neon line flux, indicating an AGN or unusually active massive star photo-ionization. Two [SIII] lines fix the characteristic electron density in the HII regions at n_e < ~200 cm^-3. Three detected H_2 rotational lines, tracing warm molecular gas, together with the observed IR continuum, are difficult to match with standard PDR models. Either additional PDR heating or shocks are required to simultaneously match lines and continuum.

preprint2003arXiv

SINGS: The SIRTF Nearby Galaxies Survey

The SIRTF Nearby Galaxy Survey is a comprehensive infrared imaging and spectroscopic survey of 75 nearby galaxies. Its primary goal is to characterize the infrared emission of galaxies and their principal infrared-emitting components, across a broad range of galaxy properties and star formation environments. SINGS will provide new insights into the physical processes connecting star formation to the interstellar medium properties of galaxies, and provide a vital foundation for understanding infrared observations of the distant universe and ultraluminous and active galaxies. The galaxy sample and observing strategy have been designed to maximize the scientific and archival value of the data set for the SIRTF user community at large. The SIRTF images and spectra will be supplemented by a comprehensive multi-wavelength library of ancillary and complementary observations, including radio continuum, HI, CO, submillimeter, BVRIJHK, H-alpha, Paschen-alpha, ultraviolet, and X-ray data. This paper describes the main astrophysical issues to be addressed by SINGS, the galaxy sample and the observing strategy, and the SIRTF and other ancillary data products.