Search

We explore the impact of retaining three-body operators within the in-medium similarity renormalization group (IMSRG), as well as various approximations schemes. After studying two toy problems, idential fermions with a contact interaction and the Lipkin-Meshkov-Glick model, we employ the valence-space formulation of the IMSRG to investigate the even-$A$ carbon isotopes with a chiral two-body potential. We find that retaining only those commutators expressions that scale as $N^7$ provides an excellent approximation of the full three-body treatment., Comment: 18 pages, 11 figures

Autonomous methods to align beamlines can decrease the amount of time spent on diagnostics, and also uncover better global optima leading to better beam quality. The alignment of these beamlines is a high-dimensional, expensive-to-sample optimization problem involving the simultaneous treatment of many optical elements with correlated and nonlinear dynamics. Bayesian optimization is a strategy of efficient global optimization that has proved successful in similar regimes in a wide variety of beamline alignment applications, though it has typically been implemented for particular beamlines and optimization tasks. In this paper, we present a basic formulation of Bayesian inference and Gaussian process models as they relate to multiobjective Bayesian optimization, as well as the practical challenges presented by beamline alignment. We show that the same general implementation of Bayesian optimization with special consideration for beamline alignment can quickly learn the dynamics of particular beamlines in an online fashion through hyperparameter fitting with no prior information. We present the implementation of a concise software framework for beamline alignment and test it on four different optimization problems for experiments at x-ray beamlines of the National Synchrotron Light Source II and the Advanced Light Source and an electron beam at the Accelerator Test Facility, along with benchmarking on a simulated digital twin. We discuss new applications of the framework, and the potential for a unified approach to beamline alignment at synchrotron facilities.

Submillimeter single-dish telescopes offer two key advantages compared to interferometers: they can efficiently map larger portions of the sky and recover larger spatial scales. Nonetheless, fluctuations in the atmosphere limit the accurate retrieval of signals from astronomical sources. Therefore, we introduce a user-friendly simulator named ${\tt maria}$ to optimize scanning strategies and instrument designs to efficiently reduce atmospheric noise and filtering effects. We further use this tool to produce synthetic time streams and maps from hydrodynamical simulations, enabling a fair comparison between theory and reality. ${\tt maria}$ has implemented a suite of telescope and instrument designs intended to mimic current and future facilities. To generate synthetic time-ordered data, each mock observatory scans through the atmosphere in a configurable pattern over the celestial object. We generate evolving and location-and-time-specific weather for each of the fiducial sites using a combination of satellite and ground-based measurements. While ${\tt maria}$ is a generic virtual telescope, this study specifically focuses on mimicking broadband bolometers observing at 100 GHz. To validate our virtual telescope, we compare the mock time streams with real MUSTANG-2 observations and find that they are quantitatively similar by conducting a k-sample Anderson-Darling test resulting in p<0.001. Subsequently, we image the time-ordered data to create noise maps and mock observations of clusters of galaxies for both MUSTANG-2 and an instrument concept for the 50m Atacama Large Aperture Submillimeter Telescope (AtLAST). Furthermore, using ${\tt maria}$, we find that a 50m dish provides the highest levels of correlation of atmospheric signals across adjacent detectors compared to smaller apertures (e.g., 42-cm and 6-m survey experiments), facilitating removal of atmospheric signal on large scales.

HiPERCAM is a portable, quintuple-beam optical imager that saw first light on the 10.4-m Gran Telescopio Canarias (GTC) in 2018. The instrument uses re-imaging optics and 4 dichroic beamsplitters to record $u_s g_s r_s i_s z_s$ ($320-1060$ nm) images simultaneously on its five CCD cameras, each of 3.1 arcmin (diagonal) field of view. The detectors in HiPERCAM are frame-transfer devices cooled thermo-electrically to 183 K, thereby allowing both long-exposure, deep imaging of faint targets, as well as high-speed (over 1000 windowed frames per second) imaging of rapidly varying targets. A comparison-star pick-off system in the telescope focal plane increases the effective field of view to 6.7 arcmin for differential photometry. Combining HiPERCAM with the world's largest optical telescope enables the detection of astronomical sources to $g_s \sim 23$ in 1 s and $g_s \sim 28$ in 1 h. In this paper we describe the scientific motivation behind HiPERCAM, present its design, report on its measured performance, and outline some planned enhancements., Comment: 18 pages, 12 figures, 1 table, accepted for publication in Monthly Notices of the Royal Astronomical Society

