Showing 1–20 of 20 results for author: Lim, O

Promise and Peril: Stellar Contamination and Strict Limits on the Atmosphere Composition of TRAPPIST-1c from JWST NIRISS Transmission Spectra

Authors: Michael Radica, Caroline Piaulet-Ghorayeb, Jake Taylor, Louis-Philippe Coulombe, Loïc Albert, Étienne Artigau, Björn Benneke, Nicolas B. Cowan, René Doyon, David Lafrenière, Alexandrine L'Heureux, Olivia Lim

Abstract: Attempts to probe the atmospheres of rocky planets around M dwarfs present both promise and peril. While their favorable planet-to-star radius ratios enable searches for even thin secondary atmospheres, their high activity levels and high-energy outputs threaten atmosphere survival. Here, we present the 0.6--2.85$μ$m transmission spectrum of the 1.1 Earth-radius, ~340K rocky planet TRAPPIST-1c obt… ▽ More Attempts to probe the atmospheres of rocky planets around M dwarfs present both promise and peril. While their favorable planet-to-star radius ratios enable searches for even thin secondary atmospheres, their high activity levels and high-energy outputs threaten atmosphere survival. Here, we present the 0.6--2.85$μ$m transmission spectrum of the 1.1 Earth-radius, ~340K rocky planet TRAPPIST-1c obtained over two JWST NIRISS/SOSS transit observations. Each of the two spectra displays 100--500 ppm signatures of stellar contamination. Despite being separated by 367 days, the retrieved spot and faculae properties are consistent between the two visits, resulting in nearly identical transmission spectra. Jointly retrieving for stellar contamination and a planetary atmosphere rules out with high confidence (>3-$σ$) not only clear hydrogen-dominated atmospheres, but even thin, 1-bar high-mean molecular weight atmospheres rich in H$_2$O, NH$_3$, or CO (at the 2-$σ$ level). We find that the only atmosphere scenarios which our spectrum cannot rule out are CH$_4$- or CO$_2$-rich atmospheres, which are both unlikely to be retained when considering the photodestruction of CH$_4$ and the susceptibility of even a CO$_2$-rich atmosphere to escape given the cumulative high-energy irradiation experienced by the planet. Our results further stress the importance of robustly accounting for stellar contamination when analyzing JWST observations of exo-Earths around M dwarfs, as well as the need for high-fidelity stellar models to search for the potential signals of thin secondary atmospheres. △ Less

Submitted 1 October, 2024; v1 submitted 28 September, 2024; originally announced September 2024.

Comments: Minor typo corrections and figure axis adjustments compared to previous version

Updated forecast for TRAPPIST-1 times of transit for all seven exoplanets incorporating JWST data

Authors: Eric Agol, Natalie H. Allen, Björn Benneke, Laetitia Delrez, René Doyon, Elsa Ducrot, Néstor Espinoza, Amélie Gressier, David Lafrenière, Olivia Lim, Jacob Lustig-Yaeger, Caroline Piaulet-Ghorayeb, Michael Radica, Zafar Rustamkulov, Kristin S. Sotzen

Abstract: The TRAPPIST-1 system has been extensively observed with JWST in the near-infrared with the goal of measuring atmospheric transit transmission spectra of these temperate, Earth-sized exoplanets. A byproduct of these observations has been much more precise times of transit compared with prior available data from Spitzer, HST, or ground-based telescopes. In this note we use 23 new timing measurement… ▽ More The TRAPPIST-1 system has been extensively observed with JWST in the near-infrared with the goal of measuring atmospheric transit transmission spectra of these temperate, Earth-sized exoplanets. A byproduct of these observations has been much more precise times of transit compared with prior available data from Spitzer, HST, or ground-based telescopes. In this note we use 23 new timing measurements of all seven planets in the near-infrared from five JWST observing programs to better forecast and constrain the future times of transit in this system. In particular, we note that the transit times of TRAPPIST-1h have drifted significantly from a prior published analysis by up to tens of minutes. Our newer forecast has a higher precision, with median statistical uncertainties ranging from 7-105 seconds during JWST Cycles 4 and 5. Our expectation is that this forecast will help to improve planning of future observations of the TRAPPIST-1 planets, whereas we postpone a full dynamical analysis to future work. △ Less

Submitted 17 September, 2024; originally announced September 2024.

Comments: Submitted to AAS journals, 4 pages, 1 figure

Transmission Spectroscopy of the Habitable Zone Exoplanet LHS 1140 b with JWST/NIRISS

Authors: Charles Cadieux, René Doyon, Ryan J. MacDonald, Martin Turbet, Étienne Artigau, Olivia Lim, Michael Radica, Thomas J. Fauchez, Salma Salhi, Lisa Dang, Loïc Albert, Louis-Philippe Coulombe, Nicolas B. Cowan, David Lafrenière, Alexandrine L'Heureux, Caroline Piaulet, Björn Benneke, Ryan Cloutier, Benjamin Charnay, Neil J. Cook, Marylou Fournier-Tondreau, Mykhaylo Plotnykov, Diana Valencia

