Search

Purpose: To present a novel generalized MR image reconstruction based on pseudoinversion of the encoding matrix (Pinv-Recon) as a simple yet powerful method, and demonstrate its computational feasibility for diverse MR imaging applications. Methods: MR image encoding constitutes a linear mapping of the unknown image to the measured k-space data mediated via an encoding matrix ($ data = Encode \times image$). Pinv-Recon addresses MR image reconstruction as a linear inverse problem ($image = Encode^{-1} \times data$), explicitly calculating the Moore-Penrose pseudoinverse of the encoding matrix using truncated singular value decomposition (tSVD). Using a discretized, algebraic notation, we demonstrate constructing a generalized encoding matrix by stacking relevant encoding mechanisms (e.g., gradient encoding, coil sensitivity encoding, chemical shift inversion) and encoding distortions (e.g., off-center positioning, B$_0$ inhomogeneity, spatiotemporal gradient imperfections, transient relaxation effects). Iterative reconstructions using the explicit generalized encoding matrix, and the computation of the spatial-response-function (SRF) and noise amplification, were demonstrated. Results: We evaluated the computation times and memory requirements (time ~ (size of the encoding matrix)$^{1.4}$). Using the Shepp-Logan phantom, we demonstrated the versatility of the method for various intertwined MR image encoding and distortion mechanisms, achieving better MSE, PSNR and SSIM metrics than conventional methods. A diversity of datasets, including the ISMRM CG-SENSE challenge, were used to validate Pinv-Recon. Conclusion: Although pseudo-inversion of large encoding matrices was once deemed computationally intractable, recent advances make Pinv-Recon feasible. It has great promise for both research and clinical applications, and for educational use., Comment: 26 pages, 8 figures, 2 tables (+ Supplementary Material). Submitted to Magnetic Resonance in Medicine

Context. Magnetic bright points (MBPs) are one of the smallest manifestations of the magnetic field in the solar atmosphere and are observed to extend from the photosphere up to the chromosphere. As such, they represent an excellent feature to use in searches for types of magnetohydrodynamic (MHD) waves and mode coupling in the solar atmosphere. Aims. In this work, we aim to study wave propagation in the lower solar atmosphere by comparing intensity oscillations in the photosphere with the chromosphere via a search for possible mode coupling, in order to establish the importance of these types of waves in the solar atmosphere, and their contribution to heating the chromosphere. Methods. Observations were conducted in July 2011 with the ROSA and the HARDCam instruments at the Dunn Solar Telescope. We used wavelet analysis to identify traveling MHD waves and derive frequencies in the G-band and H$\alpha$wave bands. We isolated a large sample of MBPs using an automated tracking algorithm throughout our observations. Two dozen of the brightest MBPs were selected from the sample for further study. Results. We find oscillations in the G-band MBPs, with frequencies between 1.5 and 3.6 mHz. Corresponding MBPs in the lower solar chromosphere observed in H$\alpha$ show a frequency range of 1.4 to 4.3 mHz. In about 38\% of the MBPs, the ratio of H$\alpha$ to G-band frequencies was near two. Thus, these oscillations show a form of mode coupling where the transverse waves in the photosphere are converted into longitudinal waves in the chromosphere. Conclusions. From simple estimates we find an energy flux of $\approx$45 $\times 10^{3}$ W m$^{-2}$ and show that the energy flowing through MBPs is enough to heat the chromosphere, and mode coupling is important in helping us understand the types of MHD waves in the lower solar atmosphere and the overall energy budget., Comment: 6 pages, 4 figures, 1 Table; Accepted in Astronomy \& Astrophysics

We study the spin-$S$ Kitaev-Heisenberg model on the honeycomb lattice for $S\!=\!1/2$, $1$ and $3/2$, by using the coupled cluster method (CCM) of microscopic quantum many-body theory. This system is one of the earliest extensions of the Kitaev model and is believed to contain two extended spin liquid phases for any value of the spin quantum number $S$. We show that the CCM delivers accurate estimates for the phase boundaries of these spin liquid phases, as well as other transition points in the phase diagram. Moreover, we find evidence of two unexpected narrow phases for $S\!=\!1/2$, one sandwiched between the zigzag and ferromagnetic phases and the other between the N\'eel and the stripy phases. The results establish the CCM as a versatile numerical technique that can capture the strong quantum-mechanical fluctuations that are inherently present in generalized Kitaev models with competing bond-dependent anisotropies., Comment: 16 pages, 10 figures. This version: figures 5,8 and 10 updated, Table I updated, references added