The powerful combination of the cutting-edge multi-object spectrograph MOSAIC with the world largest telescope, the ELT, will allow us to probe deeper into the Universe than was possible. MOSAIC is an extremely efficient instrument in providing spectra for the numerous faint sources in the Universe, including the very first galaxies and sources of cosmic reionization. MOSAIC has a high multiplex in the NIR and in the VIS, in addition to multi-Integral Field Units (Multi-IFUs) in NIR. As such it is perfectly suited to carry out an inventory of dark matter (from rotation curves) and baryons in the cool-warm gas phases in galactic haloes at z=3-4. MOSAIC will enable detailed maps of the intergalactic medium at z=3, the evolutionary history of dwarf galaxies during a Hubble time, the chemistry directly measured from stars up to several Mpc. Finally, it will measure all faint features seen in cluster gravitational lenses or in streams surrounding nearby galactic halos, providing MOSAIC to be a powerful instrument with an extremely large space of discoveries. The preliminary design of MOSAIC is expected to begin next year, and its level of readiness is already high, given the instrumental studies made by the team., Comment: 8 pages, 4 Figures, Pre-edited version, to appear in the ESO Messenger No.182 - Quarter 1 2021 - Version identical to the Edited one

We use coupled-cluster theory and nuclear interactions from chiral effective field theory to compute the nuclear matrix element for the neutrinoless double-beta decay of $^{48}$Ca. Benchmarks with the no-core shell model in several light nuclei inform us about the accuracy of our approach. For $^{48}$Ca we find a relatively small matrix element. We also compute the nuclear matrix element for the two-neutrino double-beta decay of $^{48}$Ca with a quenching factor deduced from two-body currents in recent ab-initio calculation of the Ikeda sum-rule in $^{48}$Ca [Gysbers et al., Nature Physics 15, 428-431 (2019)]., Comment: 14 pages, 13 figures; Version accepted for publication, Supplemental material also updated

The performance of tomographic adaptive optics systems is intrinsically linked to the vertical profile of optical turbulence. Firstly, a sufficient number of discrete turbulent layers must be reconstructed to model the true continuous turbulence profile. Secondly over the course of an observation, the profile as seen by the telescope changes and the tomographic reconstructor must be updated. These changes can be due to the unpredictable evolution of turbulent layers on meteorological timescales as short as minutes. Here we investigate the effect of changing atmospheric conditions on the quality of tomographic reconstruction by coupling fast analytical adaptive optics simulation to a large database of 10 691 high resolution turbulence profiles measured over two years by the Stereo-SCIDAR instrument at ESO Paranal, Chile. This work represents the first investigation of these effects with a large, statistically significant sample of turbulence profiles. The statistical nature of the study allows us to assess not only the degradation and variability in tomographic error with a set of system parameters (e.g. number of layers, temporal update period) but also the required parameters to meet some error threshold. In the most challenging conditions where the profile is rapidly changing, these parameters must be far more tightly constrained in order to meet this threshold. By providing estimates of these constraints for a wide range of system geometries as well as the impact of different temporal optimisation strategies we may assist the designers of tomographic AO for the ELT to dimension their systems., Comment: 13 pages, 13 figures. Accepted MNRAS March 2020

We present error mitigation (EM) techniques for noisy intermediate-scale quantum computers (QC) based on density matrix purification and perturbative corrections to the target energy. We incorporate this scheme into the variational quantum eigensolver (VQE) and demonstrate chemically-accurate ground state energy calculations of various alkali metal hydrides using IBM quantum computers. Both the density matrix purification improvements and the perturbative corrections require only meager classical computational resources, and are conducted exclusively as post-processing of the measured density matrix. The improved density matrix leads to better simulation accuracy at each step of the variational optimization, resulting in a better input into the next optimization step without additional measurements. Adding perturbative corrections to the resulting energies further increases the accuracy, and decreases variation between consecutive measurements. These EM schemes allow for previously unavailable levels of accuracy over remote QC resources.

For extremely large telescopes, adaptive optics will be required to correct the Earth's turbulent atmosphere. The performance of tomographic adaptive optics is strongly dependent on the vertical distribution (profile) of this turbulence. An important way in which this manifests is the tomographic error, arising from imperfect measurement and reconstruction of the turbulent phase at altitude. Conventionally, a small number of reference profiles are used to obtain this error in simulation however these profiles are not constructed to be representative in terms of tomographic error. It is therefore unknown whether these simulations are providing realistic performance estimates. Here, we employ analytical adaptive optics simulation that drastically reduces computation times to compute tomographic error for 10 691 measurements of the turbulence profile gathered by the Stereo-SCIDAR instrument at ESO Paranal. We assess for the first time the impact of the profile on tomographic error in a statistical manner. We find, in agreement with previous work, that the tomographic error is most directly linked with the distribution of turbulence into discrete, stratified layers. Reference profiles are found to provide mostly higher tomographic error than expected, which we attribute to the fact that these profiles are primarily composed of averages of many measurements resulting in unrealistic, continuous distributions of turbulence. We propose that a representative profile should be defined with respect to a particular system, and that as such simulations with a large statistical sample of profiles must be an important step in the design process., Comment: Accepted MNRAS June 2019