Abstract: LHS 1140 b is the second-closest temperate transiting planet to the Earth with an equilibrium temperature low enough to support surface liquid water. At 1.730$\pm$0.025 R$_\oplus$, LHS 1140 b falls within the radius valley separating H$_2$-rich mini-Neptunes from rocky super-Earths. Recent mass and radius revisions indicate a bulk density significantly lower than expected for an Earth-like rocky i… ▽ More LHS 1140 b is the second-closest temperate transiting planet to the Earth with an equilibrium temperature low enough to support surface liquid water. At 1.730$\pm$0.025 R$_\oplus$, LHS 1140 b falls within the radius valley separating H$_2$-rich mini-Neptunes from rocky super-Earths. Recent mass and radius revisions indicate a bulk density significantly lower than expected for an Earth-like rocky interior, suggesting that LHS 1140 b could either be a mini-Neptune with a small envelope of hydrogen ($\sim$0.1% by mass) or a water world (9--19% water by mass). Atmospheric characterization through transmission spectroscopy can readily discern between these two scenarios. Here, we present two JWST/NIRISS transit observations of LHS 1140 b, one of which captures a serendipitous transit of LHS 1140 c. The combined transmission spectrum of LHS 1140 b shows a telltale spectral signature of unocculted faculae (5.8 $σ$), covering $\sim$20% of the visible stellar surface. Besides faculae, our spectral retrieval analysis reveals tentative evidence of residual spectral features, best-fit by Rayleigh scattering from an N$_2$-dominated atmosphere (2.3 $σ$), irrespective of the consideration of atmospheric hazes. We also show through Global Climate Models (GCM) that H$_2$-rich atmospheres of various compositions (100$\times$, 300$\times$, 1000$\times$solar metallicity) are ruled out to $>$10 $σ$. The GCM calculations predict that water clouds form below the transit photosphere, limiting their impact on transmission data. Our observations suggest that LHS 1140 b is either airless or, more likely, surrounded by an atmosphere with a high mean molecular weight. Our tentative evidence of an N$_2$-rich atmosphere provides strong motivation for future transmission spectroscopy observations of LHS 1140 b. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 26 pages, 14 figures, 3 tables, accepted for publication in ApJL

NIRPS first light and early science: breaking the 1 m/s RV precision barrier at infrared wavelengths

Authors: Étienne Artigau, François Bouchy, René Doyon, Frédérique Baron, Lison Malo, François Wildi, Franceso Pepe, Neil J. Cook, Simon Thibault, Vladimir Reshetov, Xavier Dumusque, Christophe Lovis, Danuta Sosnowska, Bruno L. Canto Martins, Jose Renan De Medeiros, Xavier Delfosse, Nuno Santos, Rafael Rebolo, Manuel Abreu, Guillaume Allain, Romain Allart, Hugues Auger, Susana Barros, Luc Bazinet, Nicolas Blind , et al. (89 additional authors not shown)

Abstract: The Near-InfraRed Planet Searcher or NIRPS is a precision radial velocity spectrograph developed through collaborative efforts among laboratories in Switzerland, Canada, Brazil, France, Portugal and Spain. NIRPS extends to the 0.98-1.8 $μ$m domain of the pioneering HARPS instrument at the La Silla 3.6-m telescope in Chile and it has achieved unparalleled precision, measuring stellar radial velocit… ▽ More The Near-InfraRed Planet Searcher or NIRPS is a precision radial velocity spectrograph developed through collaborative efforts among laboratories in Switzerland, Canada, Brazil, France, Portugal and Spain. NIRPS extends to the 0.98-1.8 $μ$m domain of the pioneering HARPS instrument at the La Silla 3.6-m telescope in Chile and it has achieved unparalleled precision, measuring stellar radial velocities in the infrared with accuracy better than 1 m/s. NIRPS can be used either stand-alone or simultaneously with HARPS. Commissioned in late 2022 and early 2023, NIRPS embarked on a 5-year Guaranteed Time Observation (GTO) program in April 2023, spanning 720 observing nights. This program focuses on planetary systems around M dwarfs, encompassing both the immediate solar vicinity and transit follow-ups, alongside transit and emission spectroscopy observations. We highlight NIRPS's current performances and the insights gained during its deployment at the telescope. The lessons learned and successes achieved contribute to the ongoing advancement of precision radial velocity measurements and high spectral fidelity, further solidifying NIRPS' role in the forefront of the field of exoplanets. △ Less

Submitted 13 June, 2024; v1 submitted 12 June, 2024; originally announced June 2024.

Comments: Proceeding at the SPIE Astronomical Telescopes + Instrumentation conference [Yokohama,Japan; June 2024]

A roadmap for the atmospheric characterization of terrestrial exoplanets with JWST

Authors: TRAPPIST-1 JWST Community Initiative, :, Julien de Wit, René Doyon, Benjamin V. Rackham, Olivia Lim, Elsa Ducrot, Laura Kreidberg, Björn Benneke, Ignasi Ribas, David Berardo, Prajwal Niraula, Aishwarya Iyer, Alexander Shapiro, Nadiia Kostogryz, Veronika Witzke, Michaël Gillon, Eric Agol, Victoria Meadows, Adam J. Burgasser, James E. Owen, Jonathan J. Fortney, Franck Selsis, Aaron Bello-Arufe, Zoë de Beurs , et al. (58 additional authors not shown)