With the construction of the European Spallation Source, a remarkable opportunity has emerged to conduct high sensitivity searches for neutron oscillations, including a first search for thirty years for free neutrons converting to antineutrons. Furthermore, searches can be made for transitions of neutrons and antineutrons to sterile neutron states. Upgrades to the ESS infrastructure allow an improved HIBEAM design that would provide an increase in sensitivity by an order of magnitude compared to previous work. The HIBEAM program focuses on processes that violate baryon number by one or two units. The observation of a process satisfying one of the Sakharov conditions addresses the open question of the origin of the matter-antimatter asymmetry in the Universe. Sterile neutron states would belong to a `dark' sector of particles which may explain dark matter. As electrically neutral, meta-stable objects that can be copiously produced and studied, neutrons represent an attractive portal to a `dark' sector. The HIBEAM instrument can also be utilised for other purposes such as direct searches for ultralight axion dark matter. This paper describes the capability, design, infrastructure, and potential of the HIBEAM program. This includes a dedicated beamline, neutron optical system, magnetic shielding and control, and detectors for neutrons and antineutrons., Comment: 58 pages, 23 figures

The European Spallation Source, currently under construction in Lund, Sweden, is a multidisciplinary international laboratory. Once completed to full specifications, it will operate the world's most powerful pulsed neutron source. Supported by a 3 million Euro Research and Innovation Action within the EU Horizon 2020 program, a design study (HighNESS) has been completed to develop a second neutron source located below the spallation target. Compared to the first source, designed for high cold and thermal brightness, the new source has been optimized to deliver higher intensity, and a shift to longer wavelengths in the spectral regions of cold (CN, 2--20\,\AA), very cold (VCN, 10--120\,\AA), and ultracold (UCN, ${>}\,{500}$\,\AA) neutrons. The second source comprises a large liquid deuterium moderator designed to produce CN and support secondary VCN and UCN sources. Various options have been explored in the proposed designs, aiming for world-leading performance in neutronics. These designs will enable the development of several new instrument concepts and facilitate the implementation of a high-sensitivity neutron-antineutron oscillation experiment (NNBAR). This document serves as the Conceptual Design Report for the HighNESS project, representing its final deliverable., Comment: 269 pages, 255 figures. Volume I of the final deliverable of the HighNESS Project (HORIZON 2020 grant agreement ID: 951782)

In this paper, we study alternating links in thickened surfaces in terms of the lattices of integer flows on their Tait graphs. We use this approach to give a short proof of the first two generalised Tait conjectures. We also prove that the flow lattice is an invariant of alternating links in thickened surfaces and is further invariant under disc mutation. For classical links, the flow lattice and $d$-invariants are complete invariants of the mutation class of an alternating link. For links in thickened surfaces, we show that this is no longer the case by finding a stronger mutation invariant, namely the Gordon-Litherland linking form. In particular, we find alternating knots in thickened surfaces which have isometric flow lattices but with non-isomorphic linking forms., Comment: 34 pages, 11 figures, and 4 tables

