-
Cost-efficient Active Illumination Camera For Hyper-spectral Reconstruction
Authors:
Yuxuan Zhang,
T. M. Sazzad,
Yangyang Song,
Spencer J. Chang,
Ritesh Chowdhry,
Tomas Mejia,
Anna Hampton,
Shelby Kucharski,
Stefan Gerber,
Barry Tillman,
Marcio F. R. Resende,
William M. Hammond,
Chris H. Wilson,
Alina Zare,
Sanjeev J. Koppal
Abstract:
Hyper-spectral imaging has recently gained increasing attention for use in different applications, including agricultural investigation, ground tracking, remote sensing and many other. However, the high cost, large physical size and complicated operation process stop hyperspectral cameras from being employed for various applications and research fields. In this paper, we introduce a cost-efficient…
▽ More
Hyper-spectral imaging has recently gained increasing attention for use in different applications, including agricultural investigation, ground tracking, remote sensing and many other. However, the high cost, large physical size and complicated operation process stop hyperspectral cameras from being employed for various applications and research fields. In this paper, we introduce a cost-efficient, compact and easy to use active illumination camera that may benefit many applications. We developed a fully functional prototype of such camera. With the hope of helping with agricultural research, we tested our camera for plant root imaging. In addition, a U-Net model for spectral reconstruction was trained by using a reference hyperspectral camera's data as ground truth and our camera's data as input. We demonstrated our camera's ability to obtain additional information over a typical RGB camera. In addition, the ability to reconstruct hyperspectral data from multi-spectral input makes our device compatible to models and algorithms developed for hyperspectral applications with no modifications required.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
Scalable Optimal Motion Planning for Multi-Agent Systems by Cosserat Theory of Rods
Authors:
Amirreza Fahim Golestaneh,
Maxwell Hammond,
Venanzio Cichella
Abstract:
We address the motion planning problem for large multi-agent systems, utilizing Cosserat rod theory to model the dynamic behavior of vehicle formations. The problem is formulated as an optimal control problem over partial differential equations (PDEs) that describe the system as a continuum. This approach ensures scalability with respect to the number of vehicles, as the problem's complexity remai…
▽ More
We address the motion planning problem for large multi-agent systems, utilizing Cosserat rod theory to model the dynamic behavior of vehicle formations. The problem is formulated as an optimal control problem over partial differential equations (PDEs) that describe the system as a continuum. This approach ensures scalability with respect to the number of vehicles, as the problem's complexity remains unaffected by the size of the formation. The numerical discretization of the governing equations and problem's constraints is achieved through Bernstein surface polynomials, facilitating the conversion of the optimal control problem into a nonlinear programming (NLP) problem. This NLP problem is subsequently solved using off-the-shelf optimization software. We present several properties and algorithms related to Bernstein surface polynomials to support the selection of this methodology. Numerical demonstrations underscore the efficacy of this mathematical framework.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
Simultaneous retrieval of orbital phase resolved JWST/MIRI emission spectra of the hot Jupiter WASP-43b: evidence of water, ammonia and carbon monoxide
Authors:
Jingxuan Yang,
Mark Hammond,
Anjali A. A. Piette,
Jasmina Blecic,
Taylor J. Bell,
Patrick G. J. Irwin,
Vivien Parmentier,
Shang-Min Tsai,
Joanna K. Barstow,
Nicolas Crouzet,
Laura Kreidberg,
João M. Mendonça,
Jake Taylor,
Robin Baeyens,
Kazumasa Ohno,
Lucas Teinturier,
Matthew C. Nixon
Abstract:
Spectroscopic phase curves of hot Jupiters measure their emission spectra at multiple orbital phases, thus enabling detailed characterisation of their atmospheres. Precise constraints on the atmospheric composition of these exoplanets offer insights into their formation and evolution. We analyse four phase-resolved emission spectra of the hot Jupiter WASP-43b, generated from a phase curve observed…
▽ More
Spectroscopic phase curves of hot Jupiters measure their emission spectra at multiple orbital phases, thus enabling detailed characterisation of their atmospheres. Precise constraints on the atmospheric composition of these exoplanets offer insights into their formation and evolution. We analyse four phase-resolved emission spectra of the hot Jupiter WASP-43b, generated from a phase curve observed with the MIRI/LRS onboard the JWST, to retrieve its atmospheric properties. Using a parametric 2D temperature model and assuming a chemically homogeneous atmosphere within the observed pressure region, we simultaneously fit the four spectra to constrain the abundances of atmospheric constituents, thereby yielding more precise constraints than previous work that analysed each spectrum independently. Our analysis reveals statistically significant evidence of NH3 (4$σ$) in a hot Jupiter's emission spectra for the first time, along with evidence of H2O (6.5$σ$), CO (3.1$σ$), and a non-detection of CH4. With our abundance constraints, we tentatively estimate the metallicity of WASP-43b at 0.6-6.5$\times$solar and its C/O ratio at 0.6-0.9. Our findings offer vital insights into the atmospheric conditions and formation history of WASP-43b by simultaneously constraining the abundances of carbon, oxygen, and nitrogen-bearing species.
△ Less
Submitted 5 June, 2024;
originally announced June 2024.
-
Identifying and Fitting Eclipse Maps of Exoplanets with Cross-Validation
Authors:
Mark Hammond,
Neil T. Lewis,
Sasha Boone,
Xueqing Chen,
João M. Mendonça,
Vivien Parmentier,
Jake Taylor,
Taylor Bell,
Leonardo dos Santos,
Nicolas Crouzet,
Laura Kreidberg,
Michael Radica,
Michael Zhang
Abstract:
Eclipse mapping uses the shape of the eclipse of an exoplanet to measure its two-dimensional structure. Light curves are mostly composed of longitudinal information, with the latitudinal information only contained in the brief ingress and egress of the eclipse. This imbalance can lead to a spuriously confident map, where the longitudinal structure is constrained by out-of-eclipse data and the lati…
▽ More
Eclipse mapping uses the shape of the eclipse of an exoplanet to measure its two-dimensional structure. Light curves are mostly composed of longitudinal information, with the latitudinal information only contained in the brief ingress and egress of the eclipse. This imbalance can lead to a spuriously confident map, where the longitudinal structure is constrained by out-of-eclipse data and the latitudinal structure is wrongly determined by the priors on the map. We present a new method to address this issue. The method tests for the presence of an eclipse mapping signal by using k-fold cross-validation to compare the performance of a simple mapping model to the null hypothesis of a uniform disk. If a signal is found, the method fits a map with more degrees of freedom, optimising its information content. The information content is varied by penalising the model likelihood by a factor proportional to the spatial entropy of the map, optimised by cross-validation. We demonstrate this method for simulated datasets then apply it to three observational datasets. The method identifies an eclipse mapping signal for JWST MIRI/LRS observations of WASP-43b but does not identify a signal for JWST NIRISS/SOSS observations of WASP-18b or Spitzer Space Telescope observations of HD 189733b. It is possible to fit eclipse maps to these datasets, but we suggest that these maps are overfitting the eclipse shape. We fit a new map with more spatial freedom to the WASP-43b dataset and show a flatter east-west structure than previously derived.
△ Less
Submitted 31 May, 2024;
originally announced May 2024.
-
Two-Dimensional Eclipse Mapping of the Hot Jupiter WASP-43b with JWST MIRI/LRS
Authors:
Mark Hammond,
Taylor J. Bell,
Ryan C. Challener,
Neil T. Lewis,
Megan Weiner Mansfield,
Isaac Malsky,
Emily Rauscher,
Jacob L. Bean,
Ludmila Carone,
João M. Mendonça,
Lucas Teinturier,
Xianyu Tan,
Nicolas Crouzet,
Laura Kreidberg,
Giuseppe Morello,
Vivien Parmentier,
Jasmina Blecic,
Jean-Michel Désert,
Christiane Helling,
Pierre-Olivier Lagage,
Karan Molaverdikhani,
Matthew C. Nixon,
Benjamin V. Rackham,
Jingxuan Yang
Abstract:
We present eclipse maps of the two-dimensional thermal emission from the dayside of the hot Jupiter WASP-43b, derived from an observation of a phase curve with the JWST MIRI/LRS instrument. The observed eclipse shapes deviate significantly from those expected for a planet emitting uniformly over its surface. We fit a map to this deviation, constructed from spherical harmonics up to order…
▽ More
We present eclipse maps of the two-dimensional thermal emission from the dayside of the hot Jupiter WASP-43b, derived from an observation of a phase curve with the JWST MIRI/LRS instrument. The observed eclipse shapes deviate significantly from those expected for a planet emitting uniformly over its surface. We fit a map to this deviation, constructed from spherical harmonics up to order $\ell_{\rm max}=2$, alongside the planetary, orbital, stellar, and systematic parameters. This yields a map with a meridionally-averaged eastward hot-spot shift of $(7.75 \pm 0.36)^{\circ}$, with no significant degeneracy between the map and the additional parameters. We show the latitudinal and longitudinal contributions of the day-side emission structure to the eclipse shape, finding a latitudinal signal of $\sim$200 ppm and a longitudinal signal of $\sim$250 ppm. To investigate the sensitivity of the map to the method, we fix the non-mapping parameters and derive an "eigenmap" fitted with an optimised number of orthogonal phase curves, which yields a similar map to the $\ell_{\rm max}=2$ map. We also fit a map up to $\ell_{\rm max}=3$, which shows a smaller hot-spot shift, with a larger uncertainty. These maps are similar to those produced by atmospheric simulations. We conclude that there is a significant mapping signal which constrains the spherical harmonic components of our model up to $\ell_{\rm max}=2$. Alternative mapping models may derive different structures with smaller-scale features; we suggest that further observations of WASP-43b and other planets will drive the development of more robust methods and more accurate maps.
△ Less
Submitted 25 April, 2024;
originally announced April 2024.
-
Hot Jupiter Diversity and the Onset of TiO/VO Revealed by a Large Grid of Non-Grey Global Circulation Models
Authors:
Alexander Roth,
Vivien Parmentier,
Mark Hammond
Abstract:
The population of hot Jupiters is extremely diverse, with large variations in their irradiation, period, gravity and chemical composition. To understand the intrinsic planet diversity through the observed population level trends, we explore the a-priori scatter in the population created by the different responses of atmospheric circulation to planetary parameters. We use the SPARC/MITgcm 3D global…
▽ More
The population of hot Jupiters is extremely diverse, with large variations in their irradiation, period, gravity and chemical composition. To understand the intrinsic planet diversity through the observed population level trends, we explore the a-priori scatter in the population created by the different responses of atmospheric circulation to planetary parameters. We use the SPARC/MITgcm 3D global circulation model to simulate 345 planets spanning a wide range of instellation, metallicity, gravity and rotation periods typical for hot Jupiters, while differentiating between models with and without TiO/VO in their atmosphere. We show that the combined effect of the planetary parameters leads to a large diversity in the ability of atmospheres to transport heat from day-side to night-side at a given equilibrium temperature. We further show that the hot-spot offset is a non-monotonic function of planetary rotation period and explain our findings by a competition between the rotational and divergent parts of the circulation. As a consequence, hot-spot offset and phase curve amplitude are not necessarily correlated. Finally, we compare the observables from our grid to the population of Spitzer and Hubble observations of hot Jupiters. We find that the sudden jump in brightness temperature observed in the Spitzer secondary eclipse measurements can be naturally explained by the cold-trapping of TiO/VO at approximately 1800K. The grid of modelled spectra, phase curves and thermal structures are made available to the community, together with a python code for visualization of the grid properties, at https://doi.org/10.5281/zenodo.10785321 and http://sim3d.oca.eu/.
△ Less
Submitted 15 April, 2024;
originally announced April 2024.
-
Fabrication Tolerant Multi-Layer Integrated Photonic Topology Optimization
Authors:
Michael J. Probst,
Arjun Khurana,
Joel B. Slaby,
Alec M. Hammond,
Stephen E. Ralph
Abstract:
Optimal multi-layer device design requires consideration of fabrication uncertainties associated with inter-layer alignment and conformal layering. We present layer-restricted topology optimization (TO), a novel technique which mitigates the effects of unwanted conformal layering for multi-layer structures and enables TO in multi-etch material platforms. We explore several approaches to achieve th…
▽ More
Optimal multi-layer device design requires consideration of fabrication uncertainties associated with inter-layer alignment and conformal layering. We present layer-restricted topology optimization (TO), a novel technique which mitigates the effects of unwanted conformal layering for multi-layer structures and enables TO in multi-etch material platforms. We explore several approaches to achieve this result compatible with density-based TO projection techniques and geometric constraints. Then, we present a robust TO formulation to design devices resilient to inter-layer misalignment. The novel constraint and robust formulation are demonstrated in 2D grating couplers and a 3D polarization rotator.
△ Less
Submitted 10 April, 2024;
originally announced April 2024.
-
Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b
Authors:
Taylor J. Bell,
Nicolas Crouzet,
Patricio E. Cubillos,
Laura Kreidberg,
Anjali A. A. Piette,
Michael T. Roman,
Joanna K. Barstow,
Jasmina Blecic,
Ludmila Carone,
Louis-Philippe Coulombe,
Elsa Ducrot,
Mark Hammond,
João M. Mendonça,
Julianne I. Moses,
Vivien Parmentier,
Kevin B. Stevenson,
Lucas Teinturier,
Michael Zhang,
Natalie M. Batalha,
Jacob L. Bean,
Björn Benneke,
Benjamin Charnay,
Katy L. Chubb,
Brice-Olivier Demory,
Peter Gao
, et al. (58 additional authors not shown)
Abstract:
Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5…
▽ More
Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5-12 $μ$m with JWST's Mid-Infrared Instrument (MIRI). The spectra reveal a large day-night temperature contrast (with average brightness temperatures of 1524$\pm$35 and 863$\pm$23 Kelvin, respectively) and evidence for water absorption at all orbital phases. Comparisons with three-dimensional atmospheric models show that both the phase curve shape and emission spectra strongly suggest the presence of nightside clouds which become optically thick to thermal emission at pressures greater than ~100 mbar. The dayside is consistent with a cloudless atmosphere above the mid-infrared photosphere. Contrary to expectations from equilibrium chemistry but consistent with disequilibrium kinetics models, methane is not detected on the nightside (2$σ$ upper limit of 1-6 parts per million, depending on model assumptions).
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
Modeling, Simulation, and Maneuvering Control of a Generic Submarine
Authors:
Gage MacLin,
Maxwell Hammond,
Venanzio Cichella,
J. Ezequiel Martin
Abstract:
This work introduces two multi-level control strategies to address the problem of guidance and control of underwater vehicles. An outer-loop path-following algorithm and an outer-loop trajectory tracking algorithm are presented. Both outer-loop algorithms provide reference commands that enable the generic submarine to adhere to a three-dimensional path, and both use an inner-loop adaptive controll…
▽ More
This work introduces two multi-level control strategies to address the problem of guidance and control of underwater vehicles. An outer-loop path-following algorithm and an outer-loop trajectory tracking algorithm are presented. Both outer-loop algorithms provide reference commands that enable the generic submarine to adhere to a three-dimensional path, and both use an inner-loop adaptive controller to determine the required actuation commands. Further, a reduced order model of a generic submarine is presented. Computational fluid dynamics (CFD) results are used to create and validate a model that includes depth dependence and the effect of waves on the craft. %The model and the procedure to obtain its coefficients are discussed, and examples of the data used to obtain the model coefficients are presented. An example of operation following a complex path is presented and Results from the reduced order model for each control strategy are compared.
△ Less
Submitted 18 December, 2023;
originally announced January 2024.
-
Path Planning for Continuum Rods Using Bernstein Surfaces
Authors:
Maxwell Hammond,
Venanzio Cichella,
Amirreza F. Golestaneh,
Caterina Lamuta
Abstract:
This paper presents a method for optimal motion planning of continuum robots by employing Bernstein surfaces to approximate the system's dynamics and impose complex constraints, including collision avoidance. The main contribution is the approximation of infinite-dimensional continuous problems into their discrete counterparts, facilitating their solution using standard optimization solvers. This…
▽ More
This paper presents a method for optimal motion planning of continuum robots by employing Bernstein surfaces to approximate the system's dynamics and impose complex constraints, including collision avoidance. The main contribution is the approximation of infinite-dimensional continuous problems into their discrete counterparts, facilitating their solution using standard optimization solvers. This discretization leverages the unique properties of Bernstein surface, providing a framework that extends previous works which focused on ODEs approximated by Bernstein polynomials. Numerical validations are conducted through several numerical scenarios. The presented methodology offers a promising direction for solving complex optimal control problems in the realm of soft robotics.
△ Less
Submitted 19 December, 2023;
originally announced December 2023.
-
Bioinspired Soft Robotics: state of the art, challenges, and future directions
Authors:
Maxwell Hammond,
Venanzio Cichella,
Caterina Lamuta
Abstract:
Purpose of Review: This review provides an overview of the state of the art in bioinspired soft robotics with by examining advancements in actuation, functionality, modeling, and control. Recent Findings: Recent research into actuation methods, such as artificial muscles, have expanded the functionality and potential use of bioinspired soft robots. Additionally, the application of finite dimension…
▽ More
Purpose of Review: This review provides an overview of the state of the art in bioinspired soft robotics with by examining advancements in actuation, functionality, modeling, and control. Recent Findings: Recent research into actuation methods, such as artificial muscles, have expanded the functionality and potential use of bioinspired soft robots. Additionally, the application of finite dimensional models has improved computational efficiency for modeling soft continuum systems, and garnered interest as a basis for controller formulation. Summary: Bioinspiration in the field of soft robotics has led to diverse approaches to problems in a range of task spaces. In particular, new capabilities in system simplification, miniaturization, and untethering have each contributed to the field's growth. There is still significant room for improvement in the streamlining of design and manufacturing for these systems, as well as in their control.
△ Less
Submitted 19 December, 2023;
originally announced December 2023.
-
Global Chemical Transport on Hot Jupiters: Insights from 2D VULCAN photochemical model
Authors:
Shang-Min Tsai,
Vivien Parmentier,
João M. Mendonça,
Xianyu Tan,
Russell Deitrick,
Mark Hammond,
Arjun B. Savel,
Xi Zhang,
Raymond T. Pierrehumbert,
Edward W. Schwieterman
Abstract:
The atmospheric dynamics of tidally-locked hot Jupiters is characterized by strong equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photoc…
▽ More
The atmospheric dynamics of tidally-locked hot Jupiters is characterized by strong equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photochemical (horizontal and vertical) transport model, VULCAN 2D, which improves on the pseudo-2D approaches by allowing for non-uniform zonal winds. We extensively validate our VULCAN 2D with analytical solutions and benchmark comparisons. Applications to HD 189733 b and HD 209458 b reveal a transition in mixing regimes: horizontal transport predominates below $\sim$0.1 mbar while vertical mixing is more important at higher altitudes above 0.1 mbar. Motivated by the previously inferred carbon-rich atmosphere, we find that HD 209458 b with super-solar carbon-to-oxygen ratio (C/O) exhibits pronounced C$_2$H$_4$ absorption on the morning limb but not on the evening limb, owing to horizontal transport from the nightside. We discuss when a pseudo-2D approach is a valid assumption and its inherent limitations. Finally, we demonstrate the effect of horizontal transport in transmission observations and its impact on the morning-evening limb asymmetry with synthetic spectra, highlighting the need to consider global transport when interpreting exoplanet atmospheres.
△ Less
Submitted 25 March, 2024; v1 submitted 26 October, 2023;
originally announced October 2023.
-
An analytic theory for the resolution attainable using eclipse mapping of exoplanets
Authors:
Sasha Boone,
David Grant,
Mark Hammond
Abstract:
We present an analytic theory for the resolution attainable via eclipse mapping of exoplanets, based on the Fourier components of the brightness distribution on the planetary disk. We find that the impact parameter determines which features can and cannot be seen, via the angle of the stellar edge relative to the axis of the orbit during the eclipse. We estimate the signal-to-noise ratio as a func…
▽ More
We present an analytic theory for the resolution attainable via eclipse mapping of exoplanets, based on the Fourier components of the brightness distribution on the planetary disk. We find that the impact parameter determines which features can and cannot be seen, via the angle of the stellar edge relative to the axis of the orbit during the eclipse. We estimate the signal-to-noise ratio as a function of mapping resolution, and use this to determine the attainable resolution for a given star-planet system. We test this theory against numerical simulations and find good agreement; in particular, our predictions for the resolution as a function of stellar edge angle are accurate to the simulated data to within 10% over a wide range of angles. Our prediction for the number of spatial modes that can be constrained given a light curve error is similarly accurate. Finally, we give a list of exoplanets with the best expected resolution for observations with the NIRISS SOSS, NIRSpec G395H, and MIRI LRS instruments on JWST.
△ Less
Submitted 22 October, 2023;
originally announced October 2023.
-
Fourier modal method for inverse design of metasurface-enhanced micro-LEDs
Authors:
Martin F. Schubert,
Alec M. Hammond
Abstract:
We present a simulation capability for micro-scale light-emitting diodes (uLEDs) that achieves comparable accuracy to CPU-based finite-difference time-domain simulation but is more than 10^7 times faster. Our approach is based on the Fourier modal method (FMM) -- which, as we demonstrate, is well suited to modeling thousands of incoherent sources -- with extensions that allow rapid convergence for…
▽ More
We present a simulation capability for micro-scale light-emitting diodes (uLEDs) that achieves comparable accuracy to CPU-based finite-difference time-domain simulation but is more than 10^7 times faster. Our approach is based on the Fourier modal method (FMM) -- which, as we demonstrate, is well suited to modeling thousands of incoherent sources -- with extensions that allow rapid convergence for uLED structures that are challenging to model with standard approaches. The speed of our method makes the inverse design of uLEDs tractable, which we demonstrate by designing a metasurface-enhanced uLED that doubles the light extraction efficiency of an unoptimized device.
△ Less
Submitted 15 August, 2023;
originally announced August 2023.
-
Interior dynamics of super-Earth 55 Cancri e
Authors:
Tobias G. Meier,
Dan J. Bower,
Tim Lichtenberg,
Mark Hammond,
Paul J. Tackley
Abstract:
The ultra-short-period super-Earth 55 Cancri e has a measured radius of 1.8 Earth radii. Previous thermal phase curve observations suggest a strong temperature contrast between the dayside and nightside of around 1000 K with the hottest point shifted $41\pm12$ degrees east from the substellar point, indicating some degree of heat circulation. The dayside (and potentially even the nightside) is hot…
▽ More
The ultra-short-period super-Earth 55 Cancri e has a measured radius of 1.8 Earth radii. Previous thermal phase curve observations suggest a strong temperature contrast between the dayside and nightside of around 1000 K with the hottest point shifted $41\pm12$ degrees east from the substellar point, indicating some degree of heat circulation. The dayside (and potentially even the nightside) is hot enough to harbour a magma ocean. We use results from general circulation models (GCMs) of atmospheres to constrain the surface temperature contrasts. There is still a large uncertainty on the vigour and style of mantle convection in super-Earths, especially those that experience stellar irradiation large enough to harbour a magma ocean. In this work, we aim to constrain the mantle dynamics of the tidally locked lava world 55 Cancri e. Using the surface temperature contrasts as boundary condition, we model the mantle flow of 55 Cancri e using 2D mantle convection simulations and investigate how the convection regimes are affected by the different climate models. We find that large super-plumes form on the dayside if that hemisphere is covered by a magma ocean and the nightside remains solid or only partially molten. Cold material descends into the deep interior on the nightside, but no strong downwellings form. In some cases, the super-plume also moves several tens of degrees towards the terminator. A convective regime where the upwelling is preferentially on the dayside might lead to preferential outgassing on that hemisphere which could lead to the build-up of atmospheric species that could be chemically distinct from the nightside.
△ Less
Submitted 1 August, 2023;
originally announced August 2023.
-
Examining NHD vs QHD in the GCM THOR with non-grey radiative transfer for the hot Jupiter regime
Authors:
Pascal A. Noti,
Elspeth K. H. Lee,
Russell Deitrick,
Mark Hammond
Abstract:
Global circulation models (GCMs) play an important role in contemporary investigations of exoplanet atmospheres. Different GCMs evolve various sets of dynamical equations which can result in obtaining different atmospheric properties between models. In this study, we investigate the effect of different dynamical equation sets on the atmospheres of hot Jupiter exoplanets. We compare GCM simulations…
▽ More
Global circulation models (GCMs) play an important role in contemporary investigations of exoplanet atmospheres. Different GCMs evolve various sets of dynamical equations which can result in obtaining different atmospheric properties between models. In this study, we investigate the effect of different dynamical equation sets on the atmospheres of hot Jupiter exoplanets. We compare GCM simulations using the quasi-primitive dynamical equations (QHD) and the deep Navier-Stokes equations (NHD) in the GCM THOR. We utilise a two-stream non-grey "picket-fence" scheme to increase the realism of the radiative transfer calculations. We perform GCM simulations covering a wide parameter range grid of system parameters in the population of exoplanets. Our results show significant differences between simulations with the NHD and QHD equation sets at lower gravity, higher rotation rates or at higher irradiation temperatures. The chosen parameter range shows the relevance of choosing dynamical equation sets dependent on system and planetary properties. Our results show the climate states of hot Jupiters seem to be very diverse, where exceptions to prograde superrotation can often occur. Overall, our study shows the evolution of different climate states which arise just due to different selections of Navier-Stokes equations and approximations. We show the divergent behaviour of approximations used in GCMs for Earth, but applied for non Earth-like planets.
△ Less
Submitted 3 July, 2023;
originally announced July 2023.
-
A Mini-Chemical Scheme with Net Reactions for 3D GCMs II. 3D thermochemical modelling of WASP-39b and HD 189733b
Authors:
Elspeth K. H. Lee,
Shang-Min Tsai,
Mark Hammond,
Xianyu Tan
Abstract:
The chemical inventory of hot Jupiter (HJ) exoplanets atmospheres continue to be observed by various ground and space based instruments in increasing detail and precision. It is expected that some HJs will exhibit strong non-equilibrium chemistry characteristics in their atmospheres, which might be inferred from spectral observations. We aim to model the three dimensional thermochemical non-equili…
▽ More
The chemical inventory of hot Jupiter (HJ) exoplanets atmospheres continue to be observed by various ground and space based instruments in increasing detail and precision. It is expected that some HJs will exhibit strong non-equilibrium chemistry characteristics in their atmospheres, which might be inferred from spectral observations. We aim to model the three dimensional thermochemical non-equilibrium chemistry in the atmospheres of the HJs WASP-39b and HD 189733b. We couple a lightweight, reduced chemical network `mini-chem' that utilises net reaction rate tables to the Exo-FMS General Circulation Model (GCM). We perform GCM models of the exoplanets WASP-39b and HD 189733b as case studies of the coupled mini-chem scheme. The GCM results are then post-processed using the 3D radiative-transfer model gCMCRT to produce transmission and emission spectra to assess the impact of non-equilibrium chemistry on their observable properties. Both simulations show significant departures from chemical equilibrium (CE) due to the dynamical motions of the atmosphere. The spacial distribution of species generally follows closely the dynamical features of the atmosphere rather than the temperature field. Each molecular species exhibits a different quench level in the simulations, also dependent on the latitude of the planet. Major differences are seen in the transmission and emission spectral features between the CE and kinetic models. Our simulations indicate that considering the 3D kinetic chemical structures of HJ atmospheres has an important impact on physical interpretation of observational data. Drawing bulk atmospheric parameters from fitting feature strengths may lead to inaccurate interpretation of chemical conditions in the atmosphere of HJs. Our open source mini-chem module is simple to couple with contemporary HJ GCM models without substantially increasing required computational resources.
△ Less
Submitted 19 February, 2023;
originally announced February 2023.
-
Proximity-induced quasi-one-dimensional superconducting quantum anomalous Hall state: a promising scalable top-down approach towards localized Majorana modes
Authors:
Omargeldi Atanov,
Wai Ting Tai,
Ying-Ming Xie,
Yat Hei Ng,
Molly A. Hammond,
Tin Seng Manfred Ho,
Tsin Hei Koo,
Hui Li,
Sui Lun Ho,
Jian Lyu,
Sukong Chong,
Peng Zhang,
Lixuan Tai,
Jiannong Wang,
Kam Tuen Law,
Kang L. Wang,
Rolf Lortz
Abstract:
In this work, ~100 nm wide quantum anomalous Hall insulator (QAHI) nanoribbons are etched from a two-dimensional QAHI film. One part of the nanoribbon is covered with superconducting Nb, while the other part is connected to an Au lead via two-dimensional QAHI regions. Andreev reflection spectroscopy measurements were performed, and multiple in-gap conductance peaks were observed in three different…
▽ More
In this work, ~100 nm wide quantum anomalous Hall insulator (QAHI) nanoribbons are etched from a two-dimensional QAHI film. One part of the nanoribbon is covered with superconducting Nb, while the other part is connected to an Au lead via two-dimensional QAHI regions. Andreev reflection spectroscopy measurements were performed, and multiple in-gap conductance peaks were observed in three different devices. In the presence of an increasing magnetic field perpendicular to the QAHI film, the multiple in-gap peak structure evolves into a single zero-bias conductance peak (ZBCP). Theoretical simulations suggest that the measurements are consistent with the scenario that the increasing magnetic field drives the nanoribbons from a multi-channel occupied regime to a single channel occupied regime, and that the ZBCP may be induced by zero energy Majorana modes as previously predicted [24]. Although further experiments are needed to clarify the nature of the ZBCP, we provide initial evidence that quasi-1D QAHI nanoribbon/superconductor heterostructures are new and promising platforms for realizing zero-energy Majorana modes.
△ Less
Submitted 13 February, 2023;
originally announced February 2023.
-
A broadband thermal emission spectrum of the ultra-hot Jupiter WASP-18b
Authors:
Louis-Philippe Coulombe,
Björn Benneke,
Ryan Challener,
Anjali A. A. Piette,
Lindsey S. Wiser,
Megan Mansfield,
Ryan J. MacDonald,
Hayley Beltz,
Adina D. Feinstein,
Michael Radica,
Arjun B. Savel,
Leonardo A. Dos Santos,
Jacob L. Bean,
Vivien Parmentier,
Ian Wong,
Emily Rauscher,
Thaddeus D. Komacek,
Eliza M. -R. Kempton,
Xianyu Tan,
Mark Hammond,
Neil T. Lewis,
Michael R. Line,
Elspeth K. H. Lee,
Hinna Shivkumar,
Ian J. M. Crossfield
, et al. (51 additional authors not shown)
Abstract:
Close-in giant exoplanets with temperatures greater than 2,000 K (''ultra-hot Jupiters'') have been the subject of extensive efforts to determine their atmospheric properties using thermal emission measurements from the Hubble and Spitzer Space Telescopes. However, previous studies have yielded inconsistent results because the small sizes of the spectral features and the limited information conten…
▽ More
Close-in giant exoplanets with temperatures greater than 2,000 K (''ultra-hot Jupiters'') have been the subject of extensive efforts to determine their atmospheric properties using thermal emission measurements from the Hubble and Spitzer Space Telescopes. However, previous studies have yielded inconsistent results because the small sizes of the spectral features and the limited information content of the data resulted in high sensitivity to the varying assumptions made in the treatment of instrument systematics and the atmospheric retrieval analysis. Here we present a dayside thermal emission spectrum of the ultra-hot Jupiter WASP-18b obtained with the NIRISS instrument on JWST. The data span 0.85 to 2.85 $μ$m in wavelength at an average resolving power of 400 and exhibit minimal systematics. The spectrum shows three water emission features (at $>$6$σ$ confidence) and evidence for optical opacity, possibly due to H$^-$, TiO, and VO (combined significance of 3.8$σ$). Models that fit the data require a thermal inversion, molecular dissociation as predicted by chemical equilibrium, a solar heavy element abundance (''metallicity'', M/H = 1.03$_{-0.51}^{+1.11}$ $\times$ solar), and a carbon-to-oxygen (C/O) ratio less than unity. The data also yield a dayside brightness temperature map, which shows a peak in temperature near the sub-stellar point that decreases steeply and symmetrically with longitude toward the terminators.
△ Less
Submitted 20 January, 2023; v1 submitted 19 January, 2023;
originally announced January 2023.
-
Reduced Order Model of a Generic Submarine for Maneuvering Near the Surface
Authors:
J. Ezequiel Martin,
Maxwell Hammond,
Nicholas Rober,
Yakin Kim,
Venanzio Cichella,
Pablo Carrica
Abstract:
A reduced order model of a generic submarine is presented. Computational fluid dynamics (CFD) results are used to create and validate a model that includes depth dependence and the effect of waves on the craft. The model and the procedure to obtain its coefficients are discussed, and examples of the data used to obtain the model coefficients are presented. An example of operation following a compl…
▽ More
A reduced order model of a generic submarine is presented. Computational fluid dynamics (CFD) results are used to create and validate a model that includes depth dependence and the effect of waves on the craft. The model and the procedure to obtain its coefficients are discussed, and examples of the data used to obtain the model coefficients are presented. An example of operation following a complex path is presented and results from the reduced order model are compared to those from an equivalent CFD calculation. The controller implemented to complete these maneuvers is also presented.
△ Less
Submitted 19 December, 2022;
originally announced December 2022.
-
Inverse-designed lithium niobate nanophotonics
Authors:
Chengfei Shang,
Jingwei Yang,
Alec M. Hammond,
Zhaoxi Chen,
Mo Chen,
Zin Lin,
Steven G. Johnson,
Cheng Wang
Abstract:
Lithium niobate-on-insulator (LNOI) is an emerging photonic platform that exhibits favorable material properties (such as low optical loss, strong nonlinearities, and stability) and enables large-scale integration with stronger optical confinement, showing promise for future optical networks, quantum processors, and nonlinear optical systems. However, while photonics engineering has entered the er…
▽ More
Lithium niobate-on-insulator (LNOI) is an emerging photonic platform that exhibits favorable material properties (such as low optical loss, strong nonlinearities, and stability) and enables large-scale integration with stronger optical confinement, showing promise for future optical networks, quantum processors, and nonlinear optical systems. However, while photonics engineering has entered the era of automated "inverse design" via optimization in recent years, the design of LNOI integrated photonic devices still mostly relies on intuitive models and inefficient parameter sweeps, limiting the accessible parameter space, performance, and functionality. Here, we develop and implement a 3D gradient-based inverse-design model tailored for topology optimization of the LNOI platform, which not only could efficiently search a large parameter space but also takes into account practical fabrication constraints, including minimum feature sizes and etched sidewall angles. We experimentally demonstrate a spatial-mode multiplexer, a waveguide crossing, and a compact waveguide bend, all with low insertion losses, tiny footprints, and excellent agreement between simulation and experimental results. The devices, together with the design methodology, represent a crucial step towards the variety of advanced device functionalities needed in future LNOI photonics, and could provide compact and cost-effective solutions for future optical links, quantum technologies and nonlinear optics.
△ Less
Submitted 19 December, 2022;
originally announced December 2022.
-
CAMEMBERT: A Mini-Neptunes GCM Intercomparison, Protocol Version 1.0. A CUISINES Model Intercomparison Project
Authors:
Duncan A. Christie,
Elspeth K. H. Lee,
Hamish Innes,
Pascal A. Noti,
Benjamin Charnay,
Thomas J. Fauchez,
Nathan J. Mayne,
Russell Deitrick,
Feng Ding,
Jennifer J. Greco,
Mark Hammond,
Isaac Malsky,
Avi Mandell,
Emily Rauscher,
Michael T. Roman,
Denis E. Sergeev,
Linda Sohl,
Maria E. Steinrueck,
Martin Turbet,
Eric T. Wolf,
Maria Zamyatina,
Ludmila Carone
Abstract:
With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the…
▽ More
With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the exoplanetary science community to simulate the atmospheres of mini-Neptunes. We focus on two targets well studied both observationally and theoretically with planned JWST Cycle 1 observations: the warm GJ~1214b and the cooler K2-18b. For each target, we consider a temperature-forced case, a clear sky dual-grey radiative transfer case, and a clear sky multi band radiative transfer case, covering a range of complexities and configurations where we know differences exist between GCMs in the literature. This paper presents all the details necessary to participate in the intercomparison, with the intention of presenting the results in future papers. Currently, there are eight GCMs participating (ExoCAM, Exo-FMS, FMS PCM, Generic PCM, MITgcm, RM-GCM, THOR, and the UM), and membership in the project remains open. Those interested in participating are invited to contact the authors.
△ Less
Submitted 8 November, 2022;
originally announced November 2022.
-
Temperature Structures Associated with Different Components of the Atmospheric Circulation on Tidally Locked Exoplanets
Authors:
Neil T. Lewis,
Mark Hammond
Abstract:
Observations of time-resolved thermal emission from tidally locked exoplanets can tell us about their atmospheric temperature structure. Telescopes such as JWST and ARIEL will improve the quality and availability of these measurements. This motivates an improved understanding of the processes that determine atmospheric temperature structure, particularly atmospheric circulation. The circulation is…
▽ More
Observations of time-resolved thermal emission from tidally locked exoplanets can tell us about their atmospheric temperature structure. Telescopes such as JWST and ARIEL will improve the quality and availability of these measurements. This motivates an improved understanding of the processes that determine atmospheric temperature structure, particularly atmospheric circulation. The circulation is important in determining atmospheric temperatures, not only through its ability to transport heat, but also because any circulation pattern needs to be balanced by horizontal pressure contrasts, therefore implying a particular temperature structure. In this work, we show how the global temperature field on a tidally locked planet can be decomposed into contributions that are balanced by different components of the atmospheric circulation. These are the superrotating jet, stationary Rossby waves, and the divergent circulation. To achieve this, we partition the geopotential field into components balanced by the divergent circulation and the rotational circulation, with the latter comprising the jet and Rossby waves. The partitioned geopotential then implies a corresponding partitioning of the temperature via the hydrostatic relation. We apply these diagnostics to idealised general circulation model simulations, to show how the separate rotational and divergent circulations together make up the total three-dimensional atmospheric temperature structure. We also show how each component contributes distinct signatures to the thermal phase curve of a tidally locked planet. We conclude that this decomposition is a physically meaningful separation of the temperature field that explains its global structure, and can be used to fit observations of thermal emission.
△ Less
Submitted 13 April, 2022;
originally announced April 2022.
-
3D Path Following and L1 Adaptive Control for Underwater Vehicles
Authors:
Nicholas Rober,
Maxwell Hammond,
Venanzio Cichella,
Juan E. Martin,
Pablo Carrica
Abstract:
This paper addresses the problem of guidance and control of underwater vehicles. A multi-level control strategy is used to determine (1) outer-loop path-following commands and (2) inner-loop actuation commands. Specifically, a line-of-sight path-following algorithm is used to guide the vehicle along a three-dimensional path, and an $\mathcal{L}_1$ adaptive control algorithm is used to determine th…
▽ More
This paper addresses the problem of guidance and control of underwater vehicles. A multi-level control strategy is used to determine (1) outer-loop path-following commands and (2) inner-loop actuation commands. Specifically, a line-of-sight path-following algorithm is used to guide the vehicle along a three-dimensional path, and an $\mathcal{L}_1$ adaptive control algorithm is used to determine the low-level rudder commands to accomplish path following. The performance bounds of these outer- and inner-loop control algorithms are presented. Numerical results obtained using a physics-based Simulink model are used to aid in visualization of the control algorithm's performance.
△ Less
Submitted 2 April, 2022;
originally announced April 2022.
-
K2 and Spitzer phase curves of the rocky ultra-short-period planet K2-141 b hint at a tenuous rock vapor atmosphere
Authors:
S. Zieba,
M. Zilinskas,
L. Kreidberg,
T. G. Nguyen,
Y. Miguel,
N. B. Cowan,
R. Pierrehumbert,
L. Carone,
L. Dang,
M. Hammond,
T. Louden,
R. Lupu,
L. Malavolta,
K. B. Stevenson
Abstract:
K2-141 b is a transiting, small (1.5 Re) ultra-short-period (USP) planet discovered by Kepler orbiting a K-dwarf host star every 6.7 hours. The planet's high surface temperature makes it an excellent target for thermal emission observations. Here we present 65 hours of continuous photometric observations of K2-141 b collected with Spitzer's IRAC Channel 2 at 4.5 micron spanning 10 full orbits of t…
▽ More
K2-141 b is a transiting, small (1.5 Re) ultra-short-period (USP) planet discovered by Kepler orbiting a K-dwarf host star every 6.7 hours. The planet's high surface temperature makes it an excellent target for thermal emission observations. Here we present 65 hours of continuous photometric observations of K2-141 b collected with Spitzer's IRAC Channel 2 at 4.5 micron spanning 10 full orbits of the planet. We measure an infrared eclipse depth of 143 +/- 39 ppm and a peak to trough amplitude variation of 121 +/- 43 ppm. The best fit model to the Spitzer data shows no significant thermal hotspot offset, in contrast to the previously observed offset for the well-studied USP planet 55 Cnc e. We also jointly analyze the new Spitzer observations with the photometry collected by Kepler during two separate K2 campaigns. We model the planetary emission with a range of toy models that include a reflective and a thermal contribution. With a two-temperature model, we measure a dayside temperature of 2049 +/- 361 K and a night-side temperature that is consistent with zero (Tp,n < 1712 K at 2 sigma). Models with a steep dayside temperature gradient provide a better fit to the data than a uniform dayside temperature (DeltaBIC = 22.2). We also find evidence for a non-zero geometric albedo Ag = 0.28 +/- 0.07. We also compare the data to a physically motivated, pseudo-2D rock vapor model and a 1D turbulent boundary layer model. Both models fit the data well. Notably, we find that the optical eclipse depth can be explained by thermal emission from a hot inversion layer, rather than reflected light. A thermal inversion may also be responsible for the deep optical eclipse observed for another USP, Kepler-10 b. Finally, we significantly improve the ephemerides for K2-141 b and c, which will facilitate further follow-up observations of this interesting system with state-of-the-art observatories like JWST.
△ Less
Submitted 1 March, 2022;
originally announced March 2022.
-
Automatic Correction of Syntactic Dependency Annotation Differences
Authors:
Andrew Zupon,
Andrew Carnie,
Michael Hammond,
Mihai Surdeanu
Abstract:
Annotation inconsistencies between data sets can cause problems for low-resource NLP, where noisy or inconsistent data cannot be as easily replaced compared with resource-rich languages. In this paper, we propose a method for automatically detecting annotation mismatches between dependency parsing corpora, as well as three related methods for automatically converting the mismatches. All three meth…
▽ More
Annotation inconsistencies between data sets can cause problems for low-resource NLP, where noisy or inconsistent data cannot be as easily replaced compared with resource-rich languages. In this paper, we propose a method for automatically detecting annotation mismatches between dependency parsing corpora, as well as three related methods for automatically converting the mismatches. All three methods rely on comparing an unseen example in a new corpus with similar examples in an existing corpus. These three methods include a simple lexical replacement using the most frequent tag of the example in the existing corpus, a GloVe embedding-based replacement that considers a wider pool of examples, and a BERT embedding-based replacement that uses contextualized embeddings to provide examples fine-tuned to our specific data. We then evaluate these conversions by retraining two dependency parsers -- Stanza (Qi et al. 2020) and Parsing as Tagging (PaT) (Vacareanu et al. 2020) -- on the converted and unconverted data. We find that applying our conversions yields significantly better performance in many cases. Some differences observed between the two parsers are observed. Stanza has a more complex architecture with a quadratic algorithm, so it takes longer to train, but it can generalize better with less data. The PaT parser has a simpler architecture with a linear algorithm, speeding up training time but requiring more training data to reach comparable or better performance.
△ Less
Submitted 15 January, 2022;
originally announced January 2022.
-
General Circulation Model Errors are Variable across Exoclimate Parameter Spaces
Authors:
Pushkar Kopparla,
Russell Deitrick,
Kevin Heng,
João M. Mendonça,
Mark Hammond
Abstract:
General circulation models are often used to explore exoclimate parameter spaces and classify atmospheric circulation regimes. Models are tuned to give reasonable climate states for standard test cases, such as the Held-Suarez test, and then used to simulate diverse exoclimates by varying input parameters such as rotation rates, instellation, atmospheric optical properties, frictional timescales a…
▽ More
General circulation models are often used to explore exoclimate parameter spaces and classify atmospheric circulation regimes. Models are tuned to give reasonable climate states for standard test cases, such as the Held-Suarez test, and then used to simulate diverse exoclimates by varying input parameters such as rotation rates, instellation, atmospheric optical properties, frictional timescales and so on. In such studies, there is an implicit assumption that the model which works reasonably well for the standard test case will be credible at all points in an arbitrarily wide parameter space. Here, we test this assumption using the open-source general circulation model THOR to simulate atmospheric circulation on tidally locked Earth-like planets with rotation periods of 0.1 to 100 days. We find that the model error, as quantified by the ratio between physical and spurious numerical contributions to the angular momentum balance, is extremely variable across this range of rotation periods with some cases where numerical errors are the dominant component. Increasing model grid resolution does improve errors but using a higher-order numerical diffusion scheme can sometimes magnify errors for finite-volume dynamical solvers. We further show that to minimize error and make the angular momentum balance more physical within our model, the surface friction timescale must be smaller than the rotational timescale.
△ Less
Submitted 21 October, 2021;
originally announced October 2021.
-
Distributed Control of Truss Robots Using Consensus Alternating Direction Method of Multipliers
Authors:
Nathan S. Usevitch,
Trevor Halsted,
Zachary M. Hammond,
Allison M. Okamura,
Mac Schwager
Abstract:
Truss robots, or robots that consist of extensible links connected at universal joints, are often designed with modular physical components but require centralized control techniques. This paper presents a distributed control technique for truss robots. The truss robot is viewed as a collective, where each individual node of the robot is capable of measuring the lengths of the neighboring edges, c…
▽ More
Truss robots, or robots that consist of extensible links connected at universal joints, are often designed with modular physical components but require centralized control techniques. This paper presents a distributed control technique for truss robots. The truss robot is viewed as a collective, where each individual node of the robot is capable of measuring the lengths of the neighboring edges, communicating with a subset of the other nodes, and computing and executing its own control actions with its connected edges. Through an iterative distributed optimization, the individual members utilize local information to converge on a global estimate of the robot's state, and then coordinate their planned motion to achieve desired global behavior. This distributed optimization is based on a consensus alternating direction method of multipliers framework. This distributed algorithm is then adapted to control an isoperimetric truss robot, and the distributed algorithm is used in an experimental demonstration. The demonstration allows a user to broadcast commands to a single node of the robot, which then ensures the coordinated motion of all other nodes to achieve the desired global motion.
△ Less
Submitted 14 August, 2021;
originally announced August 2021.
-
Simulating gas giant exoplanet atmospheres with Exo-FMS: Comparing semi-grey, picket fence and correlated-k radiative-transfer schemes
Authors:
Elspeth K. H. Lee,
Vivien Parmentier,
Mark Hammond,
Simon L. Grimm,
Daniel Kitzmann,
Xianyu Tan,
Shang-Min Tsai,
Raymond T. Pierrehumbert
Abstract:
Radiative-transfer (RT) is a fundamental part of modelling exoplanet atmospheres with general circulation models (GCMs). An accurate RT scheme is required for estimates of the atmospheric energy transport and for gaining physical insight from model spectra. We implement three RT schemes for Exo-FMS: semi-grey, non-grey `picket fence', and real gas with correlated-k. We benchmark the Exo-FMS GCM us…
▽ More
Radiative-transfer (RT) is a fundamental part of modelling exoplanet atmospheres with general circulation models (GCMs). An accurate RT scheme is required for estimates of the atmospheric energy transport and for gaining physical insight from model spectra. We implement three RT schemes for Exo-FMS: semi-grey, non-grey `picket fence', and real gas with correlated-k. We benchmark the Exo-FMS GCM using these RT schemes to hot Jupiter simulation results from the literature. We perform a HD 209458b-like simulation with the three schemes and compare their results. These simulations are then post-processed to compare their observable differences. The semi-grey scheme results show qualitative agreement with previous studies in line with variations seen between GCM models. The real gas model reproduces well the temperature and dynamical structures from other studies. After post-processing our non-grey picket fence scheme compares very favourably with the real gas model, producing similar transmission spectra, emission spectra and phase curve behaviours. Exo-FMS is able to reliably reproduce the essential features of contemporary GCM models in the hot gas giant regime. Our results suggest the picket fence approach offers a simple way to improve upon RT realism beyond semi-grey schemes.
△ Less
Submitted 28 June, 2021; v1 submitted 22 June, 2021;
originally announced June 2021.
-
The rotational and divergent components of atmospheric circulation on tidally locked planets
Authors:
Mark Hammond,
Neil T. Lewis
Abstract:
Tidally locked exoplanets likely host global atmospheric circulations with a superrotating equatorial jet, planetary-scale stationary waves and thermally-driven overturning circulation. In this work, we show that each of these features can be separated from the total circulation by using a Helmholtz decomposition, which splits the circulation into rotational (divergence free) and divergent (vortic…
▽ More
Tidally locked exoplanets likely host global atmospheric circulations with a superrotating equatorial jet, planetary-scale stationary waves and thermally-driven overturning circulation. In this work, we show that each of these features can be separated from the total circulation by using a Helmholtz decomposition, which splits the circulation into rotational (divergence free) and divergent (vorticity free) components. This technique is applied to the simulated circulation of a terrestrial planet and a gaseous hot Jupiter. For both planets, the rotational component comprises the equatorial jet and stationary waves, and the divergent component contains the overturning circulation. Separating out each component allows us to evaluate their spatial structure and relative contribution to the total flow. In contrast with previous work, we show that divergent velocities are not negligible when compared with rotational velocities, and that divergent, overturning circulation takes the form of a single, roughly isotropic cell that ascends on the day-side and descends on the night-side. These conclusions are drawn for both the terrestrial case and the hot Jupiter. To illustrate the utility of the Helmholtz decomposition for studying atmospheric processes, we compute the contribution of each of the circulation components to heat transport from day- to night-side. Surprisingly, we find that the divergent circulation dominates day-night heat transport in the terrestrial case and accounts for around half of the heat transport for the hot Jupiter. The relative contributions of the rotational and divergent components to day-night heat transport are likely sensitive to multiple planetary parameters and atmospheric processes, and merit further study.
△ Less
Submitted 23 February, 2021;
originally announced February 2021.
-
Vertically resolved magma ocean-protoatmosphere evolution: H$_2$, H$_2$O, CO$_2$, CH$_4$, CO, O$_2$, and N$_2$ as primary absorbers
Authors:
Tim Lichtenberg,
Dan J. Bower,
Mark Hammond,
Ryan Boukrouche,
Patrick Sanan,
Shang-Min Tsai,
Raymond T. Pierrehumbert
Abstract:
The earliest atmospheres of rocky planets originate from extensive volatile release during magma ocean epochs that occur during assembly of the planet. These establish the initial distribution of the major volatile elements between different chemical reservoirs that subsequently evolve via geological cycles. Current theoretical techniques are limited in exploring the anticipated range of compositi…
▽ More
The earliest atmospheres of rocky planets originate from extensive volatile release during magma ocean epochs that occur during assembly of the planet. These establish the initial distribution of the major volatile elements between different chemical reservoirs that subsequently evolve via geological cycles. Current theoretical techniques are limited in exploring the anticipated range of compositional and thermal scenarios of early planetary evolution, even though these are of prime importance to aid astronomical inferences on the environmental context and geological history of extrasolar planets. Here, we present a coupled numerical framework that links an evolutionary, vertically-resolved model of the planetary silicate mantle with a radiative-convective model of the atmosphere. Using this method we investigate the early evolution of idealized Earth-sized rocky planets with end-member, clear-sky atmospheres dominated by either H$_2$, H$_2$O, CO$_2$, CH$_4$, CO, O$_2$, or N$_2$. We find central metrics of early planetary evolution, such as energy gradient, sequence of mantle solidification, surface pressure, or vertical stratification of the atmosphere, to be intimately controlled by the dominant volatile and outgassing history of the planet. Thermal sequences fall into three general classes with increasing cooling timescale: CO, N$_2$, and O$_2$ with minimal effect, H$_2$O, CO$_2$, and CH$_4$ with intermediate influence, and H$_2$ with several orders of magnitude increase in solidification time and atmosphere vertical stratification. Our numerical experiments exemplify the capabilities of the presented modeling framework and link the interior and atmospheric evolution of rocky exoplanets with multi-wavelength astronomical observations.
△ Less
Submitted 26 January, 2021;
originally announced January 2021.
-
The Equatorial Jet Speed on Tidally Locked Planets: I -- Terrestrial Planets
Authors:
Mark Hammond,
Shang-Min Tsai,
Raymond T. Pierrehumbert
Abstract:
The atmospheric circulation of tidally locked planets is dominated by a superrotating eastward equatorial jet. We develop a predictive theory for the formation of this jet, proposing a mechanism in which the three-dimensional stationary waves induced by the day-night forcing gradient produce an equatorial acceleration. This is balanced in equilibrium by an interaction between the resulting jet and…
▽ More
The atmospheric circulation of tidally locked planets is dominated by a superrotating eastward equatorial jet. We develop a predictive theory for the formation of this jet, proposing a mechanism in which the three-dimensional stationary waves induced by the day-night forcing gradient produce an equatorial acceleration. This is balanced in equilibrium by an interaction between the resulting jet and the vertical motion of the atmosphere. The three-dimensional structure of the zonal acceleration is vital to this mechanism.
We demonstrate this mechanism in a hierarchy of models. We calculate the three-dimensional stationary waves induced by the forcing on these planets, and show the vertical structure of the zonal acceleration produced by these waves, which we use to suggest a mechanism for how the jet forms. GCM simulations are used to confirm the equilibrium state predicted by this mechanism, where the acceleration from these waves is balanced by an interaction between the zonal-mean vertical velocity and the jet. We derive a simple model of this using the "Weak Temperature Gradient" approximation, which gives an estimate of the jet speed on a terrestrial tidally locked planet.
We conclude that the proposed mechanism is a good description of the formation of an equatorial jet on a terrestrial tidally locked planet, and should be useful for interpreting observations and simulations of these planets. The mechanism requires assumptions such as a large equatorial Rossby radius and weak acceleration due to transient waves, and a different mechanism may produce the equatorial jets on gaseous tidally locked planets.
△ Less
Submitted 1 September, 2020;
originally announced September 2020.
-
Factorized Machine Learning for Performance Modeling of Massively Parallel Heterogeneous Physical Simulations
Authors:
Ardavan Oskooi,
Christopher Hogan,
Alec M. Hammond,
M. T. Homer Reid,
Steven G. Johnson
Abstract:
We demonstrate neural-network runtime prediction for complex, many-parameter, massively parallel, heterogeneous-physics simulations running on cloud-based MPI clusters. Because individual simulations are so expensive, it is crucial to train the network on a limited dataset despite the potentially large input space of the physics at each point in the spatial domain. We achieve this using a two-part…
▽ More
We demonstrate neural-network runtime prediction for complex, many-parameter, massively parallel, heterogeneous-physics simulations running on cloud-based MPI clusters. Because individual simulations are so expensive, it is crucial to train the network on a limited dataset despite the potentially large input space of the physics at each point in the spatial domain. We achieve this using a two-part strategy. First, we perform data-driven static load balancing using regression coefficients extracted from small simulations, which both improves parallel performance and reduces the dependency of the runtime on the precise spatial layout of the heterogeneous physics. Second, we divide the execution time of these load-balanced simulations into computation and communication, factoring crude asymptotic scalings out of each term, and training neural nets for the remaining factor coefficients. This strategy is implemented for Meep, a popular and complex open-source electrodynamics simulation package, and are validated for heterogeneous simulations drawn from published engineering models.
△ Less
Submitted 6 October, 2020; v1 submitted 9 March, 2020;
originally announced March 2020.
-
Simplified 3D GCM modelling of the irradiated brown dwarf WD0137-349B
Authors:
Elspeth K. H. Lee,
Sarah L. Casewell,
Katy L. Chubb,
Mark Hammond,
Xianyu Tan,
Shang-Min Tsai,
Raymond T. Pierrehumbert
Abstract:
Context: White dwarf - Brown dwarf short period binaries (P$_{\rm orb}$ $\lesssim$ 2 hours) are some of the most extreme irradiated atmospheric environments known. These systems offer an opportunity to explore theoretical and modelling efforts of irradiated atmospheres different to typical hot Jupiter systems. Aims: We aim to investigate the three dimensional atmospheric structural and dynamical…
▽ More
Context: White dwarf - Brown dwarf short period binaries (P$_{\rm orb}$ $\lesssim$ 2 hours) are some of the most extreme irradiated atmospheric environments known. These systems offer an opportunity to explore theoretical and modelling efforts of irradiated atmospheres different to typical hot Jupiter systems. Aims: We aim to investigate the three dimensional atmospheric structural and dynamical properties of the Brown dwarf WD0137-349B. Methods: We use the three dimensional GCM model Exo-FMS, with a dual-band grey radiative-transfer scheme to model the atmosphere of WD0137-349B. The results of the GCM model are post-processed using the three dimensional Monte Carlo radiative-transfer model \textsc{cmcrt}. Results: Our results suggest inefficient day-night energy transport and a large day-night temperature contrast for WD0137-349B. Multiple flow patterns are present, shifting energy asymmetrically eastward or westward depending on their zonal direction and latitude. Regions of overturning are produced on the western terminator. We are able to reproduce the start of the system near-IR emission excess at $\gtrsim$ 1.95 $μ$m as observed by the GNIRS instrument. Our model over predicts the IR phase curve fluxes by factors of $\approx$1-3, but generally fits the shape of the phase curves well. Chemical kinetic modelling using \textsc{vulcan} suggests a highly ionised region at high altitudes can form on the dayside of the Brown dwarf. Conclusions: We present a first attempt at simulating the atmosphere of a short period White dwarf - Brown dwarf binary in a 3D setting. Further studies into the radiative and photochemical heating from the UV irradiation is required to more accurately capture the energy balance inside the Brown dwarf atmosphere. Cloud formation may also play an important role in shaping the emission spectra of the Brown dwarf.
△ Less
Submitted 25 June, 2020; v1 submitted 17 January, 2020;
originally announced January 2020.
-
Accelerating Silicon Photonic Parameter Extraction using Artificial Neural Networks
Authors:
Alec M. Hammond,
Easton Potokar,
Ryan M. Camacho
Abstract:
We present a novel silicon photonic parameter extraction tool that uses artificial neural networks. While other parameter extraction methods are restricted to relatively simple devices whose responses are easily modeled by analytic transfer functions, this method is capable of extracting parameters for any device with a discrete number of design parameters. To validate the method, we design and fa…
▽ More
We present a novel silicon photonic parameter extraction tool that uses artificial neural networks. While other parameter extraction methods are restricted to relatively simple devices whose responses are easily modeled by analytic transfer functions, this method is capable of extracting parameters for any device with a discrete number of design parameters. To validate the method, we design and fabricate integrated chirped Bragg gratings. We then estimate the actual device parameters by iteratively fitting the simultaneously measured group delay and reflection profiles to the artificial neural network output. The method is fast, accurate, and capable of modeling the complicated chirping and index contrast.
△ Less
Submitted 23 January, 2019;
originally announced January 2019.
-
Designing Silicon Photonic Devices using Artificial Neural Networks
Authors:
Alec M. Hammond,
Ryan M. Camacho
Abstract:
We develop and experimentally validate a novel neural network design framework for silicon photonics devices that is both practical and intuitive. The framework is applicable to nearly all known integrated photonics devices, but as case studies we consider simple waveguides and chirped Bragg Gratings. By using artificial neural networks, we decrease the computational cost relative to traditional d…
▽ More
We develop and experimentally validate a novel neural network design framework for silicon photonics devices that is both practical and intuitive. The framework is applicable to nearly all known integrated photonics devices, but as case studies we consider simple waveguides and chirped Bragg Gratings. By using artificial neural networks, we decrease the computational cost relative to traditional design methodologies by more than 4 orders of magnitude. We also demonstrate the abstraction of the device models to a few simple input and output parameters relevant to designers. We then apply the results to various design problems and experimentally compare fabricated devices to the neural network's predictions.
△ Less
Submitted 10 January, 2019; v1 submitted 1 November, 2018;
originally announced December 2018.
-
Wave-Mean Flow Interactions in the Atmospheric Circulation of Tidally Locked Planets
Authors:
Mark Hammond,
Raymond T. Pierrehumbert
Abstract:
We use a linear shallow-water model to investigate the global circulation of the atmospheres of tidally locked planets. Simulations, observations, and simple models show that if these planets are sufficiently rapidly rotating, their atmospheres have an eastward equatorial jet and a hot-spot east of the substellar point. We linearize the shallow-water model about this eastward flow and its associat…
▽ More
We use a linear shallow-water model to investigate the global circulation of the atmospheres of tidally locked planets. Simulations, observations, and simple models show that if these planets are sufficiently rapidly rotating, their atmospheres have an eastward equatorial jet and a hot-spot east of the substellar point. We linearize the shallow-water model about this eastward flow and its associated height perturbation. The forced solutions of this system show that the shear flow explains the form of the global circulation, particularly the hot-spot shift and the positions of the cold standing waves on the night-side.
We suggest that the eastward hot-spot shift seen in observations and 3D simulations of these atmospheres is caused by the zonal flow Doppler-shifting the stationary wave response eastwards, summed with the height perturbation from the flow itself. This differs from other studies which explained the hot-spot shift as pure advection of heat from air flowing eastward from the substellar point, or as equatorial waves travelling eastwards. We compare our solutions to simulations in our climate model Exo-FMS, and show that the height fields and wind patterns match.
We discuss how planetary properties affect the global circulation, and how they change observables such as the hot-spot shift or day-night contrast. We conclude that the wave-mean flow interaction between the stationary planetary waves and the equatorial jet is a vital part of the equilibrium circulation on tidally locked planets.
△ Less
Submitted 26 October, 2018;
originally announced October 2018.
-
Error Correction in Structured Optical Receivers
Authors:
Alec M. Hammond,
Ian W. Frank,
Ryan M. Camacho
Abstract:
Integrated optics Green Machines enable better communication in photon-starved environments, but fabrication inconsistencies induce unpredictable internal phase errors, making them difficult to construct. We describe and experimentally demonstrate a new method to compensate for arbitrary phase errors by deriving a convex error space and implementing an algorithm to learn a unique codebook of codew…
▽ More
Integrated optics Green Machines enable better communication in photon-starved environments, but fabrication inconsistencies induce unpredictable internal phase errors, making them difficult to construct. We describe and experimentally demonstrate a new method to compensate for arbitrary phase errors by deriving a convex error space and implementing an algorithm to learn a unique codebook of codewords corresponding to each matrix.
△ Less
Submitted 27 February, 2018;
originally announced February 2018.
-
Linking the Climate and Thermal Phase Curve of 55 Cancri e
Authors:
Mark Hammond,
Raymond Pierrehumbert
Abstract:
The thermal phase curve of 55 Cancri e is the first measurement of the temperature distribution of a tidally locked Super-Earth, but raises a number of puzzling questions about the planet's climate. The phase curve has a high amplitude and peak offset, suggesting that it has a significant eastward hot-spot shift as well as a large day-night temperature contrast. We use a general circulation model…
▽ More
The thermal phase curve of 55 Cancri e is the first measurement of the temperature distribution of a tidally locked Super-Earth, but raises a number of puzzling questions about the planet's climate. The phase curve has a high amplitude and peak offset, suggesting that it has a significant eastward hot-spot shift as well as a large day-night temperature contrast. We use a general circulation model to model potential climates, and investigate the relation between bulk atmospheric composition and the magnitude of these seemingly contradictory features. We confirm theoretical models of tidally locked circulation are consistent with our numerical model of 55 Cnc e, and rule out certain atmospheric compositions based on their thermodynamic properties. Our best-fitting atmosphere has a significant hot-spot shift and day-night contrast, although these are not as large as the observed phase curve. We discuss possible physical processes which could explain the observations, and show that night-side cloud formation from species such as SiO from a day-side magma ocean could potentially increase the phase curve amplitude and explain the observations. We conclude that the observations could be explained by an optically thick atmosphere with a low mean molecular weight, a surface pressure of several bar and a strong eastward circulation, with night-side cloud formation a possible explanation for the difference between our model and the observations.
△ Less
Submitted 10 October, 2017;
originally announced October 2017.
-
Creating Causal Embeddings for Question Answering with Minimal Supervision
Authors:
Rebecca Sharp,
Mihai Surdeanu,
Peter Jansen,
Peter Clark,
Michael Hammond
Abstract:
A common model for question answering (QA) is that a good answer is one that is closely related to the question, where relatedness is often determined using general-purpose lexical models such as word embeddings. We argue that a better approach is to look for answers that are related to the question in a relevant way, according to the information need of the question, which may be determined throu…
▽ More
A common model for question answering (QA) is that a good answer is one that is closely related to the question, where relatedness is often determined using general-purpose lexical models such as word embeddings. We argue that a better approach is to look for answers that are related to the question in a relevant way, according to the information need of the question, which may be determined through task-specific embeddings. With causality as a use case, we implement this insight in three steps. First, we generate causal embeddings cost-effectively by bootstrapping cause-effect pairs extracted from free text using a small set of seed patterns. Second, we train dedicated embeddings over this data, by using task-specific contexts, i.e., the context of a cause is its effect. Finally, we extend a state-of-the-art reranking approach for QA to incorporate these causal embeddings. We evaluate the causal embedding models both directly with a casual implication task, and indirectly, in a downstream causal QA task using data from Yahoo! Answers. We show that explicitly modeling causality improves performance in both tasks. In the QA task our best model achieves 37.3% P@1, significantly outperforming a strong baseline by 7.7% (relative).
△ Less
Submitted 26 September, 2016;
originally announced September 2016.
-
Hydraulic External Pre-Isolator System for LIGO
Authors:
S. Wen,
R. Mittleman,
K. Mason,
J. Giaime,
R. Abbott,
J. Kern,
B. O'Reilly,
R. Bork,
M. Hammond,
C. Hardham,
B. Lantz,
W. Hua,
D. Coyne,
G. Traylor,
H. Overmier,
T. Evans,
J. Hanson,
O. Spjeld,
M. Macinnis,
K. Mailand,
D. Sellers,
K. Carter,
P. Sarin
Abstract:
The Hydraulic External Pre-Isolator (HEPI) is the first 6 degrees of freedom active seismic isolation system implemented at the Laser Interferometer Gravitational Wave Observatory (LIGO). Implementation was first completed at the LIGO Livingston Observatory (LLO) prior to LIGO's 5th science run, successfully cutting down the disturbance seen by LLO's suspended optics in the two most prominent seis…
▽ More
The Hydraulic External Pre-Isolator (HEPI) is the first 6 degrees of freedom active seismic isolation system implemented at the Laser Interferometer Gravitational Wave Observatory (LIGO). Implementation was first completed at the LIGO Livingston Observatory (LLO) prior to LIGO's 5th science run, successfully cutting down the disturbance seen by LLO's suspended optics in the two most prominent seismic disturbance bands, the microseism (0.1-0.3Hz) and the anthropogenic (1-3Hz) bands, by a factor of a few to tens. The improvement in seismic isolation contributed directly to LLO's much improved duty cycle of 66.7% and LIGO's triple coincident duty cycle of 53%. We report the design, control scheme, and isolation performance of HEPI at LLO in this paper. Aided with this success, funding for incorporating HEPI into the LIGO Hanford Observatory was approved and installation is currently underway.
△ Less
Submitted 17 October, 2013; v1 submitted 22 September, 2013;
originally announced September 2013.
-
A New Catalog of Faraday Rotation Measures and Redshifts for Extragalactic Radio Sources
Authors:
Alison M. Hammond,
Timothy Robishaw,
B. M. Gaensler
Abstract:
We present a catalog of Faraday rotation measures (RMs) and redshifts for 4003 extragalactic radio sources detected at 1.4 GHz, derived by identifying optical counterparts and spectroscopic redshifts for linearly polarized radio sources from the NRAO VLA Sky Survey. This catalog is more than an order of magnitude larger than any previous sample of RM vs. redshift, and covers the redshift range 0 <…
▽ More
We present a catalog of Faraday rotation measures (RMs) and redshifts for 4003 extragalactic radio sources detected at 1.4 GHz, derived by identifying optical counterparts and spectroscopic redshifts for linearly polarized radio sources from the NRAO VLA Sky Survey. This catalog is more than an order of magnitude larger than any previous sample of RM vs. redshift, and covers the redshift range 0 < z < 5.3 ; the median redshift of the catalog is z = 0.70, and there are more than 1500 sources at redshifts z > 1. For 3650 of these sources at Galactic latitudes |b| >= 20 degrees, we present a second catalog in which we have corrected for the foreground Faraday rotation of the Milky Way, resulting in an estimate of the residual rotation measure (RRM) that aims to isolate the contribution from extragalactic magnetic fields. We find no significant evolution of RRM with redshift, but observe a strong anti-correlation between RRM and fractional polarization, p, that we argue is the result of beam depolarization from small-scale fluctuations in the foreground magnetic field or electron density. We suggest that the observed variance in RRM and the anti-correlation of RRM with p both require a population of magnetized intervening objects that lie outside the Milky Way but in the foreground to the emitting sources.
△ Less
Submitted 8 May, 2013; v1 submitted 6 September, 2012;
originally announced September 2012.
-
Parsing syllables: modeling OT computationally
Authors:
Michael Hammond
Abstract:
In this paper, I propose to implement syllabification in OT as a parser. I propose several innovations that result in a finite and small candidate set. The candidate set problem is handled with several moves: i) MAX and DEP violations are not hypothesized by the parser, ii) candidates are encoded locally, and iii) EVAL is applied constraint by constraint.
The parser I propose is implemented in…
▽ More
In this paper, I propose to implement syllabification in OT as a parser. I propose several innovations that result in a finite and small candidate set. The candidate set problem is handled with several moves: i) MAX and DEP violations are not hypothesized by the parser, ii) candidates are encoded locally, and iii) EVAL is applied constraint by constraint.
The parser I propose is implemented in Prolog. It has a number of desirable consequences. First, it runs and thus provides an existence proof that syllabification can be implemented in OT. There are a number of other desirable consequences as well. First, constraints are implemented as finite-state transducers. Second, the parser makes several interesting claims about the phonological properties of so-called nonrecoverable insertions and deletions. Third, the implementation suggests some particular reformulations of some of the benchmark constraints in the OT arsenal, e.g. *COMPLEX, PARSE, ONSET, and NOCODA.
△ Less
Submitted 14 October, 1997;
originally announced October 1997.
-
Syllable parsing in English and French
Authors:
Michael Hammond
Abstract:
In this paper I argue that Optimality Theory provides for an explanatory model of syllabic parsing in English and French. The argument is based on psycholinguistic facts that have been mysterious up to now. This argument is further buttressed by the computational implementation developed here. This model is important for several reasons. First, it provides a demonstration of how OT can be used i…
▽ More
In this paper I argue that Optimality Theory provides for an explanatory model of syllabic parsing in English and French. The argument is based on psycholinguistic facts that have been mysterious up to now. This argument is further buttressed by the computational implementation developed here. This model is important for several reasons. First, it provides a demonstration of how OT can be used in a performance domain. Second, it suggests a new relationship between phonological theory and psycholinguistics. (Code in Perl is included and a WWW-interface is running at http://mayo.douglass.arizona.edu.)
△ Less
Submitted 2 June, 1995;
originally announced June 1995.