Abstract: Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets. While JWST Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a bet… ▽ More Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets. While JWST Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a better understanding of their host star. Here, we propose a roadmap to characterize the TRAPPIST-1 system -- and others like it -- in an efficient and robust manner. We notably recommend that -- although more challenging to schedule -- multi-transit windows be prioritized to mitigate the effects of stellar activity and gather up to twice more transits per JWST hour spent. We conclude that, for such systems, planets cannot be studied in isolation by small programs, but rather need large-scale, jointly space- and ground-based initiatives to fully exploit the capabilities of JWST for the exploration of terrestrial planets. △ Less

Submitted 22 July, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

Journal ref: Nature Astronomy (2024) 8, 810-818

Near-Infrared Transmission Spectroscopy of HAT-P-18$\,$b with NIRISS: Disentangling Planetary and Stellar Features in the Era of JWST

Authors: Marylou Fournier-Tondreau, Ryan J. MacDonald, Michael Radica, David Lafrenière, Luis Welbanks, Caroline Piaulet, Louis-Philippe Coulombe, Romain Allart, Kim Morel, Étienne Artigau, Loïc Albert, Olivia Lim, René Doyon, Björn Benneke, Jason F. Rowe, Antoine Darveau-Bernier, Nicolas B. Cowan, Nikole K. Lewis, Neil James Cook, Laura Flagg, Frédéric Genest, Stefan Pelletier, Doug Johnstone, Lisa Dang, Lisa Kaltenegger , et al. (2 additional authors not shown)

Abstract: The JWST Early Release Observations (ERO) included a NIRISS/SOSS (0.6-2.8$\,μ$m) transit of the $\sim\,$850$\,$K Saturn-mass exoplanet HAT-P-18$\,$b. Initial analysis of these data reported detections of water, escaping helium, and haze. However, active K dwarfs like HAT-P-18 possess surface heterogeneities $-$ starspots and faculae $-$ that can complicate the interpretation of transmission spectr… ▽ More The JWST Early Release Observations (ERO) included a NIRISS/SOSS (0.6-2.8$\,μ$m) transit of the $\sim\,$850$\,$K Saturn-mass exoplanet HAT-P-18$\,$b. Initial analysis of these data reported detections of water, escaping helium, and haze. However, active K dwarfs like HAT-P-18 possess surface heterogeneities $-$ starspots and faculae $-$ that can complicate the interpretation of transmission spectra, and indeed, a spot-crossing event is present in HAT-P-18$\,$b's NIRISS/SOSS light curves. Here, we present an extensive reanalysis and interpretation of the JWST ERO transmission spectrum of HAT-P-18$\,$b, as well as HST/WFC3 and $\textit{Spitzer}$/IRAC transit observations. We detect H$_2$O (12.5$\,σ$), CO$_2$ (7.3$\,σ$), a cloud deck (7.4$\,σ$), and unocculted starspots (5.8$\,σ$), alongside hints of Na (2.7$\,σ$). We do not detect the previously reported CH$_4$ ($\log$ CH$_4$ $<$ -6 to 2$\,σ$). We obtain excellent agreement between three independent retrieval codes, which find a sub-solar H$_2$O abundance ($\log$ H$_2$O $\approx -4.4 \pm 0.3$). However, the inferred CO$_2$ abundance ($\log$ CO$_2$ $\approx -4.8 \pm 0.4$) is significantly super-solar and requires further investigation into its origin. We also introduce new stellar heterogeneity considerations by fitting for the active regions' surface gravities $-$ a proxy for the effects of magnetic pressure. Finally, we compare our JWST inferences to those from HST/WFC3 and $\textit{Spitzer}$/IRAC. Our results highlight the exceptional promise of simultaneous planetary atmosphere and stellar heterogeneity constraints in the era of JWST and demonstrate that JWST transmission spectra may warrant more complex treatments of the transit light source effect. △ Less

Submitted 18 December, 2023; v1 submitted 23 October, 2023; originally announced October 2023.

Characterizing the Near-infrared Spectra of Flares from TRAPPIST-1 During JWST Transit Spectroscopy Observations

Authors: Ward S. Howard, Adam F. Kowalski, Laura Flagg, Meredith A. MacGregor, Olivia Lim, Michael Radica, Caroline Piaulet, Pierre-Alexis Roy, David Lafrenière, Björn Benneke, Alexander Brown, Néstor Espinoza, René Doyon, Louis-Philippe Coulombe, Doug Johnstone, Nicolas B. Cowan, Ray Jayawardhana, Jake D. Turner, Lisa Dang