Feedback from active galactic nuclei (AGN) has long been invoked to explain the correlation between black hole mass and stellar velocity dispersion (M-$\sigma$) discovered in low redshift galaxies. We describe the time evolution of AGN in the M-$\sigma$ plane based on our gap model (Garofalo, Evans \& Sambruna 2010) for black hole accretion and jet formation illustrating a fundamental difference between jetted and non-jetted AGN. While the latter tend to evolve diagonally upward with black hole mass increasing along with stellar dispersion, we show that jetted AGN tend on average to move initially more upwards because their effect on velocity dispersion is weaker than for non-jetted AGN. But this initial phase is followed by a shift in the nature of the feedback, from positive to negative, a transition that is more dramatic on average in denser cluster environments. The feedback gets its kick from tilted jets which shut down star formation but increase velocity dispersion values. As this change in the nature of the feedback takes tens of millions to hundreds of millions of years, these cluster, merger-triggered jetted AGN, will evolve more upwards for up to order $10^{8}$ years, followed by an extremely long phase in which low excitation progressively slows black hole growth but dramatically affects stellar dispersion. As a result, powerful jetted AGN evolve for most of their lives almost horizontally on the M-$\sigma$ plane. The prediction is that strongest AGN feedback on stellar dispersion is a late universe phenomenon with M87 a good example. We show how jetted and non-jetted AGN parallel the Sersic and core-Sersic galaxy paths in the M-$\sigma$ plane found by Sahu et al (2019) and to a prediction that jetted quasars are not core-Sersic galaxies as found for lower redshift jetted AGN., Comment: Open Journal of Astrophysics

Presently under construction in Lund, Sweden, the European Spallation Source (ESS) will be the world's brightest neutron source. As such, it has the potential for a particle physics program with a unique reach and which is complementary to that available at other facilities. This paper describes proposed particle physics activities for the ESS. These encompass the exploitation of both the neutrons and neutrinos produced at the ESS for high precision (sensitivity) measurements (searches)., Comment: 121 pages, updated version after referee comments

The NNBAR experiment for the European Spallation Source will search for free neutrons converting to antineutrons with a sensitivity improvement of three orders of magnitude compared to the last such search. This paper describes progress towards a conceptual design report for NNBAR. The design of a moderator, neutron reflector, beamline, shielding and annihilation detector is reported. The simulations used form part of a model which will be used for optimisation of the experiment design and quantification of its sensitivity., Comment: 30 pages, 26 figures, accepted for publication in Journal of Instrumentation (JINST)

We present spectroscopic neutron imaging (SNI), a bridge between imaging and scattering techniques, for the analysis of hydrogenated molecules in lithium-ion cells. The scattering information of CHn-based organic solvents and electrolytes was mapped in two-dimensional space by investigating the wavelength-dependent property of hydrogen atoms through time-of-flight imaging. Our investigation demonstrates a novel approach to detect physical and chemical changes in hydrogenated liquids, which extends, but not limits, the use of SNI to relevant applications in electrochemical devices, e.g., the study of electrolytes in Li-ion batteries.

The European Spallation Source (ESS), presently under construction in Lund, Sweden, is a multidisciplinary international laboratory that will operate the world's most powerful pulsed neutron source. Supported by a 3M Euro Research and Innovation Action within the EU Horizon 2020 program, a design study (HighNESS) is now underway to develop a second neutron source below the spallation target. Compared to the first source, located above the spallation target and designed for high cold and thermal brightness, the new source will provide higher intensity, and a shift to longer wavelengths in the spectral regions of cold (2 /- 20 {\AA}), very cold (VCN, 10 /- 120 {\AA}), and ultra cold (UCN, > 500 {\AA}) neutrons. The core of the second source will consist of a large liquid deuterium moderator to deliver a high flux of cold neutrons and to serve secondary VCN and UCN sources, for which different options are under study. The features of these new sources will boost several areas of condensed matter research and will provide unique opportunities in fundamental physics. Part of the HighNESS project is also dedicated to the development of future instruments that will make use of the new source and will complement the initial suite of instruments in construction at ESS. The HighNESS project started in October 2020. In this paper, the ongoing developments and the results obtained in the first year are described., Comment: 10 pages, 10 figures, 14th International Topical Meeting on Nuclear Applications of Accelerators, November 30 to December 4, 2021, Washington, DC

Monte Carlo radiation transport simulations have clearly contributed to improve the design of nuclear systems. When performing in-beam or shielding simulations a complexity arises due to the fact that particles must be tracked to regions far from the original source or behind the shielding, often lacking sufficient statistics. Different possibilities to overcome this problem such as using particle lists or generating synthetic sources have already been reported. In this work we present a new approach by using the adaptive multivariate kernel density estimator (KDE) method. This concept was implemented in KDSource, a general tool for modelling, optimizing and sampling KDE sources, which provides a convenient user interface. The basic properties of the method were studied in an analytical problem with a known density distribution. Furthermore, the tool was used in two Monte Carlo simulations that modelled neutron beams, which showed good agreement with experimental results.