We present a quantum chemistry benchmark for noisy intermediate-scale quantum computers that leverages the variational quantum eigensolver, active space reduction, a reduced unitary coupled cluster ansatz, and reduced density purification as error mitigation. We demonstrate this benchmark on the 20 qubit IBM Tokyo and 16 qubit Rigetti Aspen processors via the simulation of alkali metal hydrides (NaH, KH, RbH),with accuracy of the computed ground state energy serving as the primary benchmark metric. We further parameterize this benchmark suite on the trial circuit type, the level of symmetry reduction, and error mitigation strategies. Our results demonstrate the characteristically high noise level present in near-term superconducting hardware, but provide a relevant baseline for future improvement of the underlying hardware, and a means for comparison across near-term hardware types. We also demonstrate how to reduce the noise in post processing with specific error mitigation techniques. Particularly, the adaptation of McWeeny purification of noisy density matrices dramatically improves accuracy of quantum computations, which, along with adjustable active space, significantly extends the range of accessible molecular systems. We demonstrate that for specific benchmark settings, the accuracy metric can reach chemical accuracy when computing over the cloud on certain quantum computers.

$\beta$-decay, a process that changes a neutron into a proton (and vice versa), is the dominant decay mode of atomic nuclei. This decay offers a unique window to physics beyond the standard model, and is at the heart of microphysical processes in stellar explosions and the synthesis of the elements in the Universe. For 50 years, a central puzzle has been that observed $\beta$-decay rates are systematically smaller than theoretical predictions. This was attributed to an apparent quenching of the fundamental coupling constant $g_A \simeq $ 1.27 in the nucleus by a factor of about 0.75 compared to the $\beta$-decay of a free neutron. The origin of this quenching is controversial and has so far eluded a first-principles theoretical understanding. Here we address this puzzle and show that this quenching arises to a large extent from the coupling of the weak force to two nucleons as well as from strong correlations in the nucleus. We present state-of-the-art computations of $\beta$-decays from light to heavy nuclei. Our results are consistent with experimental data, including the pioneering measurement for $^{100}$Sn. These theoretical advances are enabled by systematic effective field theories of the strong and weak interactions combined with powerful quantum many-body techniques. This work paves the way for systematic theoretical predictions for fundamental physics problems. These include the synthesis of heavy elements in neutron star mergers and the search for neutrino-less double-$\beta$-decay, where an analogous quenching puzzle is a major source of uncertainty in extracting the neutrino mass scale., Comment: 20 pages, 18 figures

A thick neutron skin emerges from the first determination of root mean square radii of the proton distributions for $^{17-22}$N from charge changing cross section measurements around 900$A$ MeV at GSI. Neutron halo effects are signaled for $^{22}$N from an increase in the proton and matter radii. The radii suggest an unconventional shell gap at $N$ = 14 arising from the attractive proton-neutron tensor interaction, in good agreement with shell model calculations. $Ab$ $initio$, in-medium similarity re-normalization group, calculations with a state-of-the-art chiral nucleon-nucleon and three-nucleon interaction reproduce well the data approaching the neutron drip-line isotopes but are challenged in explaining the complete isotopic trend of the radii., Comment: 18 pages, 3 figures

The first $\gamma$-ray spectroscopy of $^{52}$Ar, with the neutron number N = 34, was measured using the $^{53}$K(p,2p) one-proton removal reaction at $\sim$210 MeV/u at the RIBF facility. The 2$^{+}_{1}$ excitation energy is found at 1656(18) keV, the highest among the Ar isotopes with N $>$ 20. This result is the first experimental signature of the persistence of the N = 34 subshell closure beyond $^{54}$Ca, i.e., below the magic proton number Z = 20. Shell-model calculations with phenomenological and chiral-effective-field-theory interactions both reproduce the measured 2$^{+}_{1}$ systematics of neutron-rich Ar isotopes, and support a N = 34 subshell closure in $^{52}$Ar.