Abstract: We present the first analysis of JWST near-infrared spectroscopy of stellar flares from TRAPPIST-1 during transits of rocky exoplanets. Four flares were observed from 0.6--2.8 $μ$m with NIRISS and 0.6--3.5 $μ$m with NIRSpec during transits of TRAPPIST-1b, f, and g. We discover P$α$ and Br$β$ line emission and characterize flare continuum at wavelengths from 1--3.5 $μ$m for the first time. Observed… ▽ More We present the first analysis of JWST near-infrared spectroscopy of stellar flares from TRAPPIST-1 during transits of rocky exoplanets. Four flares were observed from 0.6--2.8 $μ$m with NIRISS and 0.6--3.5 $μ$m with NIRSpec during transits of TRAPPIST-1b, f, and g. We discover P$α$ and Br$β$ line emission and characterize flare continuum at wavelengths from 1--3.5 $μ$m for the first time. Observed lines include H$α$, P$α$-P$ε$, Br$β$, He I $λ$0.7062$μ$m, two Ca II infrared triplet (IRT) lines, and the He I IRT. We observe a reversed Paschen decrement from P$α$-P$γ$ alongside changes in the light curve shapes of these lines. The continuum of all four flares is well-described by blackbody emission with an effective temperature below 5300 K, lower than temperatures typically observed at optical wavelengths. The 0.6--1 $μ$m spectra were convolved with the TESS response, enabling us to measure the flare rate of TRAPPIST-1 in the TESS bandpass. We find flares of 10$^{30}$ erg large enough to impact transit spectra occur at a rate of 3.6$\substack{+2.1 \\ -1.3}$ flare d$^{-1}$, $\sim$10$\times$ higher than previous predictions from K2. We measure the amount of flare contamination at 2 $μ$m for the TRAPPIST-1b and f transits to be 500$\pm$450 and 2100$\pm$400 ppm, respectively. We find up to 80% of flare contamination can be removed, with mitigation most effective from 1.0--2.4 $μ$m. These results suggest transits affected by flares may still be useful for atmospheric characterization efforts. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 29 pages, 17 figures, 3 tables, accepted to The Astrophysical Journal

Atmospheric Reconnaissance of TRAPPIST-1 b with JWST/NIRISS: Evidence for Strong Stellar Contamination in the Transmission Spectra

Authors: Olivia Lim, Björn Benneke, René Doyon, Ryan J. MacDonald, Caroline Piaulet, Étienne Artigau, Louis-Philippe Coulombe, Michael Radica, Alexandrine L'Heureux, Loïc Albert, Benjamin V. Rackham, Julien de Wit, Salma Salhi, Pierre-Alexis Roy, Laura Flagg, Marylou Fournier-Tondreau, Jake Taylor, Neil J. Cook, David Lafrenière, Nicolas B. Cowan, Lisa Kaltenegger, Jason F. Rowe, Néstor Espinoza, Lisa Dang, Antoine Darveau-Bernier

Abstract: TRAPPIST-1 is a nearby system of seven Earth-sized, temperate, rocky exoplanets transiting a Jupiter-sized M8.5V star, ideally suited for in-depth atmospheric studies. Each TRAPPIST-1 planet has been observed in transmission both from space and from the ground, confidently rejecting cloud-free, hydrogen-rich atmospheres. Secondary eclipse observations of TRAPPIST-1 b with JWST/MIRI are consistent… ▽ More TRAPPIST-1 is a nearby system of seven Earth-sized, temperate, rocky exoplanets transiting a Jupiter-sized M8.5V star, ideally suited for in-depth atmospheric studies. Each TRAPPIST-1 planet has been observed in transmission both from space and from the ground, confidently rejecting cloud-free, hydrogen-rich atmospheres. Secondary eclipse observations of TRAPPIST-1 b with JWST/MIRI are consistent with little to no atmosphere given the lack of heat redistribution. Here we present the first transmission spectra of TRAPPIST-1 b obtained with JWST/NIRISS over two visits. The two transmission spectra show moderate to strong evidence of contamination from unocculted stellar heterogeneities, which dominates the signal in both visits. The transmission spectrum of the first visit is consistent with unocculted starspots and the second visit exhibits signatures of unocculted faculae. Fitting the stellar contamination and planetary atmosphere either sequentially or simultaneously, we confirm the absence of cloud-free hydrogen-rich atmospheres, but cannot assess the presence of secondary atmospheres. We find that the uncertainties associated with the lack of stellar model fidelity are one order of magnitude above the observation precision of 89 ppm (combining the two visits). Without affecting the conclusion regarding the atmosphere of TRAPPIST-1 b, this highlights an important caveat for future explorations, which calls for additional observations to characterize stellar heterogeneities empirically and/or theoretical works to improve model fidelity for such cool stars. This need is all the more justified as stellar contamination can affect the search for atmospheres around the outer, cooler TRAPPIST-1 planets for which transmission spectroscopy is currently the most efficient technique. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Comments: 26 pages, 11 figures, accepted for publication in The Astrophysical Journal Letters

Vanadium oxide and a sharp onset of cold-trapping on a giant exoplanet

Authors: Stefan Pelletier, Björn Benneke, Mohamad Ali-Dib, Bibiana Prinoth, David Kasper, Andreas Seifahrt, Jacob L. Bean, Florian Debras, Baptiste Klein, Luc Bazinet, H. Jens Hoeijmakers, Aurora Y. Kesseli, Olivia Lim, Andres Carmona, Lorenzo Pino, Núria Casasayas-Barris, Thea Hood, Julian Stürmer