If popular online platforms systematically expose their users to partisan and unreliable news, they could potentially contribute to societal issues like rising political polarization. This concern is central to the echo chamber and filter bubble debates, which critique the roles that user choice and algorithmic curation play in guiding users to different online information sources. These roles can be measured in terms of exposure, the URLs seen while using an online platform, and engagement, the URLs selected while on that platform or browsing the web more generally. However, due to the challenges of obtaining ecologically valid exposure data--what real users saw during their regular platform use--studies in this vein often only examine engagement data, or estimate exposure via simulated behavior or inference. Despite their centrality to the contemporary information ecosystem, few such studies have focused on web search, and even fewer have examined both exposure and engagement on any platform. To address these gaps, we conducted a two-wave study pairing surveys with ecologically valid measures of exposure and engagement on Google Search during the 2018 and 2020 US elections. We found that participants' partisan identification had a small and inconsistent relationship with the amount of partisan and unreliable news they were exposed to on Google Search, a more consistent relationship with the search results they chose to follow, and the most consistent relationship with their overall engagement. That is, compared to the news sources our participants were exposed to on Google Search, we found more identity-congruent and unreliable news sources in their engagement choices, both within Google Search and overall. These results suggest that exposure and engagement with partisan or unreliable news on Google Search are not primarily driven by algorithmic curation, but by users' own choices., Comment: Updated to latest version of manuscript. Abstract preview is trimmed to fit arXiv 1,920 character limit

The HIBEAM/NNBAR program is a proposed two-stage experiment at the European Spallation Source focusing on searches for baryon number violation via processes in which neutrons convert to antineutrons. This paper outlines the computing and detector simulation framework for the HIBEAM/NNBAR program. The simulation is based on predictions of neutron flux and neutronics together with signal and background generation. A range of diverse simulation packages are incorporated, including Monte Carlo transport codes, neutron ray-tracing simulation packages, and detector simulation software. The common simulation package in which these elements are interfaced together is discussed. Data management plans and triggers are also described., Comment: Contribution to CHEP2021. Accepted for publication in the European Physical Journal (EPJ) Web of Conferences

The experimental thermal neutron cross sections of the twenty proteinogenic amino acids have been measured over the incident-neutron energy range spanning from 1 meV to 10 keV and data have been interpreted using the multi-phonon expansion based on first-principles calculations. The scattering cross section, dominated by the incoherent inelastic contribution from the hydrogen atoms, can be rationalised in terms of the average contributions of different functional groups, thus neglecting their correlation. These results can be used for modelling the total neutron cross sections of complex organic systems like proteins, muscles, or human tissues from a limited number of starting input functions. This simplification is of crucial importance for fine-tuning of transport simulations used in medical applications, including boron neutron capture therapy as well as secondary neutrons-emission induced during proton therapy. Moreover, the parametrized neutron cross sections allow a better treatment of neutron scattering experiments, providing detailed sample self-attenuation corrections for a variety of biological and soft-matter systems.

