Showing 1–31 of 31 results for author: Bonev, B

Coupled Ocean-Atmosphere Dynamics in a Machine Learning Earth System Model

Authors: Chenggong Wang, Michael S. Pritchard, Noah Brenowitz, Yair Cohen, Boris Bonev, Thorsten Kurth, Dale Durran, Jaideep Pathak

Abstract: Seasonal climate forecasts are socioeconomically important for managing the impacts of extreme weather events and for planning in sectors like agriculture and energy. Climate predictability on seasonal timescales is tied to boundary effects of the ocean on the atmosphere and coupled interactions in the ocean-atmosphere system. We present the Ocean-linked-atmosphere (Ola) model, a high-resolution (… ▽ More Seasonal climate forecasts are socioeconomically important for managing the impacts of extreme weather events and for planning in sectors like agriculture and energy. Climate predictability on seasonal timescales is tied to boundary effects of the ocean on the atmosphere and coupled interactions in the ocean-atmosphere system. We present the Ocean-linked-atmosphere (Ola) model, a high-resolution (0.25°) Artificial Intelligence/ Machine Learning (AI/ML) coupled earth-system model which separately models the ocean and atmosphere dynamics using an autoregressive Spherical Fourier Neural Operator architecture, with a view towards enabling fast, accurate, large ensemble forecasts on the seasonal timescale. We find that Ola exhibits learned characteristics of ocean-atmosphere coupled dynamics including tropical oceanic waves with appropriate phase speeds, and an internally generated El Niño/Southern Oscillation (ENSO) having realistic amplitude, geographic structure, and vertical structure within the ocean mixed layer. We present initial evidence of skill in forecasting the ENSO which compares favorably to the SPEAR model of the Geophysical Fluid Dynamics Laboratory. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Pretraining Codomain Attention Neural Operators for Solving Multiphysics PDEs

Authors: Md Ashiqur Rahman, Robert Joseph George, Mogab Elleithy, Daniel Leibovici, Zongyi Li, Boris Bonev, Colin White, Julius Berner, Raymond A. Yeh, Jean Kossaifi, Kamyar Azizzadenesheli, Anima Anandkumar

Abstract: Existing neural operator architectures face challenges when solving multiphysics problems with coupled partial differential equations (PDEs), due to complex geometries, interactions between physical variables, and the lack of large amounts of high-resolution training data. To address these issues, we propose Codomain Attention Neural Operator (CoDA-NO), which tokenizes functions along the codomain… ▽ More Existing neural operator architectures face challenges when solving multiphysics problems with coupled partial differential equations (PDEs), due to complex geometries, interactions between physical variables, and the lack of large amounts of high-resolution training data. To address these issues, we propose Codomain Attention Neural Operator (CoDA-NO), which tokenizes functions along the codomain or channel space, enabling self-supervised learning or pretraining of multiple PDE systems. Specifically, we extend positional encoding, self-attention, and normalization layers to the function space. CoDA-NO can learn representations of different PDE systems with a single model. We evaluate CoDA-NO's potential as a backbone for learning multiphysics PDEs over multiple systems by considering few-shot learning settings. On complex downstream tasks with limited data, such as fluid flow simulations and fluid-structure interactions, we found CoDA-NO to outperform existing methods on the few-shot learning task by over $36\%$. The code is available at https://github.com/ashiq24/CoDA-NO. △ Less

Submitted 5 April, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Neural Operators with Localized Integral and Differential Kernels

Authors: Miguel Liu-Schiaffini, Julius Berner, Boris Bonev, Thorsten Kurth, Kamyar Azizzadenesheli, Anima Anandkumar

Abstract: Neural operators learn mappings between function spaces, which is practical for learning solution operators of PDEs and other scientific modeling applications. Among them, the Fourier neural operator (FNO) is a popular architecture that performs global convolutions in the Fourier space. However, such global operations are often prone to over-smoothing and may fail to capture local details. In cont… ▽ More Neural operators learn mappings between function spaces, which is practical for learning solution operators of PDEs and other scientific modeling applications. Among them, the Fourier neural operator (FNO) is a popular architecture that performs global convolutions in the Fourier space. However, such global operations are often prone to over-smoothing and may fail to capture local details. In contrast, convolutional neural networks (CNN) can capture local features but are limited to training and inference at a single resolution. In this work, we present a principled approach to operator learning that can capture local features under two frameworks by learning differential operators and integral operators with locally supported kernels. Specifically, inspired by stencil methods, we prove that we obtain differential operators under an appropriate scaling of the kernel values of CNNs. To obtain local integral operators, we utilize suitable basis representations for the kernels based on discrete-continuous convolutions. Both these approaches preserve the properties of operator learning and, hence, the ability to predict at any resolution. Adding our layers to FNOs significantly improves their performance, reducing the relative L2-error by 34-72% in our experiments, which include a turbulent 2D Navier-Stokes and the spherical shallow water equations. △ Less

Submitted 8 June, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

Comments: Accepted at 2024 International Conference on Machine Learning

A Practical Probabilistic Benchmark for AI Weather Models

Authors: Noah D. Brenowitz, Yair Cohen, Jaideep Pathak, Ankur Mahesh, Boris Bonev, Thorsten Kurth, Dale R. Durran, Peter Harrington, Michael S. Pritchard

Abstract: Since the weather is chaotic, forecasts aim to predict the distribution of future states rather than make a single prediction. Recently, multiple data driven weather models have emerged claiming breakthroughs in skill. However, these have mostly been benchmarked using deterministic skill scores, and little is known about their probabilistic skill. Unfortunately, it is hard to fairly compare AI wea… ▽ More Since the weather is chaotic, forecasts aim to predict the distribution of future states rather than make a single prediction. Recently, multiple data driven weather models have emerged claiming breakthroughs in skill. However, these have mostly been benchmarked using deterministic skill scores, and little is known about their probabilistic skill. Unfortunately, it is hard to fairly compare AI weather models in a probabilistic sense, since variations in choice of ensemble initialization, definition of state, and noise injection methodology become confounding. Moreover, even obtaining ensemble forecast baselines is a substantial engineering challenge given the data volumes involved. We sidestep both problems by applying a decades-old idea -- lagged ensembles -- whereby an ensemble can be constructed from a moderately-sized library of deterministic forecasts. This allows the first parameter-free intercomparison of leading AI weather models' probabilistic skill against an operational baseline. The results reveal that two leading AI weather models, i.e. GraphCast and Pangu, are tied on the probabilistic CRPS metric even though the former outperforms the latter in deterministic scoring. We also reveal how multiple time-step loss functions, which many data-driven weather models have employed, are counter-productive: they improve deterministic metrics at the cost of increased dissipation, deteriorating probabilistic skill. This is confirmed through ablations applied to a spherical Fourier Neural Operator (SFNO) approach to AI weather forecasting. Separate SFNO ablations modulating effective resolution reveal it has a useful effect on ensemble dispersion relevant to achieving good ensemble calibration. We hope these and forthcoming insights from lagged ensembles can help guide the development of AI weather forecasts and have thus shared the diagnostic code. △ Less

Submitted 27 January, 2024; originally announced January 2024.

Comments: 15 pages, 5 figures

Advancing Parsimonious Deep Learning Weather Prediction using the HEALPix Mesh

Authors: Matthias Karlbauer, Nathaniel Cresswell-Clay, Dale R. Durran, Raul A. Moreno, Thorsten Kurth, Boris Bonev, Noah Brenowitz, Martin V. Butz

Abstract: We present a parsimonious deep learning weather prediction model to forecast seven atmospheric variables with 3-h time resolution for up to one-year lead times on a 110-km global mesh using the Hierarchical Equal Area isoLatitude Pixelization (HEALPix). In comparison to state-of-the-art (SOTA) machine learning (ML) weather forecast models, such as Pangu-Weather and GraphCast, our DLWP-HPX model us… ▽ More We present a parsimonious deep learning weather prediction model to forecast seven atmospheric variables with 3-h time resolution for up to one-year lead times on a 110-km global mesh using the Hierarchical Equal Area isoLatitude Pixelization (HEALPix). In comparison to state-of-the-art (SOTA) machine learning (ML) weather forecast models, such as Pangu-Weather and GraphCast, our DLWP-HPX model uses coarser resolution and far fewer prognostic variables. Yet, at one-week lead times, its skill is only about one day behind both SOTA ML forecast models and the SOTA numerical weather prediction model from the European Centre for Medium-Range Weather Forecasts. We report several improvements in model design, including switching from the cubed sphere to the HEALPix mesh, inverting the channel depth of the U-Net, and introducing gated recurrent units (GRU) on each level of the U-Net hierarchy. The consistent east-west orientation of all cells on the HEALPix mesh facilitates the development of location-invariant convolution kernels that successfully propagate weather patterns across the globe without requiring separate kernels for the polar and equatorial faces of the cube sphere. Without any loss of spectral power after the first two days, the model can be unrolled autoregressively for hundreds of steps into the future to generate realistic states of the atmosphere that respect seasonal trends, as showcased in one-year simulations. △ Less

Submitted 19 June, 2024; v1 submitted 11 September, 2023; originally announced November 2023.

Comments: Submitted to Journal of Advances in Modeling Earth Systems (JAMES)

ACE: A fast, skillful learned global atmospheric model for climate prediction

Authors: Oliver Watt-Meyer, Gideon Dresdner, Jeremy McGibbon, Spencer K. Clark, Brian Henn, James Duncan, Noah D. Brenowitz, Karthik Kashinath, Michael S. Pritchard, Boris Bonev, Matthew E. Peters, Christopher S. Bretherton

Abstract: Existing ML-based atmospheric models are not suitable for climate prediction, which requires long-term stability and physical consistency. We present ACE (AI2 Climate Emulator), a 200M-parameter, autoregressive machine learning emulator of an existing comprehensive 100-km resolution global atmospheric model. The formulation of ACE allows evaluation of physical laws such as the conservation of mass… ▽ More Existing ML-based atmospheric models are not suitable for climate prediction, which requires long-term stability and physical consistency. We present ACE (AI2 Climate Emulator), a 200M-parameter, autoregressive machine learning emulator of an existing comprehensive 100-km resolution global atmospheric model. The formulation of ACE allows evaluation of physical laws such as the conservation of mass and moisture. The emulator is stable for 100 years, nearly conserves column moisture without explicit constraints and faithfully reproduces the reference model's climate, outperforming a challenging baseline on over 90% of tracked variables. ACE requires nearly 100x less wall clock time and is 100x more energy efficient than the reference model using typically available resources. Without fine-tuning, ACE can stably generalize to a previously unseen historical sea surface temperature dataset. △ Less

Submitted 6 December, 2023; v1 submitted 3 October, 2023; originally announced October 2023.

Comments: Accepted at Tackling Climate Change with Machine Learning: workshop at NeurIPS 2023

Guaranteed Approximation Bounds for Mixed-Precision Neural Operators

Authors: Renbo Tu, Colin White, Jean Kossaifi, Boris Bonev, Nikola Kovachki, Gennady Pekhimenko, Kamyar Azizzadenesheli, Anima Anandkumar

Abstract: Neural operators, such as Fourier Neural Operators (FNO), form a principled approach for learning solution operators for PDEs and other mappings between function spaces. However, many real-world problems require high-resolution training data, and the training time and limited GPU memory pose big barriers. One solution is to train neural operators in mixed precision to reduce the memory requirement… ▽ More Neural operators, such as Fourier Neural Operators (FNO), form a principled approach for learning solution operators for PDEs and other mappings between function spaces. However, many real-world problems require high-resolution training data, and the training time and limited GPU memory pose big barriers. One solution is to train neural operators in mixed precision to reduce the memory requirement and increase training speed. However, existing mixed-precision training techniques are designed for standard neural networks, and we find that their direct application to FNO leads to numerical overflow and poor memory efficiency. Further, at first glance, it may appear that mixed precision in FNO will lead to drastic accuracy degradation since reducing the precision of the Fourier transform yields poor results in classical numerical solvers. We show that this is not the case; in fact, we prove that reducing the precision in FNO still guarantees a good approximation bound, when done in a targeted manner. Specifically, we build on the intuition that neural operator learning inherently induces an approximation error, arising from discretizing the infinite-dimensional ground-truth input function, implying that training in full precision is not needed. We formalize this intuition by rigorously characterizing the approximation and precision errors of FNO and bounding these errors for general input functions. We prove that the precision error is asymptotically comparable to the approximation error. Based on this, we design a simple method to optimize the memory-intensive half-precision tensor contractions by greedily finding the optimal contraction order. Through extensive experiments on different state-of-the-art neural operators, datasets, and GPUs, we demonstrate that our approach reduces GPU memory usage by up to 50% and improves throughput by 58% with little or no reduction in accuracy. △ Less

Submitted 5 May, 2024; v1 submitted 27 July, 2023; originally announced July 2023.

Comments: ICLR 2024

Spherical Fourier Neural Operators: Learning Stable Dynamics on the Sphere

Authors: Boris Bonev, Thorsten Kurth, Christian Hundt, Jaideep Pathak, Maximilian Baust, Karthik Kashinath, Anima Anandkumar

Abstract: Fourier Neural Operators (FNOs) have proven to be an efficient and effective method for resolution-independent operator learning in a broad variety of application areas across scientific machine learning. A key reason for their success is their ability to accurately model long-range dependencies in spatio-temporal data by learning global convolutions in a computationally efficient manner. To this… ▽ More Fourier Neural Operators (FNOs) have proven to be an efficient and effective method for resolution-independent operator learning in a broad variety of application areas across scientific machine learning. A key reason for their success is their ability to accurately model long-range dependencies in spatio-temporal data by learning global convolutions in a computationally efficient manner. To this end, FNOs rely on the discrete Fourier transform (DFT), however, DFTs cause visual and spectral artifacts as well as pronounced dissipation when learning operators in spherical coordinates since they incorrectly assume a flat geometry. To overcome this limitation, we generalize FNOs on the sphere, introducing Spherical FNOs (SFNOs) for learning operators on spherical geometries. We apply SFNOs to forecasting atmospheric dynamics, and demonstrate stable auto\-regressive rollouts for a year of simulated time (1,460 steps), while retaining physically plausible dynamics. The SFNO has important implications for machine learning-based simulation of climate dynamics that could eventually help accelerate our response to climate change. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Gas Sources from the Coma and Nucleus of Comet 46P/Wirtanen Observed Using ALMA

Authors: M. A. Cordiner, N. X. Roth, S. N. Milam, G. Villanueva, D. Bockelee-Morvan, A. J. Remijan, S. B. Charnley, N. Biver, D. C. Lis, C. Qi, B. Bonev, J. Crovisier, J. Boissier

Abstract: Gas-phase molecules in cometary atmospheres (comae) originate primarily from (1) outgassing by the nucleus, (2) sublimation of icy grains in the near-nucleus coma, and (3) coma (photo-)chemical processes. However, the majority of cometary gases observed at radio wavelengths have yet to be mapped, so their production/release mechanisms remain uncertain. Here we present observations of six molecular… ▽ More Gas-phase molecules in cometary atmospheres (comae) originate primarily from (1) outgassing by the nucleus, (2) sublimation of icy grains in the near-nucleus coma, and (3) coma (photo-)chemical processes. However, the majority of cometary gases observed at radio wavelengths have yet to be mapped, so their production/release mechanisms remain uncertain. Here we present observations of six molecular species towards comet 46P/Wirtanen, obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) during the comet's unusually close (~0.1 au) approach to Earth in December 2018. Interferometric maps of HCN, CH3OH, CH3CN, H2CO, CS and HNC were obtained at an unprecedented sky-projected spatial resolution of up to 25 km, enabling the nucleus and coma sources of these molecules to be accurately quantified. The HCN, CH3OH and CH3CN spatial distributions are consistent with production by direct outgassing from (or very near to) the nucleus, with a significant proportion of the observed CH3OH originating from sublimation of icy grains in the near-nucleus coma (at a scale-length $L_p=36\pm7$ km). On the other hand, H2CO, CS and HNC originate primarily from distributed coma sources (with $L_p$ values in the range 550-16,000 km), the identities of which remain to be established. The HCN, CH3OH and HNC abundances in 46P are consistent with the average values previously observed in comets, whereas the H2CO, CH3CN and CS abundances are relatively low. △ Less