Abstract: The abundance of refractory elements in giant planets can provide key insights into their formation histories. Due to the Solar System giants' low temperatures, refractory elements condense below the cloud deck limiting sensing capabilities to only highly volatile elements. Recently, ultra-hot giant exoplanets have allowed for some refractory elements to be measured showing abundances broadly cons… ▽ More The abundance of refractory elements in giant planets can provide key insights into their formation histories. Due to the Solar System giants' low temperatures, refractory elements condense below the cloud deck limiting sensing capabilities to only highly volatile elements. Recently, ultra-hot giant exoplanets have allowed for some refractory elements to be measured showing abundances broadly consistent with the solar nebula with titanium likely condensed out of the photosphere. Here we report precise abundance constraints of 14 major refractory elements on the ultra-hot giant planet WASP-76b that show distinct deviations from proto-solar, and a sharp onset in condensation temperature. In particular, we find nickel to be enriched, a possible sign of the accretion of a differentiated object's core during the planet's evolution. Elements with condensation temperatures below 1,550 K otherwise closely match those of the Sun before sharply transitioning to being strongly depleted above 1,550 K, well explained by nightside cold-trapping. We further unambiguously detect vanadium oxide on WASP-76b, a molecule long hypothesized to drive atmospheric thermal inversions, and also observe a global east-west asymmetry in its absorption signals. Overall, our findings indicate that giant planets have a mostly stellar-like refractory elemental content and suggest that temperature sequences of hot Jupiter spectra can show abrupt transitions wherein a mineral species is either present, or completely absent if a cold-trap exists below its condensation temperature. △ Less

Submitted 14 June, 2023; originally announced June 2023.

Comments: Published online in Nature on June 14, 2023

ATOCA: an algorithm to treat order contamination. Application to the NIRISS SOSS mode

Authors: Antoine Darveau-Bernier, Loïc Albert, Geert Jan Talens, David Lafrenière, Michael Radica, René Doyon, Neil J. Cook, Jason F. Rowe, Étienne Artigau, Björn Benneke, Nicolas Cowan, Lisa Dang, Néstor Espinoza, Doug Johnstone, Lisa Kaltenegger, Olivia Lim, Stefan Pelletier, Caroline Piaulet, Arpita Roy, Pierre-Alexis Roy, Jared Splinter, Jake Taylor, Jake D. Turner

Abstract: After a successful launch, the James Webb Space Telescope is preparing to undertake one of its principal missions, the characterization of the atmospheres of exoplanets. The Single Object Slitless Spectroscopy (SOSS) mode of the Near Infrared Imager and Slitless Spectrograph (NIRISS) is the only observing mode that has been specifically designed for this objective. It features a wide simultaneous… ▽ More After a successful launch, the James Webb Space Telescope is preparing to undertake one of its principal missions, the characterization of the atmospheres of exoplanets. The Single Object Slitless Spectroscopy (SOSS) mode of the Near Infrared Imager and Slitless Spectrograph (NIRISS) is the only observing mode that has been specifically designed for this objective. It features a wide simultaneous spectral range (0.6--2.8\,\micron) through two spectral diffraction orders. However, due to mechanical constraints, these two orders overlap slightly over a short range, potentially introducing a ``contamination'' signal in the extracted spectrum. We show that for a typical box extraction, this contaminating signal amounts to 1\% or less over the 1.6--2.8\,\micron\ range (order 1), and up to 1\% over the 0.85--0.95\,\micron\ range (order 2). For observations of exoplanet atmospheres (transits, eclipses or phase curves) where only temporal variations in flux matter, the contamination signal typically biases the results by order of 1\% of the planetary atmosphere spectral features strength. To address this problem, we developed the Algorithm to Treat Order ContAmination (ATOCA). By constructing a linear model of each pixel on the detector, treating the underlying incident spectrum as a free variable, ATOCA is able to perform a simultaneous extraction of both orders. We show that, given appropriate estimates of the spatial trace profiles, the throughputs, the wavelength solutions, as well as the spectral resolution kernels for each order, it is possible to obtain an extracted spectrum accurate to within 10\,ppm over the full spectral range. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: Submitted to PASP. 22 pages, 12 figures

A Review on MR Based Human Brain Parcellation Methods

Authors: Pantea Moghimi, Anh The Dang, Theoden I. Netoff, Kelvin O. Lim, Gowtham Atluri