Purpose: Phosphorus saturation-transfer experiments can quantify metabolic fluxes non-invasively. Typically, the forward flux through the creatine-kinase reaction is investigated by observing the decrease in phosphocreatine (PCr) after saturation of $\gamma$-ATP. The quantification of total ATP utilisation is currently under-explored, as it requires simultaneous saturation of inorganic phosphate (Pi) and PCr. This is challenging, as currently available saturation pulses reduce the already-low $\gamma$-ATP signal present. Methods: Using a hybrid optimal-control and Shinnar-Le-Roux method, a quasi-adiabatic RF pulse was designed for the dual-saturation of PCr and Pi to enable determination of total ATP utilisation. The pulses were evaluated in Bloch equation simulations, compared with a conventional hard-cosine DANTE saturation sequence, before application to perfused rat hearts at 11.7 Tesla. Results: The quasi-adiabatic pulse was insensitive to a $>2.5$-fold variation in $B_1$, producing equivalent saturation with a 53% reduction in delivered pulse power and a 33-fold reduction in spillover at the minimum effective $B_1$. This enabled the complete quantification of the synthesis and degradation fluxes for ATP in 30-45 minutes in the perfused rat heart. While the net synthesis flux ($4.24\pm0.8$ mM/s, SEM) was not significantly different from degradation flux ($6.88\pm2$ mM/s, $p=0.06$) and both measures are consistent with prior work, nonlinear error analysis highlights uncertainties in the Pi-to-ATP measurement that may explain a trend suggesting a possible imbalance. Conclusion: This work demonstrates a novel quasi-adiabatic dual-saturation RF pulse with significantly improved performance that can be used to measure ATP turnover in the heart in vivo., Comment: 26 pages, Accepted at Magnetic Resonance in Medicine, 24/11/2020 [This version post reviews]

Purpose: Hyperpolarized imaging experiments have conflicting requirements of high spatial, temporal, and spectral resolution. Spectral-Spatial RF excitation has been shown to form an attractive magnetization-efficient method for hyperpolarized imaging, but the optimum readout strategy is not yet known. Methods: In this work we propose a novel 3D hybrid-shot spiral sequence which features two constant density regions that permit the retrospective reconstruction of either high spatial or high temporal resolution images post hoc, (adaptive spatiotemporal imaging) allowing greater flexibility in acquisition and reconstruction. Results: We have implemented this sequence, both via simulation and on a pre-clinical scanner, to demonstrate its feasibility, in both a 1H phantom and with hyperpolarized 13C pyruvate in vivo. Conclusion: This sequence forms an attractive method for acquiring hyperpolarized imaging datasets, providing adaptive spatiotemporal imaging to ameliorate the conflict of spatial and temporal resolution, with significant potential for clinical translation.

By exploring more than sixty thousand quasars from the Sloan Digital Sky Survey Data Release 5, Steinhardt & Elvis discovered a sub-Eddington boundary and a redshift-dependent drop-off at higher black hole mass, possible clues to the growth history of massive black holes. Our contribution to this special issue of Universe amounts to an application of a model for black hole accretion and jet formation to these observations. For illustrative purposes we include about 100,000 data points from the Sloan Digital Sky Survey Data Release 7 where the sub-Eddington boundary is also visible, and propose a theoretical picture that explains these features. By appealing to thin disk theory and both the lower accretion efficiency and the time evolution of jetted quasars compared to non-jetted quasars in our 'gap paradigm', we explain two features of the sub-Eddington boundary. First, we show that a drop-off on the quasar mass-luminosity plane for larger black hole mass occurs at all redshifts. But the fraction of jetted quasars is directly related to the merger function in this paradigm, which means the jetted quasar fraction drops with decrease in redshift, which allows us to explain a second feature of the sub-Eddington boundary, namely a redshift dependence of the slope of the quasar mass-luminosity boundary at high black hole mass stemming from a change in radiative efficiency with time. We are able to reproduce the mass dependence of, as well as the oscillating behavior in, the slope of the sub-Eddington boundary as a function of time. The basic physical idea involves retrograde accretion occurring only for a subset of the more massive black holes which implies that most spinning black holes in our model are prograde accretors. In short, this paper amounts to a qualitative overview of how a sub-Eddington boundary naturally emerges in the gap paradigm., Comment: Special Issue of Universe: Recent Progress in Relativistic Astrophysics

Over the last two decades, the discovery of exoplanets has fundamentally changed our perception of the universe and humanity's place within it. Recent work indicates that a solar system's X-ray and high energy particle environment is of fundamental importance to the formation and development of the atmospheres of close-in planets such as hot Jupiters, and Earth-like planets around M stars. X-ray imaging and spectroscopy provide powerful and unique windows into the high energy flux that an exoplanet experiences, and X-ray photons also serve as proxies for potentially transfigurative coronal mass ejections. Finally, if the host star is a bright enough X-ray source, transit measurements akin to those in the optical and infrared are possible and allow for direct characterization of the upper atmospheres of exoplanets. In this brief white paper, we discuss contributions to the study of exoplanets and their environs which can be made by X-ray data of increasingly high quality that are achievable in the next 10--15 years., Comment: 7 pages with 2 figures. A 2020 Decadal Survey White Paper