Knowledge of the optical turbulence profile is important in adaptive optics (AO) systems, particularly tomographic AO systems such as those to be employed by the next generation of 40 m class extremely large telescopes (ELTs). Site characterisation and monitoring campaigns have produced large quantities of turbulence profiling data for sites around the world. However AO system design and performance characterisation is dependent on Monte-Carlo simulations that cannot make use of these large datasets due to long computation times. Here we address the question of how to reduce these large datasets into small sets of profiles that can feasibly be used in such Monte-Carlo simulations, whilst minimising the loss of information inherent in this effective compression of the data. We propose hierarchical clustering to partition the dataset according to the structure of the turbulence profiles and extract a single profile from each cluster. This method is applied to the Stereo-SCIDAR dataset from ESO Paranal containing over 10000 measurements of the turbulence profile from 83 nights. We present two methods of extracting turbulence profiles from the clusters, resulting in two sets of 18 profiles providing subtly different descriptions of the variability across the entire dataset. For generality we choose integrated parameters of the turbulence to measure the representativeness of our profiles and compare to others. Using these criterion we also show that such variability is difficult to capture with small sets of profiles associated with integrated turbulence parameters such as seeing., Comment: Accepted MNRAS September 2018

We present an approach to derive effective shell-model interactions from microscopic nuclear forces. The similarity-transformed coupled-cluster Hamiltonian decouples the single-reference state of a closed-shell nucleus and provides us with a core for the shell model. We use a second similarity transformation to decouple a shell-model space from the excluded space. We show that the three-body terms induced by both similarity transformations are crucial for an accurate computation of ground and excited states. As a proof of principle we use a nucleon-nucleon interaction from chiral effective field theory, employ a $^4$He core, and compute low-lying states of $^{6-8}$He and $^{6-8}$Li in $p$-shell model spaces. Our results agree with benchmarks from full configuration interaction., Comment: 12 pages, 16 figures

HARMONI is a visible and NIR integral field spectrograph, providing the E-ELT's core spectroscopic capability at first light. HARMONI will work at the diffraction limit of the E-ELT, thanks to a Classical and a Laser Tomographic AO system. In this paper, we present the system choices that have been made for these SCAO and LTAO modules. In particular, we describe the strategy developed for the different Wave-Front Sensors: pyramid for SCAO, the LGSWFS concept, the NGSWFS path, and the truth sensor capabilities. We present first potential implementations. And we asses the first system performance., Comment: SPIE proceeding, 2016

Image slicing is a powerful technique in astronomy. It allows the instrument designer to reduce the slit width of the spectrograph, increasing spectral resolving power whilst retaining throughput. Conventionally this is done using bulk optics, such as mirrors and prisms, however more recently astrophotonic components known as photonic lanterns (PLs) and photonic reformatters have also been used. These devices reformat the multi-mode (MM) input light from a telescope into single-mode (SM) outputs, which can then be re-arranged to suit the spectrograph. The photonic dicer (PD) is one such device, designed to reduce the dependence of spectrograph size on telescope aperture and eliminate modal noise. We simulate the PD, by optimising the throughput and geometrical design using Soapy and BeamProp. The simulated device shows a transmission between 8 and 20 %, depending upon the type of adaptive optics (AO) correction applied, matching the experimental results well. We also investigate our idealised model of the PD and show that the barycentre of the slit varies only slightly with time, meaning that the modal noise contribution is very low when compared to conventional fibre systems. We further optimise our model device for both higher throughput and reduced modal noise. This device improves throughput by 6.4 % and reduces the movement of the slit output by 50%, further improving stability. This shows the importance of properly simulating such devices, including atmospheric effects. Our work complements recent work in the field and is essential for optimising future photonic reformatters.

We present a quantum-classical algorithm to study the dynamics of the two-spatial-site Schwinger model on IBM's quantum computers. Using rotational symmetries, total charge, and parity, the number of qubits needed to perform computation is reduced by a factor of $\sim 5$, removing exponentially-large unphysical sectors from the Hilbert space. Our work opens an avenue for exploration of other lattice quantum field theories, such as quantum chromodynamics, where classical computation is used to find symmetry sectors in which the quantum computer evaluates the dynamics of quantum fluctuations., Comment: 5 pages, 4 figures, 23 pages supplemental, 8 figures supplemental

We report a quantum simulation of the deuteron binding energy on quantum processors accessed via cloud servers. We use a Hamiltonian from pionless effective field theory at leading order. We design a low-depth version of the unitary coupled-cluster ansatz, use the variational quantum eigensolver algorithm, and compute the binding energy to within a few percent. Our work is the first step towards scalable nuclear structure computations on a quantum processor via the cloud, and it sheds light on how to map scientific computing applications onto nascent quantum devices.