Abstract: Brain parcellations play a ubiquitous role in the analysis of magnetic resonance imaging (MRI) datasets. Over 100 years of research has been conducted in pursuit of an ideal brain parcellation. Different methods have been developed and studied for constructing brain parcellations using different imaging modalities. More recently, several data-driven parcellation methods have been adopted from data… ▽ More Brain parcellations play a ubiquitous role in the analysis of magnetic resonance imaging (MRI) datasets. Over 100 years of research has been conducted in pursuit of an ideal brain parcellation. Different methods have been developed and studied for constructing brain parcellations using different imaging modalities. More recently, several data-driven parcellation methods have been adopted from data mining, machine learning, and statistics communities. With contributions from different scientific fields, there is a rich body of literature that needs to be examined to appreciate the breadth of existing research and the gaps that need to be investigated. In this work, we review the large body of in vivo brain parcellation research spanning different neuroimaging modalities and methods. A key contribution of this work is a semantic organization of this large body of work into different taxonomies, making it easy to understand the breadth and depth of the brain parcellation literature. Specifically, we categorized the existing parcellations into three groups: Anatomical parcellations, functional parcellations, and structural parcellations which are constructed using T1-weighted MRI, functional MRI (fMRI), and diffusion-weighted imaging (DWI) datasets, respectively. We provide a multi-level taxonomy of different methods studied in each of these categories, compare their relative strengths and weaknesses, and highlight the challenges currently faced for the development of brain parcellations. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: 31 pages, 3 figures, 2 tables

Interoception Underlies The Therapeutic Effects of Mindfulness Meditation for Post-Traumatic Stress Disorder: A Randomized Clinical Trial

Authors: Seung Suk Kang, Ph. D., Scott R. Sponheim, Ph. D., Kelvin O. Lim, M. D

Abstract: Mindfulness-based interventions have proven its efficacy in treating post-traumatic stress disorder (PTSD), but the underlying neurobiological mechanism is unknown. To determine the neurobiological mechanism of action of mindfulness-based stress reduction (MBSR) treating PTSD, we conducted a randomized clinical trial (RCT) in which 98 veterans with PTSD were randomly assigned to receive MBSR thera… ▽ More Mindfulness-based interventions have proven its efficacy in treating post-traumatic stress disorder (PTSD), but the underlying neurobiological mechanism is unknown. To determine the neurobiological mechanism of action of mindfulness-based stress reduction (MBSR) treating PTSD, we conducted a randomized clinical trial (RCT) in which 98 veterans with PTSD were randomly assigned to receive MBSR therapy (n = 47) or present-centered group therapy (PCGT; n = 51; an active-control condition). Pre- and post-intervention measures of PTSD symptom severity (PTSD Checklist) and brain activity measures of electroencephalography (EEG) were assessed, including spectral power of spontaneous neural oscillatory activities during resting and meditation periods, time-frequency (TF) power of cognitive task-related brain responses, and TF power of heartbeat-evoked brain responses (HEBR) that reflect cardiac interoceptive brain responses during resting and meditation. Compared to controls, the MBSR group had greater improvements in PTSD symptoms, spontaneous EEG alpha (8-13 Hz) power in posterior sites, task-related frontal theta power (4-7 Hz in 140-220 ms post-stimulus), and frontal theta HEBR (3-5 Hz and 265-328 ms post-R-peak). Latent difference score modeling found that only the changes in the frontal theta HEBR mediated the MBSR treatment effect. Brain source-level analysis found that the theta HEBR changes in the anterior cingulate cortex, anterior insular cortex, and the lateral prefrontal cortex predicted PTSD symptom improvements. These results indicated that mindfulness meditation improves spontaneous brain activities reflecting internally oriented relaxation and brain functions of attentional control. However, interoceptive brain capacity enhanced by MBSR appears to be the primary cerebral mechanism that regulates emotional disturbances and improves anxiety symptoms of PTSD. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 17 pages, 5 figures, 1 table

Construction and Evaluation of Hierarchical Parcellation of the Brain using fMRI with Prewhitening

Authors: Pantea Moghimi, Kelvin O. Lim, Theoden I. Netoff

Abstract: Brain atlases are a ubiquitous tool used for analyzing and interpreting brain imaging datasets. Traditionally, brain atlases divided the brain into regions separated by anatomical landmarks. In the last decade, several attempts have been made to parcellate the brain into regions with distinct functional activity using fMRI. To construct a brain atlas using fMRI, data driven algorithms are used to… ▽ More Brain atlases are a ubiquitous tool used for analyzing and interpreting brain imaging datasets. Traditionally, brain atlases divided the brain into regions separated by anatomical landmarks. In the last decade, several attempts have been made to parcellate the brain into regions with distinct functional activity using fMRI. To construct a brain atlas using fMRI, data driven algorithms are used to group voxels with similar functional activity together to form regions. Hierarchical clustering is one parcellation method that has been used for functional parcellation of the brain, resulting in parcellations that align well with cytoarchitectonic divisions of the brain. However, few rigorous data driven evaluations of the method have been performed. Moreover, the effect of removing autocorrelation trends from fMRI time series (prewhitening) on the structure of the resultant atlas has not been previously explored. In this paper, we use hierarchical clustering to produce functional parcellations of the brain using hierarchical clustering. We use both prewhitened and raw fMRI time series to construct the atlas. The resultant atlases were then evaluated for their homogeneity, separation between regions, reproducibility across subjects, and reproducibility across scans. △ Less

Submitted 7 February, 2018; v1 submitted 21 December, 2017; originally announced December 2017.

Comments: 16 pages, 7 figures