Submitted 19 June, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: Accepted for publication in the Astrophysical Journal on 2023-06-17

Molecular Outgassing in Centaur 29P/Schwassmann-Wachmann 1 During Its Exceptional 2021 Outburst: Coordinated Multi-Wavelength Observations Using nFLASH at APEX and iSHELL at the NASA-IRTF

Authors: Nathan X. Roth, Stefanie N. Milam, Michael A. DiSanti, Geronimo L. Villanueva, Sara Faggi, Boncho P. Bonev, Martin A. Cordiner, Anthony J. Remijan, Dominique Bockelée-Morvan, Nicolas Biver, Jacques Crovisier, Dariusz C. Lis, Steven B. Charnley, Emmanuel Jehin, Eva. S. Wirström, Adam J. McKay

Abstract: The extraordinary 2021 September-October outburst of Centaur 29P/Schwassmann-Wachmann 1 afforded an opportunity to test the composition of primitive Kuiper disk material at high sensitivity. We conducted nearly simultaneous multi-wavelength spectroscopic observations of 29P/Schwassmann-Wachmann 1 using iSHELL at the NASA Infrared Telescope Facility and nFLASH at the Atacama Pathfinder EXperiment (… ▽ More The extraordinary 2021 September-October outburst of Centaur 29P/Schwassmann-Wachmann 1 afforded an opportunity to test the composition of primitive Kuiper disk material at high sensitivity. We conducted nearly simultaneous multi-wavelength spectroscopic observations of 29P/Schwassmann-Wachmann 1 using iSHELL at the NASA Infrared Telescope Facility and nFLASH at the Atacama Pathfinder EXperiment (APEX) on 2021 October 6, with follow-up APEX/nFLASH observations on 2021 October 7 and 2022 April 3. This coordinated campaign between near-infrared and radio wavelengths enabled us to sample molecular emission from a wealth of coma molecules and to perform measurements that cannot be accomplished with either wavelength alone. We securely detected CO emission on all dates with both facilities, including velocity-resolved spectra of the CO (J=2-1) transition with APEX/nFLASH and multiple CO (v=1-0) rovibrational transitions with IRTF/iSHELL. We report rotational temperatures, coma kinematics, and production rates for CO and stringent (3-sigma) upper limits on abundance ratios relative to CO for CH4, C2H6, CH3OH, H2CO, CS, and OCS. Our upper limits for CS/CO and OCS/CO represent their first values in the literature for this Centaur. Upper limits for CH4, C2H6, CH3OH, and H2CO are the most stringent reported to date, and are most similar to values found in ultra CO-rich Oort cloud comet C/2016 R2 (PanSTARRS), which may have implications for how ices are preserved in cometary nuclei. We demonstrate the superb synergy of coordinated radio and near-infrared measurements, and advocate for future small body studies that jointly leverage the capabilities of each wavelength. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Radiative processes as diagnostics of cometary atmospheres

Authors: D. Bodewits, B. P. Bonev, M. A. Cordiner, G. L. Villanueva

Abstract: In this chapter, we provide a review of radiative processes in cometary atmospheres spanning a broad range of wavelengths, from radio to X-rays. We focus on spectral modeling, observational opportunities, and anticipated challenges in the interpretation of new observations, based on our current understanding of the atomic and molecular processes occurring in the atmospheres of small, icy bodies. C… ▽ More In this chapter, we provide a review of radiative processes in cometary atmospheres spanning a broad range of wavelengths, from radio to X-rays. We focus on spectral modeling, observational opportunities, and anticipated challenges in the interpretation of new observations, based on our current understanding of the atomic and molecular processes occurring in the atmospheres of small, icy bodies. Close to the surface, comets possess a thermalized atmosphere that traces the irregular shape of the nucleus. Gravity is too low to retain the gas, which flows out to form a large, collisionless exosphere (coma) that interacts with the heliospheric radiation environment. As such, cometary comae represent conditions that are familiar in the context of planetary atmosphere studies. However, the outer comae are tenuous, with densities lower than those found in vacuum chambers on Earth. Comets, therefore, provide us with unique natural laboratories that can be understood using state-of-the-art theoretical treatments of the relevant microphysical processes. Radiative processes offer direct diagnostics of the local physical conditions, as well as the macroscopic coma properties.These can be used to improve our understanding of comets and other astrophysical environments such as icy moons and the interstellar medium. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: Review for Comets III book

Comets in context: Comparing comet compositions with protosolar nebula models

Authors: Karen Willacy, Neal Turner, Boncho Bonev, Erika Gibb, Neil Dello Russo, Michael DiSanti, Ronald J. Vervack Jr., Nathan X. Roth

Abstract: Comets provide a valuable window into the chemical and physical conditions at the time of their formation in the young solar system. We seek insights into where and when these objects formed by comparing the range of abundances observed for nine molecules and their average values across a sample of 29 comets to the predicted midplane ice abundances from models of the protosolar nebula. Our fiducia… ▽ More Comets provide a valuable window into the chemical and physical conditions at the time of their formation in the young solar system. We seek insights into where and when these objects formed by comparing the range of abundances observed for nine molecules and their average values across a sample of 29 comets to the predicted midplane ice abundances from models of the protosolar nebula. Our fiducial model, where ices are inherited from the interstellar medium, can account for the observed mixing ratio ranges of each molecule considered, but no single location or time reproduces the abundances of all molecules simultaneously. This suggests that each comet consists of material processed under a range of conditions. In contrast, a model where the initial composition of disk material is `reset', wiping out any previous chemical history, cannot account for the complete range of abundances observed in comets. Using toy models that combine material processed under different thermal conditions we find that a combination of warm (CO-poor) and cold (CO-rich) material is required to account for both the average properties of the Jupiter-family and Oort cloud comets, and the individual comets we consider. This could occur by the transport (either radial or vertical) of ice-coated dust grains in the early solar system. Comparison of the models to the average Jupiter-family and Oort cloud comet compositions suggest the two families formed in overlapping regions of the disk, in agreement with the findings of A'Hearn et al. (2012) and with the predictions of the Nice model (Gomes et al. 2005, Tsiganis et al. 2005). △ Less

Submitted 15 April, 2022; originally announced April 2022.

Comments: 32 pages, 8 figures

A hierarchical preconditioner for wave problems in quasilinear complexity

Authors: Boris Bonev, Jan S. Hesthaven

Abstract: The paper introduces a novel, hierarchical preconditioner based on nested dissection and hierarchical matrix compression. The preconditioner is intended for continuous and discontinuous Galerkin formulations of elliptic problems. We exploit the property that Schur complements arising in such problems can be well approximated by hierarchical matrices. An approximate factorization can be computed ma… ▽ More The paper introduces a novel, hierarchical preconditioner based on nested dissection and hierarchical matrix compression. The preconditioner is intended for continuous and discontinuous Galerkin formulations of elliptic problems. We exploit the property that Schur complements arising in such problems can be well approximated by hierarchical matrices. An approximate factorization can be computed matrix-free and in a (quasi-)linear number of operations. The nested dissection is specifically designed to aid the factorization process using hierarchical matrices. We demonstrate the viability of the preconditioner on a range of 2D problems, including the Helmholtz equation and the elastic wave equation. Throughout all tests, including wave phenomena with high wavenumbers, the generalized minimal residual method (GMRES) with the proposed preconditioner converges in a very low number of iterations. We demonstrate that this is due to the hierarchical nature of our approach which makes the high wavenumber limit manageable. △ Less

Submitted 17 May, 2021; originally announced May 2021.

MSC Class: 65F08; 15A23; 65F30; 65F50

Rapidly Varying Anisotropic Methanol (CH$_3$OH) Production in the Inner Coma of Comet 46P/Wirtanen as Revealed by the ALMA Atacama Compact Array

Authors: Nathan X. Roth, Stefanie N. Milam, Martin A. Cordiner, Dominique Bockelée-Morvan, Michael A. DiSanti, Jérémie Boissier, Nicolas Biver, Jacques Crovisier, Neil Dello Russo, Boncho P. Bonev, Chunhua Qi, Anthony J. Remijan, Steven B. Charnley, Erika L. Gibb, Miguel de Val-Borro, Emmanuël Jehin

Abstract: We report the detection of CH$_3$OH emission in comet 46P/Wirtanen on UT 2018 December 8 and 9 using the Atacama Compact Array (ACA), part of the Atacama Large Millimeter/Submillimeter Array (ALMA). These interferometric measurements of CH$_3$OH along with continuum emission from dust probed the inner coma ($<$2000 km from the nucleus) of 46P/Wirtanen approximately one week before its closest appr… ▽ More We report the detection of CH$_3$OH emission in comet 46P/Wirtanen on UT 2018 December 8 and 9 using the Atacama Compact Array (ACA), part of the Atacama Large Millimeter/Submillimeter Array (ALMA). These interferometric measurements of CH$_3$OH along with continuum emission from dust probed the inner coma ($<$2000 km from the nucleus) of 46P/Wirtanen approximately one week before its closest approach to Earth ($Δ$ = 0.089 -- 0.092 au), revealing rapidly varying and anisotropic CH$_3$OH outgassing during five separate ACA executions between UT 23:57 December 7 and UT 04:55 December 9, with a clear progression in the spectral line profiles over a timescale of minutes. We present spectrally integrated flux maps, production rates, rotational temperatures, and spectral line profiles of CH$_3$OH during each ACA execution. The variations in CH$_3$OH outgassing are consistent with Wirtanen's 9 hr nucleus rotational period derived from optical and millimeter wavelength measurements and thus are likely coupled to the changing illumination of active sites on the nucleus. The consistent blue offset of the line center indicates enhanced CH$_3$OH sublimation from the sunward hemisphere of the comet, perhaps from icy grains. These results demonstrate the exceptional capabilities of the ACA for time-resolved measurements of comets such as 46P/Wirtanen. △ Less

Submitted 11 March, 2021; originally announced March 2021.

Molecular composition of comet 46P/Wirtanen from millimetre-wave spectroscopy

Authors: N. Biver, D. Bockelée-Morvan, J. Boissier, R. Moreno, J. Crovisier, D. C. Lis, P. Colom, M. Cordiner, S. Milam, N. X. Roth, B. P. Bonev, N. Dello Russo, R. Vervack, M. A. DiSanti

Abstract: We present the results of a molecular survey of comet 46P/Wirtanen undertaken with the IRAM 30-m and NOEMA radio telescopes in December 2018. Observations at IRAM 30-m during the 12-18 Dec. period comprise a 2 mm spectral survey covering 25 GHz and a 1 mm survey covering 62 GHz. The gas outflow velocity and kinetic temperature have been accurately constrained by the observations. We derive abundan… ▽ More We present the results of a molecular survey of comet 46P/Wirtanen undertaken with the IRAM 30-m and NOEMA radio telescopes in December 2018. Observations at IRAM 30-m during the 12-18 Dec. period comprise a 2 mm spectral survey covering 25 GHz and a 1 mm survey covering 62 GHz. The gas outflow velocity and kinetic temperature have been accurately constrained by the observations. We derive abundances of 11 molecules, some being identified remotely for the first time in a Jupiter-family comet, including complex organic molecules such as formamide, ethylene glycol, acetaldehyde, or ethanol. Sensitive upper limits on the abundances of 24 other molecules are obtained. The comet is found to be relatively rich in methanol (3.4 percent relative to water), but relatively depleted in CO, CS, HNC, HNCO, and HCOOH. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Comments: 35 pages, 16 figures, to be published in Astronomy & Astrophysics

Journal ref: A&A 648, A49 (2021)

The Crucial Role of Ground- and Space-Based Remote Sensing Studies of Cometary Volatiles in the Next Decade (2023-2032)

Authors: Nathan X. Roth, Dennis Bodewits, Boncho Bonev, Anita Cochran, Michael Combi, Martin Cordiner, Neil Dello Russo, Michael DiSanti, Sara Faggi, Lori Feaga, Yan Fernandez, Manuela Lippi, Adam McKay, Matthew Knight, Stefanie Milam, John W. Noonan, Anthony Remijan, Geronimo Villanueva

Abstract: The study of comets affords a unique window into the birth, infancy, and subsequent history of the solar system. There is strong evidence that comets incorporated pristine interstellar material as well as processed nebular matter, providing insights into the composition and prevailing conditions over wide swaths of the solar nebula at the time of planet formation. Dynamically new Oort cloud comets… ▽ More The study of comets affords a unique window into the birth, infancy, and subsequent history of the solar system. There is strong evidence that comets incorporated pristine interstellar material as well as processed nebular matter, providing insights into the composition and prevailing conditions over wide swaths of the solar nebula at the time of planet formation. Dynamically new Oort cloud comets harbor primitive ices that have been stored thousands of astronomical units from the Sun and have suffered minimal thermal or radiative processing since their emplacement ~4.5 Gyr ago. Periodic, more dynamically evolved comets such as the Halley-type and Jupiter-family comets reveal the effects of lives spent over a range of heliocentric distances, including perihelion passages into the very inner solar system. Systematically characterizing the information imprinted in the native ice compositions of these objects is critical to understanding the formation and evolution of the solar system, the presence of organic matter and water on the terrestrial planets, the chemistry present in protoplanetary disks around other stars, and the nature of interstellar interlopers such as 2I/Borisov. Although comet rendezvous and sample return missions can provide remarkable insights into the properties of a few short-period comets, the on-sky capacity necessary to perform population-level comet studies while simultaneously remaining sensitive to the paradigm-challenging science that individual comets can reveal can only be provided by remote sensing observations. Here we report the state-of-the-art in ground- and space-based remote sensing of cometary volatiles, review the remarkable progress of the previous decade, articulate the pressing questions that ground- and space-based work will address over the next ten years, and advocate for the technology and resources necessary to realize these aspirations. △ Less

Submitted 16 July, 2020; originally announced July 2020.

Comments: White paper for the National Academies of Sciences, Engineering, and Medicine Planetary Science and Astrobiology Decadal Survey 2023-2032

The Peculiar Volatile Composition of CO-Dominated Comet C/2016 R2 (PanSTARRS)

Authors: Adam McKay, Michael DiSanti, Michael Kelley, Matthew Knight, Maria Womack, Kacper Wierzchos, Olga Harrington-Pinto, Boncho Bonev, Geronimo Villanueva, Neil Dello Russo, Anita Cochran, Nicolas Biver, James Bauer, Ronald Vervack, Jr., Erika Gibbs, Nathan Roth, Hideyo Kawakita

Abstract: Comet C/2016 R2 (PanSTARRS) has a peculiar volatile composition, with CO being the dominant volatile as opposed to H$_2$O and one of the largest N$_2$/CO ratios ever observed in a comet. Using observations obtained with the \textit{Spitzer Space Telescope}, NASA's Infrared Telescope Facility, the 3.5-meter ARC telescope at Apache Point Observatory, the Discovery Channel Telescope at Lowell Observa… ▽ More Comet C/2016 R2 (PanSTARRS) has a peculiar volatile composition, with CO being the dominant volatile as opposed to H$_2$O and one of the largest N$_2$/CO ratios ever observed in a comet. Using observations obtained with the \textit{Spitzer Space Telescope}, NASA's Infrared Telescope Facility, the 3.5-meter ARC telescope at Apache Point Observatory, the Discovery Channel Telescope at Lowell Observatory, and the Arizona Radio Observatory 10-m Submillimeter Telescope we quantified the abundances of 12 different species in the coma of R2 PanSTARRS. We confirm the high abundances of CO and N$_2$ and heavy depletions of H$_2$O, HCN, CH$_3$OH, and H$_2$CO compared to CO reported by previous studies. We provide the first measurements (or most sensitive measurements/constraints) on H$_2$O, CO$_2$, CH$_4$, C$_2$H$_6$, OCS, C$_2$H$_2$, and NH$_3$, all of which are depleted relative to CO by at least one to two orders of magnitude compared to values commonly observed in comets. The observed species also show strong enhancements relative to H$_2$O, and even when compared to other species like CH$_4$ or CH$_3$OH most species show deviations from typical comets by at least a factor of two to three. The only mixing ratios found to be close to typical are CH$_3$OH/CO$_2$ and CH$_3$OH/CH$_4$. While R2 PanSTARRS was located at a heliocentric distance of 2.8 AU at the time of our observations in January/February 2018, we argue that this alone cannot account for the peculiar observed composition of this comet and therefore must reflect its intrinsic composition. We discuss possible implications for this clear outlier in compositional studies of comets obtained to date, and encourage future dynamical and chemical modeling in order to better understand what the composition of R2 PanSTARRS tells us about the early Solar System. △ Less

Submitted 16 July, 2019; originally announced July 2019.

Comments: 36 pages, 12 figures, 7 tables, accepted for publication in AJ

The 67P/Churyumov-Gerasimenko observation campaign in support of the Rosetta mission

Authors: C. Snodgrass, M. F. A'Hearn, F. Aceituno, V. Afanasiev, S. Bagnulo, J. Bauer, G. Bergond, S. Besse, N. Biver, D. Bodewits, H. Boehnhardt, B. P. Bonev, G. Borisov, B. Carry, V. Casanova, A. Cochran, B. C. Conn, B. Davidsson, J. K. Davies, J. de León, E. de Mooij, M. de Val-Borro, M. Delacruz, M. A. DiSanti, J. E. Drew , et al. (90 additional authors not shown)

Abstract: We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond t… ▽ More We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond the spacecraft, and a way to directly compare 67P with other comets. The observations revealed 67P to be a relatively `well behaved' comet, typical of Jupiter family comets and with activity patterns that repeat from orbit-to-orbit. Comparison between this large collection of telescopic observations and the in situ results from Rosetta will allow us to better understand comet coma chemistry and structure. This work is just beginning as the mission ends -- in this paper we present a summary of the ground-based observations and early results, and point to many questions that will be addressed in future studies. △ Less

Submitted 30 May, 2017; originally announced May 2017.

Comments: Author prepared version; final published version available at journal. 22 pages

Journal ref: Phil. Trans. R. Soc. A 375, 20160249 (2017)

Generating Liquid Simulations with Deformation-aware Neural Networks

Authors: Lukas Prantl, Boris Bonev, Nils Thuerey

Abstract: We propose a novel approach for deformation-aware neural networks that learn the weighting and synthesis of dense volumetric deformation fields. Our method specifically targets the space-time representation of physical surfaces from liquid simulations. Liquids exhibit highly complex, non-linear behavior under changing simulation conditions such as different initial conditions. Our algorithm captur… ▽ More We propose a novel approach for deformation-aware neural networks that learn the weighting and synthesis of dense volumetric deformation fields. Our method specifically targets the space-time representation of physical surfaces from liquid simulations. Liquids exhibit highly complex, non-linear behavior under changing simulation conditions such as different initial conditions. Our algorithm captures these complex phenomena in two stages: a first neural network computes a weighting function for a set of pre-computed deformations, while a second network directly generates a deformation field for refining the surface. Key for successful training runs in this setting is a suitable loss function that encodes the effect of the deformations, and a robust calculation of the corresponding gradients. To demonstrate the effectiveness of our approach, we showcase our method with several complex examples of flowing liquids with topology changes. Our representation makes it possible to rapidly generate the desired implicit surfaces. We have implemented a mobile application to demonstrate that real-time interactions with complex liquid effects are possible with our approach. △ Less

Submitted 20 February, 2019; v1 submitted 25 April, 2017; originally announced April 2017.

Comments: ICLR 2019, further information and videos at https://ge.in.tum.de/publications/2017-prantl-defonn/

PASCAL Boundaries: A Class-Agnostic Semantic Boundary Dataset

Authors: Vittal Premachandran, Boyan Bonev, Alan L. Yuille

Abstract: In this paper, we address the boundary detection task motivated by the ambiguities in current definition of edge detection. To this end, we generate a large database consisting of more than 10k images (which is 20x bigger than existing edge detection databases) along with ground truth boundaries between 459 semantic classes including both foreground objects and different types of background, and c… ▽ More In this paper, we address the boundary detection task motivated by the ambiguities in current definition of edge detection. To this end, we generate a large database consisting of more than 10k images (which is 20x bigger than existing edge detection databases) along with ground truth boundaries between 459 semantic classes including both foreground objects and different types of background, and call it the PASCAL Boundaries dataset, which will be released to the community. In addition, we propose a novel deep network-based multi-scale semantic boundary detector and name it Multi-scale Deep Semantic Boundary Detector (M-DSBD). We provide baselines using models that were trained on edge detection and show that they transfer reasonably to the task of boundary detection. Finally, we point to various important research problems that this dataset can be used for. △ Less

Submitted 25 November, 2015; originally announced November 2015.

Ground-truth dataset and baseline evaluations for image base-detail separation algorithms

Authors: Xuan Dong, Boyan Bonev, Weixin Li, Weichao Qiu, Xianjie Chen, Alan Yuille

Abstract: Base-detail separation is a fundamental computer vision problem consisting of modeling a smooth base layer with the coarse structures, and a detail layer containing the texture-like structures. One of the challenges of estimating the base is to preserve sharp boundaries between objects or parts to avoid halo artifacts. Many methods have been proposed to address this problem, but there is no ground… ▽ More Base-detail separation is a fundamental computer vision problem consisting of modeling a smooth base layer with the coarse structures, and a detail layer containing the texture-like structures. One of the challenges of estimating the base is to preserve sharp boundaries between objects or parts to avoid halo artifacts. Many methods have been proposed to address this problem, but there is no ground-truth dataset of real images for quantitative evaluation. We proposed a procedure to construct such a dataset, and provide two datasets: Pascal Base-Detail and Fashionista Base-Detail, containing 1000 and 250 images, respectively. Our assumption is that the base is piecewise smooth and we label the appearance of each piece by a polynomial model. The pieces are objects and parts of objects, obtained from human annotations. Finally, we proposed a way to evaluate methods with our base-detail ground-truth and we compared the performances of seven state-of-the-art algorithms. △ Less

Submitted 18 February, 2016; v1 submitted 20 November, 2015; originally announced November 2015.

Comments: This paper has been withdrawn by the author due to some un-proper examples

Statistical sentiment analysis performance in Opinum

Authors: Boyan Bonev, Gema Ramírez-Sánchez, Sergio Ortiz Rojas

Abstract: The classification of opinion texts in positive and negative is becoming a subject of great interest in sentiment analysis. The existence of many labeled opinions motivates the use of statistical and machine-learning methods. First-order statistics have proven to be very limited in this field. The Opinum approach is based on the order of the words without using any syntactic and semantic informati… ▽ More The classification of opinion texts in positive and negative is becoming a subject of great interest in sentiment analysis. The existence of many labeled opinions motivates the use of statistical and machine-learning methods. First-order statistics have proven to be very limited in this field. The Opinum approach is based on the order of the words without using any syntactic and semantic information. It consists of building one probabilistic model for the positive and another one for the negative opinions. Then the test opinions are compared to both models and a decision and confidence measure are calculated. In order to reduce the complexity of the training corpus we first lemmatize the texts and we replace most named-entities with wildcards. Opinum presents an accuracy above 81% for Spanish opinions in the financial products domain. In this work we discuss which are the most important factors that have impact on the classification performance. △ Less

Submitted 2 March, 2013; originally announced March 2013.

The chemical composition of CO-rich comet C/2009 P1 (Garradd) at Rh = 2.4 and 2.0 AU before perihelion

Authors: L. Paganini, M. J. Mumma, G. L. Villanueva, M. A. DiSanti, B. P. Bonev, M. Lippi, H. Boehnhardt

Abstract: We quantified ten parent volatiles in comet C/2009 P1 (Garradd) before perihelion, through high-dispersion infrared spectra acquired with CRIRES at ESO's VLT on UT 2011 August 07 (Rh = 2.4 AU) and September 17-21 (Rh = 2.0 AU). On August 07, water was searched but not detected at an upper limit (3σ) of 2.1 \times 10^28 s-1, while ethane was detected with a production rate of 6.1 \times 10^26 s-1 (… ▽ More We quantified ten parent volatiles in comet C/2009 P1 (Garradd) before perihelion, through high-dispersion infrared spectra acquired with CRIRES at ESO's VLT on UT 2011 August 07 (Rh = 2.4 AU) and September 17-21 (Rh = 2.0 AU). On August 07, water was searched but not detected at an upper limit (3σ) of 2.1 \times 10^28 s-1, while ethane was detected with a production rate of 6.1 \times 10^26 s-1 (apparent mixing ratio > 2.90%). On September 17-21, the mean production rate for water was 8.4 \times 10^28 s-1, and abundance ratios (relative to water) of detected trace species were: CO (12.51%), CH3OH (3.90%), CH4 (1.24%), C2H6 (1.01%) and HCN (0.36%). Upper limits (3σ) to abundances for four minor species were: NH3 (1.55%), C2H2 (0.13%), HDO (0.89%) and OCS (0.20%). Given the relatively large heliocentric distance, we explored the effect of water not being fully sublimated within our FOV and identified the 'missing' water fraction needed to reconcile the retrieved abundance ratios with the mean values found for "organics-normal". The individual spatial profiles of parent volatiles and the continuum displayed rather asymmetric outgassing. Indications of H2O and CO gas being released in different directions suggest different active vents and/or the possible existence of polar and apolar ice aggregates in the nucleus. The high fractional abundance of CO identifies comet C/2009 P1 as a CO-rich comet. △ Less

Submitted 1 February, 2012; originally announced February 2012.

Comments: To appear in: The Astrophysical Journal Letters (LET27048R)

The formation heritage of Jupiter Family Comet 10P/Tempel 2 as revealed by infrared spectroscopy

Authors: L. Paganini, M. J. Mumma, B. P. Bonev, G. L. Villanueva, M. A. DiSanti, J. V. Keane, K. J. Meech

Abstract: We present spectral and spatial information for major volatile species in Comet 10P/Tempel 2, based on high-dispersion infrared spectra acquired on UT 2010 July 26 (heliocentric distance Rh = 1.44 AU) and September 18 (Rh = 1.62 AU), following the comet's perihelion passage on UT 2010 July 04. The total production rate for water on July 26 was (1.90 +/- 0.12) x 10^28 molecules s-1, and abundances… ▽ More We present spectral and spatial information for major volatile species in Comet 10P/Tempel 2, based on high-dispersion infrared spectra acquired on UT 2010 July 26 (heliocentric distance Rh = 1.44 AU) and September 18 (Rh = 1.62 AU), following the comet's perihelion passage on UT 2010 July 04. The total production rate for water on July 26 was (1.90 +/- 0.12) x 10^28 molecules s-1, and abundances of six trace gases (relative to water) were: CH3OH (1.58% +/- 0.23), C2H6 (0.39% +/- 0.04), NH3 (0.83% +/- 0.20), and HCN (0.13% +/- 0.02). A detailed analysis of intensities for water emission lines provided a rotational temperature of 35 +/- 3 K. The mean OPR is consistent with nuclear spin populations in statistical equilibrium (OPR = 3.01 +/- 0.18), and the (1-sigma) lower bound corresponds to a spin temperature > 38 K. Our measurements were contemporaneous with a jet-like feature observed at optical wavelengths. The spatial profiles of four primary volatiles display strong enhancements in the jet direction, which favors release from a localized vent on the nucleus. The measured IR continuum is much more sharply peaked and is consistent with a dominant contribution from the nucleus itself. The peak intensities for H2O, CH3OH, and C2H6 are offset by ~200 km in the jet direction, suggesting the possible existence of a distributed source, such as the release of icy grains that subsequently sublimed in the coma. On UT September 18, no obvious emission lines were present in our spectra, nevertheless we obtained a 3-sigma upper limit Q(H2O) < 2.86 x 10^27 molecules s-1. △ Less

Submitted 9 February, 2012; v1 submitted 24 January, 2012; originally announced January 2012.

Temporal and Spatial Aspects of Gas Release During the 2010 Apparition of Comet 103P/Hartley-2

Authors: M. J. Mumma, B. P. Bonev, G. L. Villanueva, L. Paganini, M. A. DiSanti, E. L. Gibb, J. V. Keane, K. J. Meech, G. A. Blake, R. S. Ellis, M. Lippi, H. Bõhnhardt, K. Magee-Sauer

Abstract: We report measurements of eight primary volatiles (H2O, HCN, CH4, C2H6, CH3OH, C2H2, H2CO, and NH3) and two product species (OH and NH2) in comet 103P/Hartley-2 using high dispersion infrared spectroscopy. We quantified the long- and short-term behavior of volatile release over a three-month interval that encompassed the comet's close approach to Earth, its perihelion passage, and flyby of the com… ▽ More We report measurements of eight primary volatiles (H2O, HCN, CH4, C2H6, CH3OH, C2H2, H2CO, and NH3) and two product species (OH and NH2) in comet 103P/Hartley-2 using high dispersion infrared spectroscopy. We quantified the long- and short-term behavior of volatile release over a three-month interval that encompassed the comet's close approach to Earth, its perihelion passage, and flyby of the comet by the Deep Impact spacecraft during the EPOXI mission. We present production rates for individual species, their mixing ratios relative to water, and their spatial distributions in the coma on multiple dates. The production rates for water, ethane, HCN, and methanol vary in a manner consistent with independent measures of nucleus rotation, but mixing ratios for HCN, C2H6, & CH3OH are independent of rotational phase. Our results demonstrate that the ensemble average composition of gas released from the nucleus is well defined, and relatively constant over the three-month interval (September 18 through December 17). If individual vents vary in composition, enough diverse vents must be active simultaneously to approximate (in sum) the bulk composition of the nucleus. The released primary volatiles exhibit diverse spatial properties which favor the presence of separate polar and apolar ice phases in the nucleus, establish dust and gas release from icy clumps (and also, directly from the nucleus), and provide insights into the driver for the cyanogen (CN) polar jet. The spatial distributions of C2H6 & HCN along the near-polar jet (UT 19.5 October) and nearly orthogonal to it (UT 22.5 October) are discussed relative to the origin of CN. The ortho-para ratio (OPR) of water was 2.85 \pm 0.20; the lower bound (2.65) defines Tspin > 32 K. These values are consistent with results returned from ISO in 1997. △ Less

Submitted 24 April, 2011; originally announced April 2011.

Comments: 18 pages, 3 figures, to be published in: Astrophysical Journal Letters

A Multi-Wavelength Study of Parent Volatile Abundances in Comet C/2006 M4 (SWAN)

Authors: Michael A. Disanti, Geronimo L. Villanueva, Stefanie N. Milam, Lindsay N. Zack, Boncho P. Bonev, Michael J. Mumma, Lucy M. Ziurys, William M. Anderson

Abstract: Volatile organic emissions were detected post-perihelion in the long period comet C/2006 M4 (SWAN) in October and November 2006. Our study combines target-of-opportunity observations using the infrared Cryogenic Echelle Spectrometer (CSHELL) at the NASA-IRTF 3-m telescope, and millimeter wavelength observations using the Arizona Radio Observatory (ARO) 12-m telescope. Five parent volatiles were… ▽ More Volatile organic emissions were detected post-perihelion in the long period comet C/2006 M4 (SWAN) in October and November 2006. Our study combines target-of-opportunity observations using the infrared Cryogenic Echelle Spectrometer (CSHELL) at the NASA-IRTF 3-m telescope, and millimeter wavelength observations using the Arizona Radio Observatory (ARO) 12-m telescope. Five parent volatiles were measured with CSHELL (H2O, CO, CH3OH, CH4, and C2H6), and two additional species (HCN and CS) were measured with the ARO 12-m. These revealed highly depleted CO and somewhat enriched CH3OH compared with abundances observed in the dominant group of long-period (Oort cloud) comets in our sample and similar to those observed recently in Comet 8P/Tuttle. This may indicate highly efficient H-atom addition to CO at very low temperature (~ 10 - 20 K) on the surfaces of interstellar (pre-cometary) grains. Comet C/2006 M4 had nearly "normal" C2H6 and CH4, suggesting a processing history similar to that experienced by the dominant group. When compared with estimated water production at the time of the millimeter observations, HCN was slightly depleted compared with the normal abundance in comets based on IR observations but was consistent with the majority of values from the millimeter. The ratio CS/HCN in C/2006 M4 was within the range measured in ten comets at millimeter wavelengths. The higher apparent H-atom conversion efficiency compared with most comets may indicate that the icy grains incorporated into C/2006 M4 were exposed to higher H-atom densities, or alternatively to similar densities but for a longer period of time. △ Less

Submitted 12 May, 2009; originally announced May 2009.

Comments: Accepted for Publication in Icarus (27 pages, 4 tables, 4 figures)

Comet C/2004 Q2 (MACHHOLZ): Parent Volatiles, a Search for Deuterated Methane, and Constraint on the CH4 Spin Temperature

Authors: Boncho P. Bonev, Michael J. Mumma, Erika L. Gibb, Michael A. Disanti, Geronimo L. Villanueva, Karen Magee-Sauer, Richard S. Ellis

Abstract: High-dispersion (l/dl ~ 25,000) infrared spectra of Comet C/2004 Q2 (Machholz) were acquired on Nov. 28-29, 2004, and Jan. 19, 2005 (UT dates) with NIRSPEC at the Keck-2 telescope on Mauna Kea. We detected H2O, CH4, C2H2, C2H6, CO, H2CO, CH3OH, HCN, and NH3 and we conducted a sensitive search for CH3D. We report rotational temperatures, production rates, and mixing ratios (with respect to H2O) a… ▽ More High-dispersion (l/dl ~ 25,000) infrared spectra of Comet C/2004 Q2 (Machholz) were acquired on Nov. 28-29, 2004, and Jan. 19, 2005 (UT dates) with NIRSPEC at the Keck-2 telescope on Mauna Kea. We detected H2O, CH4, C2H2, C2H6, CO, H2CO, CH3OH, HCN, and NH3 and we conducted a sensitive search for CH3D. We report rotational temperatures, production rates, and mixing ratios (with respect to H2O) at heliocentric distances of 1.49 AU (Nov. 2004) and 1.21 AU (Jan. 2005). We highlight three principal results: (1) The mixing ratios of parent volatiles measured at 1.49 AU and 1.21 AU agree within confidence limits, consistent with homogeneous composition in the mean volatile release from the nucleus of C/2004 Q2. Notably, the relative abundance of C2H6/C2H2 is substantially higher than those measured in other comets, while the mixing ratios C2H6/H2O, CH3OH/H2O, and HCN/H2O are similar to those observed in comets, referred to as "organics-normal". (2) The spin temperature of CH4 is > 35-38 K, an estimate consistent with the more robust spin temperature found for H2O. (3) We obtained a 3s upper limit of CH3D/CH4 < 0.020 (D/H < 0.005). This limit suggests that methane released from the nucleus of C/2004 Q2 is not dominated by a component formed in extremely cold (near 10 K) environments. Formation pathways of both interstellar and nebular origin consistent with the measured D/H in methane are discussed. Evaluating the relative contributions of these pathways requires further modeling of chemistry including both gas-phase and gas-grain processes in the natal interstellar cloud and in the protoplanetary disk. △ Less

Submitted 2 May, 2009; v1 submitted 21 April, 2009; originally announced April 2009.

Comments: Accepted by The Astrophysical Journal

Journal ref: Astrophys.J.699:1563-1572,2009

The unusual volatile composition of the Halley-type comet 8P/Tuttle: Addressing the existence of an Inner Oort Cloud

Authors: H. Boehnhardt, M. J. Mumma, G. L. Villanueva, M. A. DiSanti, B. P. Bonev, M. Lippi, H. U. Kaeufl

Abstract: We measured organic volatiles (CH4, CH3OH, C2H6, H2CO), CO, and water in comet 8P/Tuttle, a comet from the Oort cloud reservoir now in a short-period Halley-type orbit. We compare its composition with two other comets in Halley-type orbits, and with comets of the "organics-normal" and "organics-depleted" classes. Chemical gradients are expected in the comet-forming region of the proto-planetary… ▽ More We measured organic volatiles (CH4, CH3OH, C2H6, H2CO), CO, and water in comet 8P/Tuttle, a comet from the Oort cloud reservoir now in a short-period Halley-type orbit. We compare its composition with two other comets in Halley-type orbits, and with comets of the "organics-normal" and "organics-depleted" classes. Chemical gradients are expected in the comet-forming region of the proto-planetary disk, and an individual comet should reflect its specific heritage. If Halley-type comets came from the inner Oort cloud as proposed, we see no common characteristics that could distinguish such comets from those that were stored in the outer Oort cloud. △ Less

Submitted 24 July, 2008; originally announced July 2008.

Comments: 14 pages, including 1 figure and 2 Tables

Journal ref: Astrophys.J.683:L71-L74,2008

Discovery of OH in Circumstellar Disks Around Young Intermediate-Mass Stars

Authors: Avi M. Mandell, Michael J. Mumma, Goeffrey A. Blake, Boncho P. Bonev, Geronimo L. Villanueva, Colette Salyk

Abstract: We detect emission from multiple low-excitation ro-vibrational transitions of OH from the two Herbig Ae stars AB Aurigae and MWC 758 in the 3.0 - 3.7 micron wavelength range (L-band), using the NIRSPEC instrument on Keck II. The inner radius for the emitting region in both stars is close to 1 AU. We compare an optically thin LTE model and a thin-wedge fluorescence model, finding rotational tempe… ▽ More We detect emission from multiple low-excitation ro-vibrational transitions of OH from the two Herbig Ae stars AB Aurigae and MWC 758 in the 3.0 - 3.7 micron wavelength range (L-band), using the NIRSPEC instrument on Keck II. The inner radius for the emitting region in both stars is close to 1 AU. We compare an optically thin LTE model and a thin-wedge fluorescence model, finding rotational temperatures of 650 - 800 K and OH abundances of 10^42 - 10^45 molecules for the two stars. Comparisons with current chemical models support the fluorescence excitation model for AB Aurigae and possibly MWC 758, but further observations and detailed modeling are necessary to improve constraints on OH emission in different disk environments. △ Less

Submitted 22 May, 2008; originally announced May 2008.

Comments: 9 pages, 1 table, 3 figures, accepted to ApJL

The Peculiar Volatile Composition of Comet 8P/Tuttle: A Contact Binary of Chemically Distinct Cometesimals?

Authors: B. P. Bonev, M. J. Mumma, Y. L. Radeva, M. A. DiSanti, E. L. Gibb, G. L. Villanueva

Abstract: We report measurements of eight native (i.e., released directly from the comet nucleus) volatiles (H2O, HCN, CH4, C2H2, C2H6, CO, H2CO, and CH3OH) in comet 8P/Tuttle using NIRSPEC at Keck 2. Comet Tuttle reveals a truly unusual composition, distinct from that of any comet observed to date at infrared wavelengths. The prominent enrichment of methanol relative to water contrasts the depletions of… ▽ More We report measurements of eight native (i.e., released directly from the comet nucleus) volatiles (H2O, HCN, CH4, C2H2, C2H6, CO, H2CO, and CH3OH) in comet 8P/Tuttle using NIRSPEC at Keck 2. Comet Tuttle reveals a truly unusual composition, distinct from that of any comet observed to date at infrared wavelengths. The prominent enrichment of methanol relative to water contrasts the depletions of other molecules, especially C2H2 and HCN. We suggest that the nucleus of 8P/Tuttle may contain two cometesimals characterized by distinct volatile composition. The relative abundances C2/CN, C2/OH, and CN/OH in 8P/Tuttle (measured at optical/near-UV wavelengths) differ substantially from the mixing ratios of their potential parents (C2H2/HCN, C2H2/H2O, and HCN/H2O) found in this work. Based on this comparison, our results do not support C2H2 and HCN being the principal precursors for respectively C2 and CN in Tuttle. The peculiar native composition observed in 8P/Tuttle (compared to other comets) provides new strong evidence for chemical diversity in the volatile materials stored in comet nuclei. We discuss the implications of this diversity for expected variations in the deuterium enrichment of water among comets. △ Less

Submitted 29 April, 2008; originally announced April 2008.

Comments: Accepted for Astrophysical Journal Letters

IRCS/Subaru Observations of Water in the Inner Coma of Comet 73P-B/Schwassmann-Wachmann 3: Spatially Resolved Rotational Temperatures and Ortho-Para Ratios

Authors: Boncho P . Bonev, Michael J. Mumma, Hideyo Kawakita, Hitomi Kobayashi, Geronimo L. Villanueva

Abstract: Comet 73P-B/Schwassmann-Wachmann 3 was observed with IRCS/Subaru at geocentric distance of 0.074 AU on UT 10 May 2006. Multiple H2O emission lines were detected in non-resonant fluorescence near 2.9 microns. No significant variation in total H2O production rate was found during the (3 hour) duration of our observations. H2O rotational temperatures and ortho-to-para abundance ratios were measured… ▽ More Comet 73P-B/Schwassmann-Wachmann 3 was observed with IRCS/Subaru at geocentric distance of 0.074 AU on UT 10 May 2006. Multiple H2O emission lines were detected in non-resonant fluorescence near 2.9 microns. No significant variation in total H2O production rate was found during the (3 hour) duration of our observations. H2O rotational temperatures and ortho-to-para abundance ratios were measured for several positions in the coma. The temperatures extracted from two different time intervals show very similar spatial distributions. For both, the rotational temperature decreased from ~110 K to ~90 K as the projected distance from the nucleus increased from ~5 to ~30 km. We see no evidence for OPR change in the coma. The H2O ortho-para ratio is consistent with the statistical equilibrium value (3.0) for all spatially resolved measurements. This implies a nuclear spin temperature higher than ~ 45 K. △ Less

Submitted 28 March, 2008; originally announced March 2008.

Comments: accepted for Icarus