The nuclear shell structure, which originates in the nearly independent motion of nucleons in an average potential, provides an important guide for our understanding of nuclear structure and the underlying nuclear forces. Its most remarkable fingerprint is the existence of the so-called `magic numbers' of protons and neutrons associated with extra stability. Although the introduction of a phenomenological spin-orbit (SO) coupling force in 1949 helped explain the nuclear magic numbers, its origins are still open questions. Here, we present experimental evidence for the smallest SO-originated magic number (subshell closure) at the proton number 6 in 13-20C obtained from systematic analysis of point-proton distribution radii, electromagnetic transition rates and atomic masses of light nuclei. Performing ab initio calculations on 14,15C, we show that the observed proton distribution radii and subshell closure can be explained by the state-of-the-art nuclear theory with chiral nucleon-nucleon and three-nucleon forces, which are rooted in the quantum chromodynamics., Comment: 7 pages, 5 figures

We link the structure of nuclei around $^{100}$Sn, the heaviest doubly magic nucleus with equal neutron and proton numbers ($N=Z=50$), to nucleon-nucleon ($NN$) and three-nucleon ($NNN$) forces constrained by data of few-nucleon systems. Our results indicate that $^{100}$Sn is doubly magic, and we predict its quadrupole collectivity. We present precise computations of $^{101}$Sn based on three-particle--two-hole excitations of $^{100}$Sn, and reproduce the small splitting between the lowest $J^\pi=7/2^+$ and $5/2^+$ states. Our results are consistent with the sparse available data., Comment: 8 pages, 4 figures

We construct a nuclear interaction in chiral effective field theory with explicit inclusion of the $\Delta$-isobar $\Delta(1232)$ degree of freedom at all orders up to next-to-next-to-leading order (NNLO). We use pion-nucleon ($\pi N$) low-energy constants (LECs) from a Roy-Steiner analysis of $\pi N$ scattering data, optimize the LECs in the contact potentials up to NNLO to reproduce low-energy nucleon-nucleon scattering phase shifts, and constrain the three-nucleon interaction at NNLO to reproduce the binding energy and point-proton radius of $^{4}$He. For heavier nuclei we use the coupled-cluster method to compute binding energies, radii, and neutron skins. We find that radii and binding energies are much improved for interactions with explicit inclusion of $\Delta(1232)$, while $\Delta$-less interactions produce nuclei that are not bound with respect to breakup into $\alpha$ particles. The saturation of nuclear matter is significantly improved, and its symmetry energy is consistent with empirical estimates., Comment: 13 pages, 8 figures

Since fall 2012, several National Centers of Academic Excellence in Cyber Defense Research (CAE-Rs) fielded a collaborative course to engage students in solving applied cybersecurity research problems. We describe our experiences with this Information Security Research and Education (INSuRE) research collaborative. We explain how we conducted our project-based research course, give examples of student projects, and discuss the outcomes and lessons learned., Comment: A shorter version of this paper has been submitted to IEEE Security and Privacy

The performance of adaptive optics systems is partially dependant on the algorithms used within the real-time control system to compute wavefront slope measurements. We demonstrate use of a matched filter algorithm for the processing of elongated laser guide star (LGS) Shack-Hartmann images, using the CANARY adaptive optics instrument on the 4.2m William Herschel Telescope and the European Southern Observatory Wendelstein LGS Unit placed 40m away. This algorithm has been selected for use with the forthcoming Thirty Meter Telescope, but until now had not been demonstrated on-sky. From the results of a first observing run, we show that the use of matched filtering improves our adaptive optics system performance, with increases in on-sky H-band Strehl measured up to about a factor of 1.1 with respect to a conventional centre of gravity approach. We describe the algorithm used, and the methods that we implemented to enable on-sky demonstration., Comment: Accepted for publication in MNRAS

As telescopes become larger, into the era of ~40 m Extremely Large Telescopes, the high- resolution vertical profile of the optical turbulence strength is critical for the validation, optimization and operation of optical systems. The velocity of atmospheric optical turbulence is an important parameter for several applications including astronomical adaptive optics systems. Here, we compare the vertical profile of the velocity of the atmospheric wind above La Palma by means of a comparison of Stereo-SCIntillation Detection And Ranging (Stereo- SCIDAR) with the Global Forecast System models and nearby balloon-borne radiosondes. We use these data to validate the automated optical turbulence velocity identification from the Stereo-SCIDAR instrument mounted on the 2.5 m Isaac Newton Telescope, La Palma. By comparing these data we infer that the turbulence velocity and the wind velocity are consistent and that the automated turbulence velocity identification of the Stereo-SCIDAR is precise. The turbulence velocities can be used to increase the sensitivity of the turbulence strength profiles, as weaker turbulence that may be misinterpreted as noise can be detected with a velocity vector. The turbulence velocities can also be used to increase the altitude resolution of a detected layer, as the altitude of the velocity vectors can be identified to a greater precision than the native resolution of the system. We also show examples of complex velocity structure within a turbulent layer caused by wind shear at the interface of atmospheric zones., Comment: Published in MNRAS, 464, 3998-4007, 2017

We present a pedagogical discussion of Similarity Renormalization Group (SRG) methods, in particular the In-Medium SRG (IMSRG) approach for solving the nuclear many-body problem. These methods use continuous unitary transformations to evolve the nuclear Hamiltonian to a desired shape. The IMSRG, in particular, is used to decouple the ground state from all excitations and solve the many-body Schr\"odinger equation. We discuss the IMSRG formalism as well as its numerical implementation, and use the method to study the pairing model and infinite neutron matter. We compare our results with those of Coupled cluster theory, Configuration-Interaction Monte Carlo, and the Self-Consistent Green's Function approach. The chapter concludes with an expanded overview of current research directions, and a look ahead at upcoming developments., Comment: 73 pages, to appear in Lecture Notes in Physics (Springer), "An advanced course in computational nuclear physics: Bridging the scales from quarks to neutron stars", edited by M. Hjorth-Jensen, M. P. Lombardo, and U. van Kolck

We present two new methods for performing ab initio calculations of excited states for closed-shell systems within the in-medium similarity renormalization group (IMSRG) framework. Both are based on combining the IMSRG with simple many-body methods commonly used to target excited states, such as the Tamm-Dancoff approximation (TDA) and equations-of-motion (EOM) techniques. In the first approach, a two-step sequential IMSRG transformation is used to drive the Hamiltonian to a form where a simple TDA calculation (i.e., diagonalization in the space of $1$p$1$h excitations) becomes exact for a subset of eigenvalues. In the second approach, equations-of-motion (EOM) techniques are applied to the ground-state-decoupled IMSRG Hamiltonian to access excited states. We perform proof-of-principle calculations for parabolic quantum dots in two-dimensions and the closed shell nuclei $^{16}$O and $^{22}$O. We find that the TDA-IMSRG approach gives better accuracy than the EOM-IMSRG when calculations converge, but is otherwise lacking the versatility and numerical stability of the latter. Our calculated spectra are in reasonable agreement with analogous EOM-coupled-cluster (EOM-CC) calculations, which paves the way for more interesting applications of the EOM-IMSRG to calculations of consistently evolved observables such as electromagnetic strength functions and nuclear matrix elements, and extensions to nuclei within 1-2 nucleons of a closed shell by generalizing the EOM ladder operator to include particle-number nonconserving terms., Comment: 16 pages, 13 figures

Canary is the Multi-Object Adaptive Optics (MOAO) pathfinder for the future MOAO-assisted Integral-Field Units (IFU) proposed for Extremely Large Telescopes (ELT). The MOAO concept relies on tomographically reconstructing the turbulence using multiple measurements along different lines of sight. Tomography requires the knowledge of the statistical turbulence parameters, commonly recovered from the system telemetry using a dedicated profiling technique. For demonstration purposes with the MOAO pathfinder Canary , this identification is performed thanks to the Learn & Apply (L&A) algorithm, that consists in model- fitting the covariance matrix of WFS measurements dependent on relevant parameters: $C_n^2(h)$ profile, outer scale profile and system mis-registration. We explore an upgrade of this algorithm, the Learn 3 Steps (L3S) approach, that allows one to dissociate the identification of the altitude layers from the ground in order to mitigate the lack of convergence of the required empirical covariance matrices therefore reducing the required length of data time-series for reaching a given accuracy. For nominal observation conditions, the L3S can reach the same level of tomographic error in using five times less data frames than the L&A approach. The L3S technique has been applied over a large amount of Canary data to characterize the turbu- lence above the William Herschel Telescope (WHT). These data have been acquired the 13th, 15th, 16th, 17th and 18th September 2013 and we find 0.67"/8.9m/3.07m/s of total seeing/outer scale/wind-speed, with 0.552"/9.2m/2.89m/s below 1.5 km and 0.263"/10.3m/5.22m/s between 1.5 and 20 km. We have also de- termined the high altitude layers above 20 km, missed by the tomographic reconstruction on Canary , have a median seeing of 0.187" and have occurred 16% of observation time.

CANARY is the multi-object adaptive optics (MOAO) on-sky pathfinder developed in the perspective of Multi-Object Spectrograph on Extremely Large Telescopes~(ELTs). In 2013, CANARY was operated on-sky at the William Herschel telescope~(WHT), using three off-axis natural guide stars~(NGS) and four off-axis Rayleigh laser guide stars~(LGS), in open-loop, with the on-axis compensated turbulence observed with a H-band imaging camera and a Truth wave-front sensor~(TS) for diagnostic purposes. Our purpose is to establish a reliable and accurate wave-front error breakdown for LGS MOAO. This will enable a comprehensive analysis of \cana on-sky results and provide tools for validating simulations of MOAO systems for ELTs. To evaluate the MOAO performance, we compared the CANARY on-sky results running in MOAO, in Single Conjugated Adaptive Optics~(SCAO) and in Ground Layer Adaptive Optics~(GLAO) modes, over a large set of data acquired in 2013. We provide a statistical study of the seeing. We also evaluated the wave-front error breakdown from both analytic computations, one based on a MOAO system modelling and the other on the measurements from the CANARY TS. We have focussed especially on the tomographic error and we detail its vertical error decomposition~(VED). We show that CANARY obtained 30.1\%, 21.4\% and 17.1\% H-band Strehl ratios in SCAO, MOAO and GLAO respectively, for median seeing conditions with 0.66" of total seeing including 0.59" at the ground. Moreover, we get 99\% of correlation over 4,500 samples, for any AO modes, between two analytic computations of residual phase variance. Based on these variances, we obtain a reasonable Strehl-ratio~(SR) estimation when compared to the measured IR image SR. We evaluate the gain in compensation for the altitude turbulence brought by MOAO when compared to GLAO.

In preparation of future Multi-Object Spectrographs (MOS) whose one of the major role is to provide an extensive statistical studies of high redshifted galaxies surveyed, the demonstrator Canary has been designed to tackle technical challenges related to open-loop Adaptive-Optics (AO) control with jointed Natural Guide Star (NGS) and Laser Guide Star (LGS) tomography. We have developed a Point Spread Function (PSF)-Reconstruction algorithm dedicated to MOAO systems using system telemetry to estimate the PSF potentially anywhere in the observed field, a prerequisite to post- process AO-corrected observations in Integral Field Spectroscopy (IFS). In this paper we show how to handle off-axis data to estimate the PSF using atmospheric tomography and compare it to a classical approach that uses on-axis residual phase from a truth sensor observing a natural bright source. We have reconstructed over 450 on-sky Canary PSFs and we get bias/1-${\sigma}$ standard-deviation (std) of 1.3/4.8 on the H-band Strehl ratio (SR) with 92.3% of correlation between reconstructed and sky SR. On the Full Width at Half Maximum (FWHM), we get respectively 2.94 mas, 19.9 mas and 88.3% for the bias, std and correlation. The reference method achieves 0.4/3.5/95% on the SR and 2.71 mas/14.9 mas/92.5% on the FWHM for the bias/std/correlation.

The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study., Comment: 12 pages, in Ground-based and Airborne Instrumentation for Astronomy VI, 2016, Proc. SPIE 9908, 23

Building on the comprehensive White Paper on the scientific case for multi-object spectroscopy on the European ELT, we present the top-level instrument requirements that are being used in the Phase A design study of the MOSAIC concept. The assembled cases span the full range of E-ELT science and generally require either 'high multiplex' or 'high definition' observations to best exploit the excellent sensitivity and spatial performance of the telescope. We highlight some of the science studies that are now being used in trade-off studies to inform the capabilities of MOSAIC and its technical design., Comment: 12 pages, in Ground-based and Airborne Instrumentation for Astronomy VI, 2016, Proc. SPIE 9908, 9J

HiPERCAM is a high-speed camera for the study of rapid variability in the Universe. The project is funded by a 3.5MEuro European Research Council Advanced Grant. HiPERCAM builds on the success of our previous instrument, ULTRACAM, with very significant improvements in performance thanks to the use of the latest technologies. HiPERCAM will use 4 dichroic beamsplitters to image simultaneously in 5 optical channels covering the u'g'r'i'z' bands. Frame rates of over 1000 per second will be achievable using an ESO CCD controller (NGC), with every frame GPS timestamped. The detectors are custom-made, frame-transfer CCDs from e2v, with 4 low-noise (2.5e-) outputs, mounted in small thermoelectrically-cooled heads operated at 180 K, resulting in virtually no dark current. The two reddest CCDs will be deep-depletion devices with anti-etaloning, providing high quantum efficiencies across the red part of the spectrum with no fringing. The instrument will also incorporate scintillation noise correction via the conjugate-plane photometry technique. The opto-mechanical chassis will make use of additive manufacturing techniques in metal to make a light-weight, rigid and temperature-invariant structure. First light is expected on the 4.2m William Herschel Telescope on La Palma in 2017 (on which the field of view will be 10' with a 0.3"/pixel scale), with subsequent use planned on the 10.4m Gran Telescopio Canarias on La Palma (on which the field of view will be 4' with a 0.11"/pixel scale) and the 3.5m New Technology Telescope in Chile., Comment: 9 pages, 8 figures, 1 table, SPIE Astronomical Telescopes and Instrumentation, Edinburgh (26 June - 1 July, 2016)

Recent advances in adaptive optics (AO) have led to the implementation of wide field-of-view AO systems. A number of wide-field AO systems are also planned for the forthcoming Extremely Large Telescopes. Such systems have multiple wavefront sensors of different types, and usually multiple deformable mirrors (DMs). Here, we report on our experience integrating cameras and DMs with the real-time control systems of two wide-field AO systems. These are CANARY, which has been operating on-sky since 2010, and DRAGON, which is a laboratory adaptive optics real-time demonstrator instrument. We detail the issues and difficulties that arose, along with the solutions we developed. We also provide recommendations for consideration when developing future wide-field AO systems., Comment: Accepted for publication in MNRAS

Laser Guide Stars (LGS) have greatly increased the sky-coverage of Adaptive Optics (AO) systems. Due to the up-link turbulence experienced by LGSs, a Natural Guide Star (NGS) is still required, preventing full sky-coverage. We present a method of obtaining partial tip-tilt information from LGSs alone in multi-LGS tomographic LGS AO systems. The method of LGS up-link tip-tilt determination is derived using a geometric approach, then an alteration to the Learn and Apply algorithm for tomographic AO is made to accommodate up-link tip-tilt. Simulation results are presented, verifying that the technique shows good performance in correcting high altitude tip-tilt, but not that from low altitudes. We suggest that the method is combined with multiple far off-axis tip-tilt NGSs to provide gains in performance and sky-coverage over current tomographic AO systems., Comment: Accepted by MNRAS, 11th March 2016

We present a comprehensive review of the In-Medium Similarity Renormalization Group (IM-SRG), a novel ab inito method for nuclei. The IM-SRG employs a continuous unitary transformation of the many-body Hamiltonian to decouple the ground state from all excitations, thereby solving the many-body problem. Starting from a pedagogical introduction of the underlying concepts, the IM-SRG flow equations are developed for systems with and without explicit spherical symmetry. We study different IM-SRG generators that achieve the desired decoupling, and how they affect the details of the IM-SRG flow. Based on calculations of closed-shell nuclei, we assess possible truncations for closing the system of flow equations in practical applications, as well as choices of the reference state. We discuss the issue of center-of-mass factorization and demonstrate that the IM-SRG ground-state wave function exhibits an approximate decoupling of intrinsic and center-of-mass degrees of freedom, similar to Coupled Cluster (CC) wave functions. To put the IM-SRG in context with other many-body methods, in particular many-body perturbation theory and non-perturbative approaches like CC, a detailed perturbative analysis of the IM-SRG flow equations is carried out. We conclude with a discussion of ongoing developments, including IM-SRG calculations with three-nucleon forces, the multi-reference IM-SRG for open-shell nuclei, first non-perturbative derivations of shell- model interactions, and the consistent evolution of operators in the IM-SRG. We dedicate this review to the memory of Gerry Brown, one of the pioneers of many-body calculations of nuclei., Comment: 92 pages, 33 figures, to appear in Physics Reports

From the Wilsonian point of view, renormalisable theories are understood as submanifolds in theory space emanating from a particular fixed point under renormalisation group evolution. We show how this picture precisely applies to their gravity duals. We investigate the Hamilton-Jacobi equation satisfied by the Wilson action and find the corresponding fixed points and their eigendeformations, which have a diagonal evolution close to the fixed points. The relevant eigendeformations are used to construct renormalised theories. We explore the relation of this formalism with holographic renormalisation. We also discuss different renormalisation schemes and show that the solutions to the gravity equations of motion can be used as renormalised couplings that parametrise the renormalised theories. This provides a transparent connection between holographic renormalisation group flows in the Wilsonian and non-Wilsonian approaches. The general results are illustrated by explicit calculations in an interacting scalar theory in AdS space., Comment: 63 pages. Minor changes and references added. Matches JHEP version

We present an improved variant of the in-medium similarity renormalization group (IM-SRG) based on the Magnus expansion. In the new formulation, one solves flow equations for the anti-hermitian operator that, upon exponentiation, yields the unitary transformation of the IM-SRG. The resulting flow equations can be solved using a first-order Euler method without any loss of accuracy, resulting in substantial memory savings and modest computational speedups. Since one obtains the unitary transformation directly, the transformation of additional operators beyond the Hamiltonian can be accomplished with little additional cost, in sharp contrast to the standard formulation of the IM-SRG. Ground state calculations of the homogeneous electron gas (HEG) and $^{16}$O nucleus are used as test beds to illustrate the efficacy of the Magnus expansion., Comment: 12 pages, 9 figures; fixed typos and added a reference