SPIRou Input Catalogue: Global properties of 440 M dwarfs observed with ESPaDOnS at CFHT

Authors: Pascal Fouqué, Claire Moutou, Lison Malo, Eder Martioli, Olivia Lim, Arvind Rajpurohit, Etienne Artigau, Xavier Delfosse, Jean-François Donati, Thierry Forveille, Julien Morin, France Allard, Raphaël Delage, René Doyon, Elodie Hébrard, Vasco Neves

Abstract: Present and future high-precision radial-velocity spectrometers dedicated to the discovery of low-mass planets orbiting low-mass dwarfs need to focus on the best selected stars to make an efficient use of telescope time. In the framework of the preparation of the SPIRou Input Catalog, the CoolSnap program aims at screening M dwarfs in the solar neighborhood against binarity, rapid rotation, activi… ▽ More Present and future high-precision radial-velocity spectrometers dedicated to the discovery of low-mass planets orbiting low-mass dwarfs need to focus on the best selected stars to make an efficient use of telescope time. In the framework of the preparation of the SPIRou Input Catalog, the CoolSnap program aims at screening M dwarfs in the solar neighborhood against binarity, rapid rotation, activity, ... To optimize the selection, the present paper describes the methods used to compute effective temperature, metallicity, projected rotation velocity of a large sample of 440 M dwarfs observed in the visible with the high-resolution spectro-polarimeter ESPaDOnS at CFHT. It also summarizes known and newly-discovered spectroscopic binaries, and stars known to belong to visual multiple systems. A calibration of the projected rotation velocity versus measured line widths for M dwarfs observed by the ESPaDOnS spectro-polarimeter is derived, and the resulting values are compared to equatorial rotation velocities deduced from rotation periods and radii. A comparison of the derived effective temperatures and metallicities with literature values is also conducted. Finally, the radial velocity uncertainty of each star in the sample is estimated, to narrow down the selection of stars to be included into the SPIRou Input Catalogue (SPIC). △ Less

Submitted 12 December, 2017; originally announced December 2017.

Comments: 34 pages, 12 figures, 13 tables, MNRAS in press December 11, 2017

Choosing Wavelet Methods, Filters, and Lengths for Functional Brain Network Construction

Authors: Zitong Zhang, Qawi K. Telesford, Chad Giusti, Kelvin O. Lim, Danielle S. Bassett

Abstract: Wavelet methods are widely used to decompose fMRI, EEG, or MEG signals into time series representing neurophysiological activity in fixed frequency bands. Using these time series, one can estimate frequency-band specific functional connectivity between sensors or regions of interest, and thereby construct functional brain networks that can be examined from a graph theoretic perspective. Despite th… ▽ More Wavelet methods are widely used to decompose fMRI, EEG, or MEG signals into time series representing neurophysiological activity in fixed frequency bands. Using these time series, one can estimate frequency-band specific functional connectivity between sensors or regions of interest, and thereby construct functional brain networks that can be examined from a graph theoretic perspective. Despite their common use, however, practical guidelines for the choice of wavelet method, filter, and length have remained largely undelineated. Here, we explicitly explore the effects of wavelet method (MODWT vs. DWT), wavelet filter (Daubechies Extremal Phase, Daubechies Least Asymmetric, and Coiflet families), and wavelet length (2 to 24) - each essential parameters in wavelet-based methods - on the estimated values of network diagnostics and in their sensitivity to alterations in psychiatric disease. We observe that the MODWT method produces less variable estimates than the DWT method. We also observe that the length of the wavelet filter chosen has a greater impact on the estimated values of network diagnostics than the type of wavelet chosen. Furthermore, wavelet length impacts the sensitivity of the method to detect differences between health and disease and tunes classification accuracy. Collectively, our results suggest that the choice of wavelet method and length significantly alters the reliability and sensitivity of these methods in estimating values of network diagnostics drawn from graph theory. They furthermore demonstrate the importance of reporting the choices utilized in neuroimaging studies and support the utility of exploring wavelet parameters to maximize classification accuracy in the development of biomarkers of psychiatric disease and neurological disorders. △ Less

Submitted 23 December, 2015; originally announced December 2015.

Comments: working paper

Resolving Structure in Human Brain Organization: Identifying Mesoscale Organization in Weighted Network Representations

Authors: Christian Lohse, Danielle S. Bassett, Kelvin O. Lim, Jean M. Carlson