It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current instrumentation. In line with theCommittee on an Exoplanet Science Strategy recognition that holistic observational approaches are needed, we point out here that a full understanding of exoplanet atmospheres, their evolution and determination of habitability requires a powerful high-resolution X-ray imaging and spectroscopic observatory. This is the only capability that can: (1) characterize by proxy the crucial, difficult to observe, EUV stellar flux, its history and its variations for planet hosting stars; (2) observe the stellar wind; (3) detect the subtle Doppler signatures of coronal mass ejections., Comment: Astro2020 Decadal Survey Science White Paper

Small-scale magnetic reconnection processes, in the form of nanoflares, have become increasingly hypothesized as important mechanisms for the heating of the solar atmosphere, for driving propagating disturbances along magnetic field lines in the Sun's corona, and for instigating rapid jet-like bursts in the chromosphere. Unfortunately, the relatively weak signatures associated with nanoflares places them below the sensitivities of current observational instrumentation. Here, we employ Monte Carlo techniques to synthesize realistic nanoflare intensity time series from a dense grid of power-law indices and decay timescales. Employing statistical techniques, which examine the modeled intensity fluctuations with more than 10^7 discrete measurements, we show how it is possible to extract and quantify nanoflare characteristics throughout the solar atmosphere, even in the presence of significant photon noise. A comparison between the statistical parameters (derived through examination of the associated intensity fluctuation histograms) extracted from the Monte Carlo simulations and SDO/AIA 171{\AA} and 94{\AA} observations of active region NOAA 11366 reveals evidence for a flaring power-law index within the range of 1.82 - 1.90, combined with e-folding timescales of 385 +/- 26 s and 262 +/- 17 s for the SDO/AIA 171{\AA} and 94{\AA} channels, respectively. These results suggest that nanoflare activity is not the dominant heating source for the active region under investigation. This opens the door for future dedicated observational campaigns to not only unequivocally search for the presence of small-scale reconnection in solar and stellar environments, but also quantify key characteristics related to such nanoflare activity., Comment: 26 pages, 6 figures, published in ApJ

We present multi-wavelength observations of an M-class flare (M3.9) that occurred on 2014 June 11. Our observations were conducted with the Dunn Solar Telescope (DST), adaptive optics, the multi-camera system ROSA (Rapid Oscillations in Solar Atmosphere) and new HARDcam (Hydrogen-Alpha Rapid Dynamics) camera in various wavelengths, such as Ca~II~K, Mg~I~b$_2$ (at 5172.7 Ang), and H$\alpha$ narrow-band, and G-band continuum filters. Images were re-constructed using the Kiepencheuer-Institut Speckle Interferometry Package (KISIP) code, to improve our image resolution. We observed intensity increases of $\approx$120-150% in the Mg, Ca~K and H$\alpha$ narrow band filters during the flare. Intensity increases for the flare observed in the SDO EUV channels were several times larger, and the GOES X-rays increased over a factor of 30 for the harder band. Only a modest delay is found between the onset of flare ribbons of a nearby sympathetic flare and the main flare ribbons observed in these narrow-band filters. The peak flare emission occurs within a few seconds for the Ca~K, Mg, and H$\alpha$ bands. Time-distance techniques find propagation velocities of $\approx$60 km/s for the main flare ribbon and as high as 300 km/s for smaller regions we attribute to filament eruptions. This result and delays and velocities observed with SDO ($\approx$100 km/s) for different coronal heights agree well with the simple model of energy propagation versus height, although a more detailed model for the flaring solar atmosphere is needed. And finally, we detected marginal quasi-periodic pulsations (QPPs) in the 40--60 second range for the Ca~K, Mg and H$\alpha$ bands, and such measurements are important for disentangling the detailed flare-physics., Comment: 16 Pages, 7 Figures, 1 Table (1 video in on-line journal); Accepted in Research in Astronomy and Astrophysics