Abstract: Human brain anatomy and function display a combination of modular and hierarchical organization, suggesting the importance of both cohesive structures and variable resolutions in the facilitation of healthy cognitive processes. However, tools to simultaneously probe these features of brain architecture require further development. We propose and apply a set of methods to extract cohesive structure… ▽ More Human brain anatomy and function display a combination of modular and hierarchical organization, suggesting the importance of both cohesive structures and variable resolutions in the facilitation of healthy cognitive processes. However, tools to simultaneously probe these features of brain architecture require further development. We propose and apply a set of methods to extract cohesive structures in network representations of brain connectivity using multi-resolution techniques. We employ a combination of soft thresholding, windowed thresholding, and resolution in community detection, that enable us to identify and isolate structures associated with different weights. One such mesoscale structure is bipartivity, which quantifies the extent to which the brain is divided into two partitions with high connectivity between partitions and low connectivity within partitions. A second, complementary mesoscale structure is modularity, which quantifies the extent to which the brain is divided into multiple communities with strong connectivity within each community and weak connectivity between communities. Our methods lead to multi-resolution curves of these network diagnostics over a range of spatial, geometric, and structural scales. For statistical comparison, we contrast our results with those obtained for several benchmark null models. Our work demonstrates that multi-resolution diagnostic curves capture complex organizational profiles in weighted graphs. We apply these methods to the identification of resolution-specific characteristics of healthy weighted graph architecture and altered connectivity profiles in psychiatric disease. △ Less

Submitted 20 December, 2013; originally announced December 2013.

Comments: Comments welcome

Optimization of gain in traveling-wave optical parametric amplifiers by tuning the offset between pump- and signal-waist locations

Authors: Gideon Alon, Oo-Kaw Lim, Amar Bhagwat, Chao-Hsiang Chen, Muthiah Annamalai, Michael Vasilyev, Prem Kumar

Abstract: We present experimental demonstration and modeling of the optimization of a phase-sensitive optical parametric amplifier by tuning the relative position between the pump- and signal-beam waists along the propagation direction. At the optimum position, the pump beam focuses after the signal beam, and this departure from co-located waists increases with increasing pump power. Such optimization leads… ▽ More We present experimental demonstration and modeling of the optimization of a phase-sensitive optical parametric amplifier by tuning the relative position between the pump- and signal-beam waists along the propagation direction. At the optimum position, the pump beam focuses after the signal beam, and this departure from co-located waists increases with increasing pump power. Such optimization leads to more than 3\,dB improvement in the measured de-amplification response of the amplifier. △ Less

Submitted 1 March, 2013; v1 submitted 4 October, 2012; originally announced October 2012.

Comments: 4 pages, 4 figures

Observation of twin beam correlations and quadrature entanglement by frequency doubling in a two-port resonator

Authors: O. -K. Lim, B. Boland, M. Saffman

Abstract: We demonstrate production of quantum correlated and entangled beams by second harmonic generation in a nonlinear resonator with two output ports. The output beams at wavelength 428.5 nm exhibit 0.9 dB of nonclassical intensity correlations and 0.3 dB of entanglement. We demonstrate production of quantum correlated and entangled beams by second harmonic generation in a nonlinear resonator with two output ports. The output beams at wavelength 428.5 nm exhibit 0.9 dB of nonclassical intensity correlations and 0.3 dB of entanglement. △ Less

Submitted 4 September, 2006; originally announced September 2006.

Comments: 5 pages, 7 figures

Journal ref: Europhys. Lett. 78, 40004 (2007)

Intensity correlations and entanglement by frequency doubling in a dual ported resonator

Authors: Oo-Kaw Lim, Mark Saffman

Abstract: We show that non-classical intensity correlations and quadrature entanglement can be generated by frequency doubling in a resonator with two output ports. We predict twin-beam intensity correlations 6 dB below the coherent state limit, and that the product of the inference variances of the quadrature fluctuations gives an Einstein-Podolsky-Rosen (EPR) correlation coefficient of V_EPR = 0.6 < 1.… ▽ More We show that non-classical intensity correlations and quadrature entanglement can be generated by frequency doubling in a resonator with two output ports. We predict twin-beam intensity correlations 6 dB below the coherent state limit, and that the product of the inference variances of the quadrature fluctuations gives an Einstein-Podolsky-Rosen (EPR) correlation coefficient of V_EPR = 0.6 < 1. Comparison with an entanglement source based on combining two frequency doublers with a beam splitter shows that the dual ported resonator provides stronger entanglement at lower levels of individual beam squeezing. Calculations are performed using a self-consistent propagation method that does not invoke a mean field approximation. Results are given for physically realistic parameters that account for the Gaussian shape of the intracavity beams, as well as intracavity losses. △ Less

Submitted 26 March, 2006; v1 submitted 9 March, 2006; originally announced March 2006.

Comments: 12 pages, 9 figures, normalization corrected, fig 3 and fig 7 updated

Journal ref: Phys. Rev. A 74, 023816 (2006)

Creation, doubling, and splitting, of vortices in intracavity second harmonic generation

Authors: O-K Lim, B Boland, M Saffman, W. Krolikowski

Abstract: We demonstrate generation and frequency doubling of unit charge vortices in a linear astigmatic resonator. Topological instability of the double charge harmonic vortices leads to well separated vortex cores that are shown to rotate, and become anisotropic, as the resonator is tuned across resonance. We demonstrate generation and frequency doubling of unit charge vortices in a linear astigmatic resonator. Topological instability of the double charge harmonic vortices leads to well separated vortex cores that are shown to rotate, and become anisotropic, as the resonator is tuned across resonance. △ Less

Submitted 6 October, 2003; originally announced October 2003.

Journal ref: J. Opt. A: Pure Appl. Opt. 6, 486-489 (2004)