Catalog footprint

What is connected

74works
21topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

74 published item(s)

preprint2026arXiv

Multi-Agent Collaborative Reward Design for Enhancing Reasoning in Reinforcement Learning

We present CRM (Multi-Agent Collaborative Reward Model), a framework that replaces a single black-box reward model with a coordinated team of specialist evaluators to improve robustness and interpretability in RLHF. Conventional reward models struggle to jointly optimize multiple, sometimes conflicting, preference dimensions (e.g., factuality, helpfulness, safety) and offer limited transparency into why a score is assigned. CRM addresses these issues by decomposing preference evaluation into domain-specific agents that each produce partial signals, alongside global evaluators such as ranker-based and embedding-similarity rewards. A centralized aggregator fuses these signals at each timestep, balancing factors like step-wise correctness, multi-agent agreement, and repetition penalties, yielding a single training reward compatible with standard RL pipelines. The policy is optimized with advantage-based updates (e.g., GAE), while a value model regresses to the aggregated reward, enabling multi-perspective reward shaping without requiring additional human annotations beyond those used to train the evaluators. To support training and assessment, we introduce rewardBench, a benchmark and training suite aligned with the collaborative structure of CRM. Together, CRM and rewardBench provide a practical, modular path to more transparent reward modeling and more stable optimization.

preprint2026arXiv

S$^2$F: Principled Hybrid Testing With Fuzzing, Symbolic Execution, and Sampling

Hybrid testing that integrates fuzzing, symbolic execution, and sampling has demonstrated superior testing efficiency compared to individual techniques. However, the state-of-the-art (SOTA) hybrid testing tools do not fully exploit the capabilities of symbolic execution and sampling in two key aspects. First, the SOTA hybrid testing tools employ tailored symbolic execution engines that tend to over-prune branches, leading to considerable time wasted waiting for seeds from the fuzzer and missing opportunities to discover crashes. Second, existing methods do not apply sampling to the appropriate branches and therefore cannot utilize the full capability of sampling. To address these two limitations, we propose a novel hybrid testing architecture that combines the precision of conventional symbolic execution with the scalability of tailored symbolic execution engines. Based on this architecture, we propose several principles for combining fuzzing, symbolic execution, and sampling. We implement our method in a hybrid testing tool S$^2$F. To evaluate its effectiveness, we conduct extensive experiments on 15 real-world programs. Experimental results demonstrate that S$^2$F outperforms the SOTA tool, achieving an average improvement of 6.14% in edge coverage and 32.6% in discovered crashes. Notably, our tool uncovers three previously unknown crashes in real-world programs.

preprint2025arXiv

CUT: Pruning Pre-Trained Multi-Task Models into Compact Models for Edge Devices

Multi-task learning has garnered widespread attention in the industry due to its efficient data utilization and strong generalization capabilities, making it particularly suitable for providing high-quality intelligent services to users. Edge devices, as the primary platforms directly serving users, play a crucial role in delivering multi-task services. However, current multi-task models are often large, and user task demands are increasingly diverse. Deploying such models directly on edge devices not only increases the burden on these devices but also leads to task redundancy. To address this issue, this paper innovatively proposes a pre-trained multi-task model pruning method specifically designed for edge computing. The goal is to utilize existing pre-trained multi-task models to construct a compact multi-task model that meets the needs of edge devices. The specific implementation steps are as follows: First, decompose the tasks within the pre-trained multi-task model and select tasks based on actual user needs. Next, while retaining the knowledge of the original pre-trained model, evaluate parameter importance and use a parameter fusion method to effectively integrate shared parameters among tasks. Finally, obtain a compact multi-task model suitable for edge devices. To validate the effectiveness of the proposed method, we conducted experiments on three public image datasets. The experimental results fully demonstrate the superiority and efficiency of this method, providing a new solution for multi-task learning on edge devices.

preprint2025arXiv

Dynamical Architectures of S-type Transiting Planets in Binaries II: A Dichotomy in Orbital Alignment of Small Planets in Close Binary Systems

Stellar multiplicity plays a crucial role in shaping planet formation and dynamical evolution. We present a survey of 54 TESS Objects of Interest (TOIs) within 300 pc that exhibit significant Hipparcos-Gaia astrometric accelerations. We identified 35 TOIs with stellar companions at projected separations between $0.1^{\prime\prime}$ to $2^{\prime\prime}$ (or $10-200$ AU). We also identified 12 TOIs that could host planetary-mass or brown dwarf companions, including 6 that are newly discovered. Furthermore, we perform three-dimensional orbital characterization for 12 binaries hosting confirmed planets or planet candidates, allowing us to constrain the line-of-sight mutual inclination, $ΔI_{\mathrm{los}}$, between the planetary and binary orbits. Combining our sample with previous measurements, we apply Bayesian hierarchical analysis to a total of 26 binary systems with S-type transiting planets ($r_p<5R_{\oplus}$). Specifically, we fit the $ΔI_{\mathrm{los}}$ distribution with both single (Rayleigh) and mixture models (two-component Rayleigh and Rayleigh-isotropic mixture). We find the mixture models are strongly favored ($\log Z\gtrsim13.9$, or $\approx$5$σ$), indicating the observed planet-binary $ΔI_{\mathrm{los}}$ values likely originate from two underlying populations: one nearly aligned ($σ_1 = 2^{\circ}.4^{+0.7}_{-0.9}$) and one with more scattered mutual inclinations ($σ_2 = 23^{\circ}.6^{+8.8}_{-7.1}$). Alternatively, the misaligned systems can be equally well described by an isotropic distribution of inclinations. This observed dichotomy likely reflects different dynamical histories. Notably, the misaligned population only emerges in systems with stellar periastron distances $>40$ AU while systems with close-in or eccentric stellar companions (periastron distances $<40$ AU) preserve planet-binary alignment.

preprint2025arXiv

Multi-task Federated Learning with Encoder-Decoder Structure: Enabling Collaborative Learning Across Different Tasks

Federated learning has been extensively studied and applied due to its ability to ensure data security in distributed environments while building better models. However, clients participating in federated learning still face limitations, as clients with different structures or tasks cannot participate in learning together. In view of this, constructing a federated learning framework that allows collaboration between clients with different model structures and performing different tasks, enabling them to share valuable knowledge to enhance model efficiency, holds significant practical implications for the widespread application of federated learning. To achieve this goal, we propose a multi-task federated learning with encoder-decoder structure (M-Fed). Specifically, given the widespread adoption of the encoder-decoder architecture in current models, we leverage this structure to share intra-task knowledge through traditional federated learning methods and extract general knowledge from the encoder to achieve cross-task knowledge sharing. The training process is similar to traditional federated learning, and we incorporate local decoder and global decoder information into the loss function. The local decoder iteratively updates and gradually approaches the global decoder until sufficient cross-task knowledge sharing is achieved. Our method is lightweight and modular, demonstrating innovation compared to previous research. It enables clients performing different tasks to share general knowledge while maintaining the efficiency of traditional federated learning systems. We conducted experiments on two widely used benchmark datasets to verify the feasibility of M-Fed and compared it with traditional methods. The experimental results demonstrate the effectiveness of M-Fed in multi-task federated learning.

preprint2025arXiv

PEPSI Investigation, Retrieval, and Atlas of Numerous Giant Atmospheres (PIRANGA). IV. High-Resolution Phased-Resolved Spectroscopy of The Ultra Hot Jupiter KELT-20 b

We present five datasets of high-resolution optical emission spectra of the ultra-hot Jupiter KELT-20 b with the PEPSI spectrograph. Using a Bayesian retrieval framework, we constrain its dayside pressure-temperature profile and abundances of Fe, Ni, and Ca, providing the first measurements for Ni and Ca for KELT-20 b in emission. We retrieve the pre- and post-eclipse datasets separately (corresponding to the evening and morning sides, respectively), and compare the constraints on their thermal structures and chemical abundances. We constrain lower abundances in the pre-eclipse datasets compared to the post-eclipse datasets. We interpret these results with an equilibrium chemistry model which suggests ~10-30x supersolar refractory abundances. Due to the well-known degeneracy between absolute abundances and continuum opacities, the abundance ratios are more precise probes of the planetary abundances. Therefore we measure the abundance ratios [Ni/Fe] and [Ca/Fe] across these datasets and find they agree within 1-sigma. We constrain [Ni/Fe] to be consistent with solar within 2-sigma, and [Ca/Fe] to be 0.001-0.01x solar, not accounting for ionization. We compare these abundance ratios with literature results for KELT-20 b in transmission, and find they agree within 2-sigma, suggesting that even though the abundances vary significantly as a function of phase, the abundance ratios of these species remain relatively constant. We find a ~100 K difference in temperature at the top of the thermal inversion, suggesting a hotter evening side than morning side and underscoring the importance of considering 3D effects when studying ultra-hot Jupiters.

preprint2022arXiv

A Super-Earth Orbiting Near the Inner Edge of the Habitable Zone around the M4.5-dwarf Ross 508

We report the near-infrared radial-velocity (RV) discovery of a super-Earth planet on a 10.77-day orbit around the M4.5 dwarf Ross 508 ($J_\mathrm{mag}=9.1$). Using precision RVs from the Subaru Telescope IRD (InfraRed Doppler) instrument, we derive a semi-amplitude of $3.92^{+0.60}_{-0.58}$ ${\rm m\,s}^{-1}$, corresponding to a planet with a minimum mass $m \sin i = 4.00^{+0.53}_{-0.55}\ M_{\oplus}$. We find no evidence of significant signals at the detected period in spectroscopic stellar activity indicators or MEarth photometry. The planet, Ross 508 b, has a semimajor-axis of $0.05366^{+0.00056}_{-0.00049}$ au. This gives an orbit-averaged insolation of $\approx$1.4 times the Earth's value, placing Ross 508 b near the inner edge of its star's habitable zone. We have explored the possibility that the planet has a high eccentricity and its host is accompanied by an additional unconfirmed companion on a wide orbit. Our discovery demonstrates that the near-infrared RV search can play a crucial role to find a low-mass planet around cool M dwarfs like Ross 508.

preprint2022arXiv

A Unified Spectroscopic and Photometric Model to Infer Surface Inhomogeneity: Application to Luhman 16B

Extremely large telescopes (ELTs) provide an opportunity to observe surface inhomogeneities for ultracool objects including M dwarfs, brown dwarfs (BDs), and gas giant planets via Doppler imaging and spectro-photometry techniques. These inhomogeneities can be caused by star spots, clouds, and vortices. Star spots and associated stellar flares play a significant role in habitability, either stifling life or catalyzing abiogenesis depending on the emission frequency, magnitude, and orientation. Clouds and vortices may be the source of spectral and photometric variability observed at the L/T transition of BDs and are expected in gas giant exoplanets. We develop a versatile analytical framework to model and infer surface inhomogeneities which can be applied to both spectroscopic and photometric data. This model is validated against a slew of numerical simulations. Using archival spectroscopic and photometric data, we infer star spot parameters (location, size, and contrast) and generate global surface maps for Luhman 16B (an early T dwarf and one of our solar system's nearest neighbors at a distance of approximately 2 pc). We confirm previous findings that Luhman 16B's atmosphere is inhomogeneous with time-varying features. In addition, we provide tentative evidence of longer timescale atmospheric structures such as dark equatorial and bright mid-latitude to polar spots. These findings are discussed in the context of atmospheric circulation and dynamics for ultracool dwarfs. Our analytical model will be valuable in assessing the feasibility of using ELTs to study surface inhomogeneities of gas giant exoplanets and other ultracool objects.

preprint2022arXiv

KELT-9 as an eclipsing double-lined spectroscopic binary: a unique and self-consistent solution to the system

Transiting hot Jupiters present a unique opportunity to measure absolute planetary masses due to the magnitude of their radial velocity signals and known orbital inclination. Measuring planet mass is critical to understanding atmospheric dynamics and escape under extreme stellar irradiation. Here, we present the ultra-hot Jupiter system, KELT-9, as a double-lined spectroscopic binary. This allows us to directly and empirically constrain the mass of the star and its planetary companion, without reference to any theoretical stellar evolutionary models or empirical stellar scaling relations. Using data from the PEPSI, HARPS-N, and TRES spectrographs across multiple epochs, we apply least-squares deconvolution to measure out-of-transit stellar radial velocities. With the PEPSI and HARPS-N datasets, we measure in-transit planet radial velocities using transmission spectroscopy. By fitting the circular orbital solution that captures these Keplerian motions, we recover a planetary dynamical mass of 2.17 $\pm$ 0.56 $\mathrm{M_J}$ and stellar dynamical mass of 2.11 $\pm$ 0.78 $\mathrm{M_\odot}$, both of which agree with the discovery paper. Furthermore, we argue that this system, as well as systems like it, are highly overconstrained, providing multiple independent avenues for empirically cross-validating model-independent solutions to the system parameters. We also discuss the implications of this revised mass for studies of atmospheric escape.

preprint2022arXiv

Measuring Elemental Abundances of JWST Target Stars for Exoplanet Characterization I. FGK Stars

With the launch of the JWST, we will obtain more precise data for exoplanets than ever before. However, this data can only inform and revolutionize our understanding of exoplanets when placed in the larger context of planet-star formation. Therefore, gaining a deeper understanding of their host stars is equally important and synergistic with the upcoming JWST data. We present detailed chemical abundance profiles of 17 FGK stars that will be observed in exoplanet-focused Cycle 1 JWST observer programs. The elements analyzed (C, N, O, Na, Mg, Si, S, K, and Fe) were specifically chosen as being informative to the composition and formation of planets. Using archival high-resolution spectra from a variety of sources, we perform an LTE equivalent width analysis to derive these abundances. We look to literature sources to correct the abundances for non-LTE effects, especially for O, S, and K, where the the corrections are large (often $> 0.2~\textrm{dex}$). With these abundances and the ratios thereof, we will begin to paint clearer pictures of the planetary systems analyzed by this work. With our analysis, we can gain insight into the composition and extent of migration of Hot Jupiters, as well as the possibility of carbon-rich terrestrial worlds.

preprint2022arXiv

Nodal Precession and Tidal Evolution of Two Hot-Jupiters: WASP-33 b and KELT-9 b

Hot Jupiters orbiting rapidly rotating stars on inclined orbits undergo tidally induced nodal precession measurable over several years of observations. The Hot Jupiters WASP-33 b and KELT-9 b are particularly interesting targets as they are among the hottest planets found to date, orbiting relatively massive stars. Here, we analyze archival and new data that span 11 and 5 years for WASP-33 b and KELT-9 b, respectively, in order to to model and improve upon their tidal precession parameters. Our work confirms the nodal precession for WASP-33 b and presents the first clear detection of the precession of KELT-9 b. We determine that WASP-33 and KELT-9 have gravitational quadrupole moments $(6.3^{+1.2}_{-0.8})\times10^{-5}$ and $(3.26^{+0.93}_{-0.80})\times10^{-4}$, respectively. We estimate the planets' precession periods to be $1460^{+170}_{-130}$ years and $890^{+200}_{-140}$ years, respectively, and that they will cease to transit their host stars around the years $2090^{+17}_{-10}$~CE and $2074^{+12}_{-10}$~CE, respectively. Additionally, we investigate both planets' tidal and orbital evolution, suggesting that a high-eccentricity tidal migration scenario is possible to produce both system architectures and that they will most likely not be engulfed by their hosts before the end of their main sequence lifetimes.

preprint2022arXiv

Reach-avoid Verification Based on Convex Optimization

In this paper we propose novel optimization-based methods for verifying reach-avoid (or, eventuality) properties of continuous-time systems modelled by ordinary differential equations. Given a system, an initial set, a safe set and a target set of states, we say that the reach-avoid property holds if for all initial conditions in the initial set, any trajectory of the system starting at them will eventually, i.e.\ in unbounded yet finite time, enter the target set while remaining inside the safe set until that first target hit. Based on a discount value function, two sets of quantified constraints are derived for verifying the reach-avoid property via the computation of exponential/asymptotic guidance-barrier functions (they form a barrier escorting the system to the target set safely at an exponential or asymptotic rate). It is interesting to find that one set of constraints whose solution is termed exponential guidance-barrier functions is just a simplified version of the existing one derived from the moment based method, while the other one whose solution is termed asymptotic guidance-barrier functions is completely new. Furthermore, built upon this new set of constraints, we derive a set of more expressive constraints, which includes the aforementioned two sets of constraints as special instances, providing more chances for verifying the reach-avoid properties successfully. When the involved datum are polynomials, i.e., the initial set, safe set and target set are semi-algebraic, and the system has polynomial dynamics, the problem of solving these sets of constraints can be framed as a semi-definite optimization problem using sum-of-squares decomposition techniques and thus can be efficiently solved in polynomial time via interior point methods. Finally, several examples demonstrate the theoretical developments and performance of proposed methods.

preprint2022arXiv

Retrieving the C and O Abundances of HR 7672~AB: a Solar-Type Primary Star with a Benchmark Brown Dwarf

A benchmark brown dwarf (BD) is a BD whose properties (e.g., mass and chemical composition) are precisely and independently measured. Benchmark BDs are valuable in testing theoretical evolutionary tracks, spectral synthesis, and atmospheric retrievals for sub-stellar objects. Here, we report results of atmospheric retrieval on a synthetic spectrum and a benchmark BD -- HR 7672~B -- with \petit. First, we test the retrieval framework on a synthetic PHOENIX BT-Settl spectrum with a solar composition. We show that the retrieved C and O abundances are consistent with solar values, but the retrieved C/O is overestimated by 0.13-0.18, which is $\sim$4 times higher than the formal error bar. Second, we perform retrieval on HR 7672~B using high spectral resolution data (R=35,000) from the Keck Planet Imager and Characterizer (KPIC) and near infrared photometry. We retrieve [C/H], [O/H], and C/O to be $-0.24\pm0.05$, $-0.19\pm0.04$, and $0.52\pm0.02$. These values are consistent with those of HR 7672~A within 1.5-$σ$. As such, HR 7672~B is among only a few benchmark BDs (along with Gl 570~D and HD 3651~B) that have been demonstrated to have consistent elemental abundances with their primary stars. Our work provides a practical procedure of testing and performing atmospheric retrieval, and sheds light on potential systematics of future retrievals using high- and low-resolution data.

preprint2022arXiv

The PEPSI Exoplanet Transit Survey (PETS) I: Investigating the presence of a silicate atmosphere on the super-Earth 55 Cnc e

The study of exoplanets and especially their atmospheres can reveal key insights on their evolution by identifying specific atmospheric species. For such atmospheric investigations, high-resolution transmission spectroscopy has shown great success, especially for Jupiter-type planets. Towards the atmospheric characterization of smaller planets, the super-Earth exoplanet 55 Cnc e is one of the most promising terrestrial exoplanets studied to date. Here, we present a high-resolution spectroscopic transit observation of this planet, acquired with the PEPSI instrument at the Large Binocular Telescope. Assuming the presence of Earth-like crust species on the surface of 55 Cnc e, from which a possible silicate-vapor atmosphere could have originated, we search in its transmission spectrum for absorption of various atomic and ionized species such as Fe , Fe+, Ca , Ca+, Mg and K , among others. Not finding absorption for any of the investigated species, we are able to set absorption limits with a median value of 1.9 x RP. In conclusion, we do not find evidence of a widely extended silicate envelope on this super-Earth reaching several planetary radii.

preprint2022arXiv

Variable and super-sonic winds in the atmosphere of an ultra-hot giant planet

Hot Jupiters receive intense irradiation from their stellar hosts. The resulting extreme environments in their atmospheres allow us to study the conditions that drive planetary atmospheric dynamics, e.g., global-scale winds. General circulation models predict day-to-nightside winds and equatorial jets with speeds on the order of a few km $\mathrm{s^{-1}}$. To test these models, we apply high-resolution transmission spectroscopy using the PEPSI spectrograph on the Large Binocular Telescope to study the atmosphere of KELT-9 b, an ultra-hot Jupiter and currently the hottest known planet. We measure $\sim$10 km $\mathrm{s^{-1}}$ day-to-nightside winds traced by Fe II features in the planet's atmosphere. This is at odds with previous literature (including data taken with PEPSI), which report no significant day-to-nightside winds on KELT-9 b. We identify the cause of this discrepancy as due to an inaccurate ephemeris for KELT-9 b in previous literature. We update the ephemeris, which shifts the mid-transit time by up to 10 minutes for previous datasets, resulting in consistent detections of blueshifts in all the datasets analyzed here. Furthermore, a comparison with archival HARPS-N datasets suggests temporal wind variability $\sim$5-8 km $\mathrm{s^{-1}}$ over timescales between weeks to years. Temporal variability of atmospheric dynamics on hot Jupiters is a phenomenon anticipated by certain general circulation models that has not been observed over these timescales until now. However, such large variability as we measure on KELT-9 b challenges general circulation models, which predict much lower amplitudes of wind variability over timescales between days to weeks.

preprint2021arXiv

DepOwl: Detecting Dependency Bugs to Prevent Compatibility Failures

Applications depend on libraries to avoid reinventing the wheel. Libraries may have incompatible changes during evolving. As a result, applications will suffer from compatibility failures. There has been much research on addressing detecting incompatible changes in libraries, or helping applications co-evolve with the libraries. The existing solution helps the latest application version work well against the latest library version as an afterthought. However, end users have already been suffering from the failures and have to wait for new versions. In this paper, we propose DepOwl, a practical tool helping users prevent compatibility failures. The key idea is to avoid using incompatible versions from the very beginning. We evaluated DepOwl on 38 known compatibility failures from StackOverflow, and DepOwl can prevent 32 of them. We also evaluated DepOwl using the software repository shipped with Ubuntu-19.10. DepOwl detected 77 unknown dependency bugs, which may lead to compatibility failures.

preprint2021arXiv

Detecting Biosignatures in the Atmospheres of Gas Dwarf Planets with the James Webb Space Telescope

Exoplanets with radii between those of Earth and Neptune have stronger surface gravity than Earth, and can retain a sizable hydrogen-dominated atmosphere. In contrast to gas giant planets, we call these planets gas dwarf planets. The James Webb Space Telescope (JWST) will offer unprecedented insight into these planets. Here, we investigate the detectability of ammonia (NH$_{3}$, a potential biosignature) in the atmospheres of seven temperate gas dwarf planets using various JWST instruments. We use petitRadTRANS and PandExo to model planet atmospheres and simulate JWST observations under different scenarios by varying cloud conditions, mean molecular weights (MMWs), and NH$_{3}$ mixing ratios. A metric is defined to quantify detection significance and provide a ranked list for JWST observations in search of biosignatures in gas dwarf planets. It is very challenging to search for the 10.3--10.8 $μ$m NH$_{3}$ feature using eclipse spectroscopy with MIRI in the presence of photon and a systemic noise floor of 12.6 ppm for 10 eclipses. NIRISS, NIRSpec, and MIRI are feasible for transmission spectroscopy to detect NH$_{3}$ features from 1.5 $μ$m to 6.1 $μ$m under optimal conditions such as a clear atmosphere and low MMWs for a number of gas dwarf planets. We provide examples of retrieval analyses to further support the detection metric that we use. Our study shows that searching for potential biosignatures such as NH$_{3}$ is feasible with a reasonable investment of JWST time for gas dwarf planets given optimal atmospheric conditions.

preprint2021arXiv

Elemental abundances of nearby M dwarfs based on high-resolution near-infrared spectra obtained by the Subaru/IRD survey: Proof of concept

Detailed chemical analyses of M dwarfs are scarce but necessary to constrain the formation environment and internal structure of planets being found around them. We present elemental abundances of 13 M dwarfs (2900 < Teff < 3500 K) observed in the Subaru/IRD planet search project. They are mid-to-late M dwarfs whose abundance of individual elements has not been well studied. We use the high-resolution (~70,000) near-infrared (970-1750 nm) spectra to measure the abundances of Na, Mg, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Sr by the line-by-line analysis based on model atmospheres, with typical errors ranging from 0.2 dex for [Fe/H] to 0.3-0.4 dex for other [X/H]. We measure radial velocities from the spectra and combine them with Gaia astrometry to calculate the Galactocentric space velocities UVW. The resulting [Fe/H] values agree with previous estimates based on medium-resolution K-band spectroscopy, showing a wide distribution of metallicity (-0.6 < [Fe/H] < +0.4). The abundance ratios of individual elements [X/Fe] are generally aligned with the solar values in all targets. While the [X/Fe] distributions are comparable to those of nearby FGK stars, most of which belong to the thin disk population, the most metal-poor object, GJ 699, could be a thick disk star. The UVW velocities also support this. The results raise the prospect that near-infrared spectra of M dwarfs obtained in the planet search projects can be used to grasp the trend of elemental abundances and Galactic stellar population of nearby M dwarfs.

preprint2021arXiv

Time-resolved rotational velocities in the upper atmosphere of WASP-33 b

While steady empirical progress has been made in understanding the structure and composition of hot planet atmospheres, direct measurements of velocity signatures, including winds, rotation, and jets, have lagged behind. Quantifying atmospheric dynamics of hot planets is critical to a complete understanding of their atmospheres and such measurements may even illuminate other planetary properties, such as magnetic field strengths. In this manuscript we present the first detection of the Balmer lines H$α$ and H$β$ in the atmosphere of the ultra-hot Jupiter WASP-33 b. Using atmospheric models which include the effects of atmospheric dynamics, we show that the shape of the average Balmer line transmission spectrum is consistent with rotational velocities in the planet's thermosphere of $v_\text{rot} = 10.1^{+0.8}_{-1.0}$ km s$^{-1}$. We also measure a low-significance day-to-night side velocity shift of $-4.6^{+3.4}_{-3.4}$ km s$^{-1}$ in the transmission spectrum which is naturally explained by a global wind across the planet's terminator. In a separate analysis the time-resolved velocity centroids of individual transmission spectra show unambiguous evidence of rotation, with a best-fit velocity of $10.0^{+2.4}_{-2.0}$ km s$^{-1}$, consistent with the value of $v_\text{rot}$ derived from the shape of the average Balmer line transmission spectrum. Our observations and analysis confirm the power of high signal-to-noise, time-resolved transmission spectra to measure the velocity structures in exoplanet atmospheres. The large rotational and wind velocities we measure highlight the need for more detailed 3D global climate simulations of the rarefied upper-atmospheres of ultra-hot gas giants.

preprint2020arXiv

Combining Symbolic Execution and Model Checking to Verify MPI Programs

Message passing is the standard paradigm of programming in high-performance computing. However, verifying Message Passing Interface (MPI) programs is challenging, due to the complex program features (such as non-determinism and non-blocking operations). In this work, we present MPI symbolic verifier (MPI-SV), the first symbolic execution based tool for automatically verifying MPI programs with non-blocking operations. MPI-SV combines symbolic execution and model checking in a synergistic way to tackle the challenges in MPI program verification. The synergy improves the scalability and enlarges the scope of verifiable properties. We have implemented MPI-SV (footnote: https://mpi-sv.github.io) and evaluated it with 111 real-world MPI verification tasks. The pure symbolic execution-based technique successfully verifies 61 out of the 111 tasks (55\%) within one hour, while in comparison, MPI-SV verifies 100 tasks (90\%). On average, compared with pure symbolic execution, MPI-SV achieves 19x speedups on verifying the satisfaction of the critical property and 5x speedups on finding violations.

preprint2020arXiv

Delay-Compensated Control of Sandwiched ODE-PDE-ODE Hyperbolic Systems for Oil Drilling and Disaster Relief

Motivated by engineering applications of subsea installation by deepwater construction vessels in oil drilling, and of aid delivery by unmanned aerial vehicles in disaster relief, we develop output-feedback boundary control of heterodirectional coupled hyperbolic PDEs sandwiched between two ODEs, where the measurement is the output state of one ODE and suffers a time delay. After rewriting the time-delay dynamics as a transport PDE of which the left boundary connects with the sandwiched system, a state observer is built to estimate the states of the overall system of ODE-heterodirectional coupled hyperbolic PDEs-ODE-transport PDE using the right boundary state of the last transport PDE. An observer-based output-feedback controller acting at the first ODE is designed to stabilize the overall system using backstepping transformations and frequency-domain designs. The exponential stability results of the closed-loop system, boundedness and exponential convergence of the control input are proved. The obtained theoretical result is applied to control of a deepwater oil drilling construction vessel as a simulation case, where the simulation results show the proposed control design reduces cable oscillations and places the oil drilling equipment to be installed in the target area on the sea floor. Performance deterioration under extreme and unmodeled disturbances is also illustrated.

preprint2020arXiv

EarthFinder Probe Mission Concept Study: Characterizing nearby stellar exoplanet systems with Earth-mass analogs for future direct imaging

EarthFinder is a NASA Astrophysics Probe mission concept selected for study as input to the 2020 Astrophysics National Academies Decadal Survey. The EarthFinder concept is based on a dramatic shift in our understanding of how PRV measurements should be made. We propose a new paradigm which brings the high precision, high cadence domain of transit photometry as demonstrated by Kepler and TESS to the challenges of PRV measurements at the cm/s level. This new paradigm takes advantage of: 1) broad wavelength coverage from the UV to NIR which is only possible from space to minimize the effects of stellar activity; 2) extremely compact, highly stable, highly efficient spectrometers (R>150,000) which require the diffraction-limited imaging possible only from space over a broad wavelength range; 3) the revolution in laser-based wavelength standards to ensure cm/s precision over many years; 4) a high cadence observing program which minimizes sampling-induced period aliases; 5) exploiting the absolute flux stability from space for continuum normalization for unprecedented line-by-line analysis not possible from the ground; and 6) focusing on the bright stars which will be the targets of future imaging missions so that EarthFinder can use a ~1.5 m telescope.

preprint2020arXiv

Keck/NIRC2 $L$'-Band Imaging of Jovian-Mass Accreting Protoplanets around PDS 70

We present $L$'-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wavefront sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semi-major axis of $20^{+3}_{-4}$~au and PDS 70 c to have a semi-major axis of $34^{+12}_{-6}$~au (95\% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2-3~$R_{Jup}$. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 $M_{\textrm{Jup}}$ and a mean mass accretion rate between $3 \times 10^{-7}$ and $8 \times 10^{-7}~M_{\textrm{Jup}}/\textrm{yr}$. For PDS 70 c, we computed a mass between 1 and 3 $M_{\textrm{Jup}}$ and mean mass accretion rate between $1 \times 10^{-7}$ and $5 \times~10^{-7} M_{\textrm{Jup}}/\textrm{yr}$. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets' SEDs.

preprint2020arXiv

Not Just Cloud Privacy: Protecting Client Privacy in Teacher-Student Learning

Ensuring the privacy of sensitive data used to train modern machine learning models is of paramount importance in many areas of practice. One recent popular approach to study these concerns is using the differential privacy via a "teacher-student" model, wherein the teacher provides the student with useful, but noisy, information, hopefully allowing the student model to perform well on a given task. However, these studies only solve the privacy concerns of the teacher by assuming the student owns a public but unlabelled dataset. In real life, the student also has privacy concerns on its unlabelled data, so as to inquire about privacy protection on any data sent to the teacher. In this work, we re-design the privacy-preserving "teacher-student" model consisting of adopting both private arbitrary masking and local differential privacy, which protects the sensitive information of each student sample. However, the traditional training of teacher model is not robust on any perturbed data. We use the adversarial learning techniques to improve the robustness of the perturbed sample that supports returning good feedback without having all private information of each student sample. The experimental results demonstrate the effectiveness of our new privacy-preserving "teacher-student" model.

preprint2020arXiv

On the Chemical Abundance of HR 8799 and the Planet c

Comparing chemical abundances of a planet and the host star reveals the origin and formation path. Stellar abundance is measured with high-resolution spectroscopy. Planet abundance, on the other hand, is usually inferred from low-resolution data. For directly imaged exoplanets, the data are available from a slew of high-contrast imaging/spectroscopy instruments. Here, we study the chemical abundance of HR 8799 and its planet c. We measure stellar abundance using LBT/PEPSI (R=120,000) and archival HARPS data: stellar [C/H], [O/H], and C/O are 0.11$\pm$0.12, 0.12$\pm$0.14, and 0.54$^{+0.12}_{-0.09}$, all consistent with solar values. We conduct atmospheric retrieval using newly obtained Subaru/CHARIS data together with archival Gemini/GPI and Keck/OSIRIS data. We model the planet spectrum with petitRADTRANS and conduct retrieval using PyMultiNest. Retrieved planetary abundance can vary by $\sim$0.5 dex, from sub-stellar to stellar C and O abundances. The variation depends on whether strong priors are chosen to ensure a reasonable planet mass. Moreover, comparison with previous works also reveals inconsistency in abundance measurements. We discuss potential issues that can cause the inconsistency, e.g., systematics in individual data sets and different assumptions in the physics and chemistry in retrieval. We conclude that no robust retrieval can be obtained unless the issues are fully resolved.

preprint2020arXiv

SIR Model-based Prediction of Infected Population of Coronavirus in Hubei Province

After the sudden outbreak of Coronavirus in Wuhan, continuous and rich data of the epidemic has been made public as the vital fact for decision support in control measures and aggressive implementation of containment strategies and plans. With the further growth and spreading of the virus, future resource allocation and planning under updated strategies and measures rely on careful study of the epidemic data and characteristics for accurate prediction and estimation. By using the SIR model and reported data, key parameters are obtained from least sum of squared errors for an accurate prediction of epidemic trend in the last four weeks.

preprint2020arXiv

Vibration Suppression for Coupled Wave PDEs in Deep-Sea Construction

A deep-sea construction vessel is used to install underwater parts of an off-shore oil drilling platform at the designated locations on the seafloor. By using extended Hamilton's principle, a nonlinear PDE system governing the lateral-longitudinal coupled vibration dynamics of the deep-sea construction vessel consisting of a time-varying-length cable with an attached item is derived, and it is linearized at the steady state generating a linear PDE model, which is extended to a more general system including two coupled wave PDEs connected with two interacting ODEs at the uncontrolled boundaries. Through a preliminary transformation, an equivalent reformulated plant is generated as a 4*4 coupled heterodirectional hyperbolic PDE-ODE system characterized by spatially-varying coefficients on a time-varying domain. To stabilize such a system, an observer-based output-feedback control design is proposed, where the measurements are only placed at the actuated boundary of the PDE, namely, at the platform at the sea surface. The exponential stability of the closed-loop system, boundedness and exponential convergence of the control inputs, are proved via Lyapunov analysis. The obtained theoretical result is tested on a nonlinear model with ocean disturbances, even though the design is developed in the absence of such real-world effects.

preprint2016arXiv

Characterization of the Companion to $μ$ Her

$μ$ Her is a nearby quadruple system with a G-subgiant primary and several low mass companions arranged in a 2+2 architecture. While the BC components have been well characterized, the Ab component has been detected astrometrically and with direct imaging but there has been some confusion over its nature, in particular whether the companion is stellar or substellar. Using near-infrared spectroscopy we are able to estimate the spectral type of the companion as a M4$\pm$1V star. In addition, we have measured the astrometry of the system for over a decade. We combined the astrometry with archival radial velocity measurements to compute an orbit of the system. From the combined orbit, we are able to compute the mass sum of the system. Using the estimated mass of the primary, we estimate the mass of the secondary as 0.32 M_sun, which agrees with the estimated spectral type. Our computed orbit is preliminary due to the incomplete orbital phase coverage, but it should be sufficient to predict ephemerides over the next decade.

preprint2016arXiv

Dependence of Small Planet Frequency on Stellar Metallicity Hidden by Their Prevalence

The dependence of gas giant planet occurrence rate on stellar metallicity has been firmly established. We extend this so-called planet-metallicity correlation to broader ranges of metallicities and planet masses/radii. In particular, we assume that the planet-metallicity correlation is a power law below some critical saturation threshold, and that the probability of hosting at least one planet is unity for stars with metallicity above the threshold. We then are able to explain the discrepancy between the tentative detection and null detection in previous studies regarding the planet-metallicity correlation for small planets. In particular, we find that the null detection of this correlation can be attributed to the combination of high planet occurrence rate and low detection efficiency. Therefore, a planet-metallicity correlation for small planets cannot be ruled out. We propose that stars with metallicities lower than the Solar value are better targets for testing the planet-metallicity correlation for small planets.

preprint2016arXiv

Friends of Hot Jupiters. IV. Stellar companions beyond 50 AU might facilitate giant planet formation, but most are unlikely to cause Kozai-Lidov migration

Stellar companions can influence the formation and evolution of planetary systems, but there are currently few observational constraints on the properties of planet-hosting binary star systems. We search for stellar companions around 77 transiting hot Jupiter systems to explore the statistical properties of this population of companions as compared to field stars of similar spectral type. After correcting for survey incompleteness, we find that $47\%\pm7\%$ of hot Jupiter systems have stellar companions with semi-major axes between 50-2000 AU. This is 2.9 times larger than the field star companion fraction in this separation range, with a significance of $4.4σ$. In the 1-50AU range, only $3.9^{+4.5}_{-2.0}\%$ of hot Jupiters host stellar companions compared to the field star value of $16.4\%\pm0.7\%$, which is a $2.7σ$ difference. We find that the distribution of mass ratios for stellar companions to hot Jupiter systems peaks at small values and therefore differs from that of field star binaries which tend to be uniformly distributed across all mass ratios. We conclude that either wide separation stellar binaries are more favorable sites for gas giant planet formation at all separations, or that the presence of stellar companions preferentially causes the inward migration of gas giant planets that formed farther out in the disk via dynamical processes such as Kozai-Lidov oscillations. We determine that less than 20% of hot Jupiters have stellar companions capable of inducing Kozai-Lidov oscillations assuming initial semi-major axes between 1-5 AU, implying that the enhanced companion occurrence is likely correlated with environments where gas giants can form efficiently.

preprint2016arXiv

Planet Hunters X: Searching for Nearby Neighbors of 75 Planet and Eclipsing Binary Candidates from the K2 Kepler Extended Mission

We present high-resolution observations of a sample of 75 K2 targets from Campaigns 1-3 using speckle interferometry on the Southern Astrophysical Research (SOAR) telescope and adaptive optics (AO) imaging at the Keck II telescope. The median SOAR $I$-band and Keck $K_s$-band detection limits at 1" were $Δm_{I}=4.4$~mag and $Δm_{K_s}=6.1$~mag, respectively. This sample includes 37 stars likely to host planets, 32 targets likely to be eclipsing binaries (EBs), and 6 other targets previously labeled as likely planetary false positives. We find nine likely physically bound companion stars within 3" of three candidate transiting exoplanet host stars and six likely EBs. Six of the nine detected companions are new discoveries; one of the six, EPIC 206061524, is associated with a planet candidate. Among the EB candidates, companions were only found near the shortest period ones ($P<3$ days), which is in line with previous results showing high multiplicity near short-period binary stars. This high-resolution data, including both the detected companions and the limits on potential unseen companions, will be useful in future planet vetting and stellar multiplicity rate studies for planets and binaries.

preprint2016arXiv

Synthesis, Optical, and Magnetic Properties of Ba$_2$Ni$_3$F$_{10}$ Nanowires

A low temperature hydrothermal route has been developed, and pure phase Ba$_2$Ni$_3$F$_{10}$ nanowires have been successfully prepared under the optimized conditions. Under the 325 nm excitation, the Ba$_2$Ni$_3$F$_{10}$ nanowires exhibit three emission bands with peak positions locating at 360 nm, 530 nm, and 700 nm, respectively. Combined with the first-principles calculations, the photoluminescence property can be explained by the electron transitions between the t2g and eg orbitals. Clear hysteresis loops observed below the temperature of 60 K demonstrates the weak ferromagnetism in Ba$_2$Ni$_3$F$_{10}$ nanowires, which has been attributed to the surface strain of nanowires. Exchange bias with blocking temperature of 55 K has been observed, which originates from the magnetization pinning under the cooling field due to antiferromagnetic core/weak ferromagnetic shell structure of Ba2Ni3F10 nanowires.

preprint2016arXiv

Very Low-Mass Stellar and Substellar Companions to Solar-like Stars From MARVELS VI: A Giant Planet and a Brown Dwarf Candidate in a Close Binary System HD 87646

We report the detections of a giant planet (MARVELS-7b) and a brown dwarf candidate (MARVELS-7c) around the primary star in the close binary system, HD 87646. It is the first close binary system with more than one substellar circum-primary companion discovered to the best of our knowledge. The detection of this giant planet was accomplished using the first multi-object Doppler instrument (KeckET) at the Sloan Digital Sky Survey (SDSS) telescope. Subsequent radial velocity observations using ET at Kitt Peak National Observatory, HRS at HET, the "Classic" spectrograph at the Automatic Spectroscopic Telescope at Fairborn Observatory, and MARVELS from SDSS-III confirmed this giant planet discovery and revealed the existence of a long-period brown dwarf in this binary. HD 87646 is a close binary with a separation of $\sim22$ AU between the two stars, estimated using the Hipparcos catalogue and our newly acquired AO image from PALAO on the 200-inch Hale Telescope at Palomar. The primary star in the binary, HD 87646A, has Teff = 5770$\pm$80K, log(g)=4.1$\pm$0.1 and [Fe/H] = $-0.17\pm0.08$. The derived minimum masses of the two substellar companions of HD 87646A are 12.4$\pm$0.7M$_{\rm Jup}$ and 57.0$\pm3.7$M$_{\rm Jup}$. The periods are 13.481$\pm$0.001 days and 674$\pm$4 days and the measured eccentricities are 0.05$\pm$0.02 and 0.50$\pm$0.02 respectively. Our dynamical simulations show the system is stable if the binary orbit has a large semi-major axis and a low eccentricity, which can be verified with future astrometry observations.

preprint2015arXiv

A Simple Probabilistic Extension of Modal Mu-calculus

Probabilistic systems are an important theme in AI domain. As the specification language, the logic PCTL is now the default logic for reasoning about probabilistic properties. In this paper, we present a natural and succinct probabilistic extension of mu-calculus, another prominent logic in the concurrency theory. We study the relationship with PCTL. Interestingly, the expressiveness is highly orthogonal with PCTL. The proposed logic captures some useful properties which cannot be expressed in PCTL. We investigate the model checking and satisfiability problem, and show that the model checking problem is in UP and co-UP, and the satisfiability checking can be decided via reducing into solving parity games. This is in contrast to PCTL as well, whose satisfiability checking is still an open problem.

preprint2015arXiv

An Analysis of Thickness-shear Vibrations of an Annular Plate with the Mindlin Plate Equations

The Mindlin plate equations with the consideration of thickness-shear deformation as an independent variable have been used for the analysis of vibrations of quartz crystal resonators of both rectangular and circular types. The Mindlin or Lee plate theories that treat thickness-shear deformation as an independent higher-order vibration mode in a coupled system of two-dimensional variables are the choice of theory for analysis. For circular plates, we derived the Mindlin plate equations in a systematic manner as demonstrated by Mindlin and others and obtained the truncated two-dimensional equations of closely coupled modes in polar coordinates. We simplified the equations for vibration modes in the vicinity of fundamental thickness-shear frequency and validated the equations and method. To explore newer structures of quartz crystal resonators, we utilized the Mindlin plate equations for the analysis of annular plates with fixed inner and free outer edges for frequency spectra. The detailed analysis of vibrations of circular plates for the normalized frequency versus dimensional parameters provide references for optimal selection of parameters based on the principle of strong thickness-shear mode and minimal presence of other modes to enhance energy trapping through maintaining the strong and pure thickness-shear vibrations insensitive to some complication factors such as thermal and initial stresses.

preprint2015arXiv

Dynamical Masses of Young M Dwarfs: Masses and Orbital Parameters of GJ 3305 AB, the Wide Binary Companion to the Imaged Exoplanet Host 51 Eri

We combine new high resolution imaging and spectroscopy from Keck/NIRC2, Discovery Channel Telescope/DSSI, and Keck/HIRES with published astrometry and radial velocities to measure individual masses and orbital elements of the GJ 3305 AB system, a young (~20 Myr) M+M binary (unresolved spectral type M0) member of the beta Pictoris moving group comoving with the imaged exoplanet host 51 Eri. We measure a total system mass of 1.11 \pm 0.04 M_sun, a period of 29.03 \pm 0.50$ yr, a semimajor axis of 9.78 \pm 0.14 AU, and an eccentricity of 0.19 \pm 0.02. The primary component has a dynamical mass of 0.67 \pm 0.05 M_sun and the secondary has a mass of 0.44 \pm 0.05 M_sun. The recently updated BHAC15 models are consistent with the masses of both stars to within 1.5 sigma. Given the observed masses the models predict an age of the GJ 3305 AB system of 37 \pm 9 Myr. Based on the the observed system architecture and our dynamical mass measurement, it is unlikely that the orbit of 51 Eri b has been significantly altered by the Kozai-Lidov mechanism.

preprint2015arXiv

Influence of Stellar Multiplicity On Planet Formation. III. Adaptive Optics Imaging of Kepler Stars With Gas Giant Planets

As hundreds of gas giant planets have been discovered, we study how these planets form and evolve in different stellar environments, specifically in multiple stellar systems. In such systems, stellar companions may have a profound influence on gas giant planet formation and evolution via several dynamical effects such as truncation and perturbation. We select 84 Kepler Objects of Interest (KOIs) with gas giant planet candidates. We obtain high-angular resolution images using telescopes with adaptive optics (AO) systems. Together with the AO data, we use archival radial velocity data and dynamical analysis to constrain the presence of stellar companions. We detect 59 stellar companions around 40 KOIs for which we develop methods of testing their physical association. These methods are based on color information and galactic stellar population statistics. We find evidence of suppressive planet formation within 20 AU by comparing stellar multiplicity. The stellar multiplicity rate for planet host stars is 0$^{+5}_{-0}$\% within 20 AU. In comparison, the stellar multiplicity rate is 18\%$\pm$2\% for the control sample, i.e., field stars in the solar neighborhood. The stellar multiplicity rate for planet host stars is 34\%$\pm$8\% for separations between 20 and 200 AU, which is higher than the control sample at 12\%$\pm$2\%. Beyond 200 AU, stellar multiplicity rates are comparable between planet host stars and the control sample. We discuss the implications of the results to gas giant planet formation and evolution.

preprint2015arXiv

Influence of Stellar Multiplicity On Planet Formation. IV. Adaptive Optics Imaging of Kepler Stars With Multiple Transiting Planet Candidates

The Kepler mission provides a wealth of multiple transiting planet systems (MTPS). The formation and evolution of multi-planet systems are likely to be influenced by companion stars given the abundance of multi stellar systems. We study the influence of stellar companions by measuring the stellar multiplicity rate of MTPS. We select 138 bright (KP < 13.5) Kepler MTPS and search for stellar companions with AO imaging data and archival radial velocity (RV) data. We obtain new AO images for 73 MTPS. Other MTPS in the sample have archival AO imaging data from the Kepler Community Follow-up Observation Program (CFOP). From these imaging data, we detect 42 stellar companions around 35 host stars. For stellar separation 1 AU < a < 100 AU, the stellar multiplicity rate is 5.2 $\pm$ 5.0% for MTPS, which is 2.8σ lower than 21.1 $\pm$ 2.8% for the control sample, i.e., the field stars in the solar neighborhood. We identify two origins for the deficit of stellar companions within 100 AU to MTPS: (1) a suppressive planet formation, and (2) the disruption of orbital coplanarity due to stellar companions. To distinguish between the two origins, we compare the stellar multiplicity rates of MTPS and single transiting planet systems (STPS). However, current data are not sufficient for this purpose. For 100 AU < a < 2000 AU, the stellar multiplicity rates are comparable for MTPS (8.0 $\pm$ 4.0%), STPS (6.4 $\pm$ 5.8%), and the control sample (12.5 $\pm$ 2.8%).

preprint2015arXiv

KELT-4Ab: An inflated Hot Jupiter transiting the bright (V~10) component of a hierarchical triple

We report the discovery of KELT-4Ab, an inflated, transiting Hot Jupiter orbiting the brightest component of a hierarchical triple stellar system. The host star is an F star with $T_{\rm eff}=6206\pm75$ K, $\log g=4.108\pm0.014$, $\left[{\rm Fe}/{\rm H}\right]=-0.116_{-0.069}^{+0.065}$, ${\rm M_*}=1.201_{-0.061}^{+0.067} \ {\rm M}_{\odot}$, and ${\rm R_*}=1.610_{-0.068}^{+0.078} \ {\rm R}_{\odot}$. The best-fit linear ephemeris is $\rm {BJD_{TDB}} = 2456193.29157 \pm 0.00021 + E\left(2.9895936 \pm 0.0000048\right)$. With a magnitude of $V\sim10$, a planetary radius of $1.699_{-0.045}^{+0.046} \ {\rm R_J}$, and a mass of $0.902_{-0.059}^{+0.060} \ {\rm M_J}$, it is the brightest host among the population of inflated Hot Jupiters ($R_P > 1.5R_J$), making it a valuable discovery for probing the nature of inflated planets. In addition, its existence within a hierarchical triple and its proximity to Earth ($210$ pc) provides a unique opportunity for dynamical studies with continued monitoring with high resolution imaging and precision radial velocities. In particular, the motion of the binary stars around each other and of both stars around the primary star relative to the measured epoch in this work should be detectable when it rises in October 2015.

preprint2015arXiv

KELT-7b: A hot Jupiter transiting a bright V=8.54 rapidly rotating F-star

We report the discovery of KELT-7b, a transiting hot Jupiter with a mass of $1.28 \pm 0.18$ MJ, radius of $1.53_{-0.047}^{+0.046}$ RJ, and an orbital period of $2.7347749 \pm 0.0000039$ days. The bright host star (HD33643; KELT-7) is an F-star with $V=8.54$, Teff $=6789_{-49}^{+50}$ K, [Fe/H] $=0.139_{-0.081}^{+0.075}$, and $\log{g}=4.149 \pm 0.019$. It has a mass of $1.535_{-0.054}^{+0.066}$ Msun, a radius of $1.732_{-0.045}^{+0.043}$ Rsun, and is the fifth most massive, fifth hottest, and the ninth brightest star known to host a transiting planet. It is also the brightest star around which KELT has discovered a transiting planet. Thus, KELT-7b is an ideal target for detailed characterization given its relatively low surface gravity, high equilibrium temperature, and bright host star. The rapid rotation of the star ($73 \pm 0.5$ km/s) results in a Rossiter-McLaughlin effect with an unusually large amplitude of several hundred m/s. We find that the orbit normal of the planet is likely to be well-aligned with the stellar spin axis, with a projected spin-orbit alignment of $λ=9.7 \pm 5.2$ degrees. This is currently the second most rapidly rotating star to have a reflex signal (and thus mass determination) due to a planetary companion measured.

preprint2015arXiv

Kepler Eclipsing Binary Stars. VI. Identification of Eclipsing Binaries in the K2 Campaign 0 Data-set

The original {\it Kepler} mission observed and characterized over 2400 eclipsing binaries in addition to its prolific exoplanet detections. Despite the mechanical malfunction and subsequent non-recovery of two reaction wheels used to stabilize the instrument, the {\it Kepler} satellite continues collecting data in its repurposed {\it K2} mission surveying a series of fields along the ecliptic plane. Here we present an analysis of the first full baseline {\it K2} data release: the Campaign 0 data-set. In the 7761 light curves, we have identified a total of 207 eclipsing binaries. Of these, 97 are new discoveries that were not previously identified. Our pixel-level analysis of these objects has also resulted in identification of several false positives (observed targets contaminated by neighboring eclipsing binaries), as well as the serendipitous discovery of two short period exoplanet candidates. We provide catalog cross-matched source identifications, orbital periods, morphologies and ephemerides for these eclipsing systems. We also describe the incorporation of the K2 sample into the Kepler Eclipsing Binary Catalog\footnote{\url{keplerebs.villanova.edu/k2}}, present spectroscopic follow-up observations for a limited selection of nine systems, and discuss prospects for upcoming {\it K2} campaigns.

preprint2015arXiv

Magnetoelectricity coupled exchange bias in BaMnF$_4$

Multiferroic BaMnF$_4$ powder were prepared by hydrothermal method. Hysteretic field dependent magnetization curve at 5 K confirms the weak ferromagnetism aroused from the canted antiferromagnetic spins by magnetoelectric coupling. The blocking temperature of 65 K for exchange bias coincides well with the peak at 65 K in the zero-field cooled temperature-dependent magnetization curve, which has been assigned to the onset temperature of two-dimensional antiferromagnetism. An upturn kink of exchange field and coercivity with decreasing temperature was observed from 40 K to 20 K, which is consistent with the two-dimensional to three-dimensional antiferromagnetic transition at Néel temperature (~26 K). In contrast to the conventional mechanism of magnetization pinned by interfacial exchange coupling in multiphases, the exchange bias in BaMnF$_4$ is argued to be a bulk effect in single phase, due to the magnetization pinned by the polarization through magnetoelectric coupling.

preprint2015arXiv

The APOGEE Spectroscopic Survey of Kepler Planet Hosts: Feasibility, Efficiency, and First Results

The Kepler mission has yielded a large number of planet candidates from among the Kepler Objects of Interest (KOIs), but spectroscopic follow-up of these relatively faint stars is a serious bottleneck in confirming and characterizing these systems. We present motivation and survey design for an ongoing project with the SDSS-III multiplexed APOGEE near-infrared spectrograph to monitor hundreds of KOI host stars. We report some of our first results using representative targets from our sample, which include current planet candidates that we find to be false positives, as well as candidates listed as false positives that we do not find to be spectroscopic binaries. With this survey, KOI hosts are observed over ~20 epochs at a radial velocity precision of 100-200 m/s. These observations can easily identify a majority of false positives caused by physically-associated stellar or substellar binaries, and in many cases, fully characterize their orbits. We demonstrate that APOGEE is capable of achieving RV precision at the 100-200 m/s level over long time baselines, and that APOGEE's multiplexing capability makes it substantially more efficient at identifying false positives due to binaries than other single-object spectrographs working to confirm KOIs as planets. These APOGEE RVs enable ancillary science projects, such as studies of fundamental stellar astrophysics or intrinsically rare substellar companions. The coadded APOGEE spectra can be used to derive stellar properties (T_eff, log(g)) and chemical abundances of over a dozen elements to probe correlations of planet properties with individual elemental abundances.

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

preprint2015arXiv

Thickness-shear Vibration Frequencies of an Infinite Plate with a Generalized Material Property Grading along the Thickness

For quartz crystal resonators of thickness-shear type, the vibration frequency and mode shapes, which are key features of resonators in circuit applications, reflect the basic material and structural properties of the quartz plate and its variation with time under various factors such as erosive gases and liquids that can cause surface and internal damages and degradation of crystal blanks. The accumulated effects eventually will change the surface conditions in terms of elastic constants and stiffness and more importantly, the gradient of such properties along the thickness. This is a typical functionally graded materials (FGM) structure and has been studied extensively for structural applications under multiple loadings such as thermal and electromagnetic fields in recent years. For acoustic wave resonators, such studies are equally important and the wave propagation in FGM structures can be used in the evaluation and assessment of performance, reliability, and life of sensors based on acoustic waves such as the quartz crystal microbalances (QCM). Now we studied the thickness-shear vibrations of FGM plates with properties of AT-cut quartz crystal varying along the thickness in a general pattern represented by a trigonometric function with both sine and cosine functions of the thickness coordinate. The solutions are obtained by using Fourier expansion of the plate deformation. We also obtained the frequency changes of the fundamental and overtone modes which are strongly coupled for the evaluation of resonator structures with property variation or design to take advantages of FGM in novel applications.

preprint2015arXiv

Tuning and Optimizing the Finite Element Analysis with Elements of Large Nodal DOF on a Linux Cluster

The finite element analysis of high frequency vibrations of quartz crystal plates is a necessary process required in the design of quartz crystal resonators of precision types for applications in filters and sensors. The anisotropic materials and extremely high frequency in radiofrequency range of resonators determine that vibration frequency spectra are complicated with strong couplings of large number of different vibration modes representing deformations which do not appear in usual structural problems. For instance, the higher-order thickness-shear vibrations usually representing the sharp deformation of thin plates in the thickness direction, expecting the analysis is to be done with refined meshing schemes along the relatively small thickness and consequently the large plane area. To be able to represent the precise vibration mode shapes, a very large number of elements are needed in the finite element analysis with either the three-dimensional theory or the higher-order plate theory, although considerable reduction of numbers of degree-of-freedom (DOF) are expected for the two-dimensional analysis without scarifying the accuracy. In this paper, we reviewed the software architecture for the analysis and demonstrated the evaluation and tuning of parameters for the improvement of the analysis with problems of elements with a large number of DOF in each node, or a problem with unusually large bandwidth of the banded stiffness and mass matrices in comparison with conventional finite element formulation. Such a problem can be used as an example for the optimization and tuning of problems from multi-physics analysis which are increasingly important in applications with excessive large number of DOF and bandwidth in engineering.

preprint2014arXiv

Accurate Atmospheric Parameters at Moderate Resolution Using Spectral Indices: Preliminary Application to the MARVELS Survey

Studies of Galactic chemical and dynamical evolution in the solar neighborhood depend on the availability of precise atmospheric parameters (Teff, [Fe/H] and log g) for solar-type stars. Many large-scale spectroscopic surveys operate at low to moderate spectral resolution for efficiency in observing large samples, which makes the stellar characterization difficult due to the high degree of blending of spectral features. While most surveys use spectral synthesis, in this work we employ an alternative method based on spectral indices to determine the atmospheric parameters of a sample of nearby FGK dwarfs and subgiants observed by the MARVELS survey at moderate resolving power (R~12,000). We have developed three codes to automatically normalize the observed spectra, measure the equivalent widths of the indices and, through the comparison of those with values calculated with pre-determined calibrations, derive the atmospheric parameters of the stars. The calibrations were built using a sample of 309 stars with precise stellar parameters obtained from the analysis of high-resolution FEROS spectra. A validation test of the method was conducted with a sample of 30 MARVELS targets that also have reliable atmospheric parameters from high-resolution spectroscopic analysis. Our approach was able to recover the parameters within 80 K for Teff, 0.05 dex for [Fe/H] and 0.15 dex for log g, values that are lower or equal to the typical external uncertainties found between different high-resolution analyzes. An additional test was performed with a subsample of 138 stars from the ELODIE stellar library and the literature atmospheric parameters were recovered within 125 K for Teff, 0.10 dex for [Fe/H] and 0.29 dex for log g. These results show that the spectral indices are a competitive tool to characterize stars with the intermediate resolution spectra.

preprint2014arXiv

Influence of Stellar Multiplicity On Planet Formation. I. Evidence of Suppressed Planet Formation Due to Stellar Companions Within 20 AU and Validation of Four Planets From the Kepler Multiple Planet Candidates

The planet occurrence rate for multiple stars is important in two aspects. First, almost half of stellar systems in the solar neighborhood are multiple systems. Second, the comparison of the planet occurrence rate for multiple stars to that for single stars sheds light on the influence of stellar multiplicity on planet formation and evolution. We developed a method of distinguishing planet occurrence rates for single and multiple stars. From a sample of 138 bright (KP <13.5) Kepler multi-planet candidate systems, we compared the stellar multiplicity rate of these planet host stars to that of field stars. Using dynamical stability analyses and archival Doppler measurements, we find that the stellar multiplicity rate of planet host stars is significantly lower than field stars for semi-major axes less than 20 AU, suggesting that planet formation and evolution are suppressed by the presence of a close-in companion star at these separations. The influence of stellar multiplicity at larger separations is uncertain because of search in- completeness due to a limited Doppler observation time baseline and a lack of high resolution imaging observation. We calculated the planet confidence for the sample of mutlti-planet candidates, and find that the planet confidences for KOI 82.01, KOI 115.01, KOI 282.01 and KOI 1781.02 are higher than 99.7% and thus validate the planetary nature of these four planet candidates. This sample of bright Kepler multi-planet candidates with refined stellar and orbital parameters, planet confidence estimation, and nearby stellar companion identification offers a well-characterized sample for future theoretical and observational study.

preprint2014arXiv

Influence of Stellar Multiplicity On Planet Formation. II. Planets Are Less Common in Multiple-Star Systems with Separations Smaller than 1500 AU

Almost half of the stellar systems in the solar neighborhood are made up of multiple stars. In multiple-star systems, planet formation is under the dynamical influence of stellar companions, and the planet occurrence rate is expected to be different from that for single stars. There have been numerous studies on the planet occurrence rate of single star systems. However, to fully understand planet formation, the planet occurrence rate in multiple-star systems needs to be addressed. In this work, we {infer} the planet occurrence rate in multiple-star systems by measuring the stellar multiplicity rate for planet host stars. For a sub-sample of 56 $Kepler$ planet host stars, we use adaptive optics (AO) imaging and the radial velocity (RV) technique to search for stellar companions. The combination of these two techniques results in high search completeness for stellar companions. We detect 59 visual stellar companions to 25 planet host stars with AO data. {{Three stellar companions are within 2$^{\prime\prime}$, and 27 within 6$^{\prime\prime}$. We also detect 2 possible stellar companions (KOI 5 and KOI 69) showing long-term RV acceleration.}} After correcting for a bias against planet detection in multiple-star systems due to flux contamination, we find that planet formation is suppressed in multiple-star systems with separations smaller than 1500 AU. Specifically, we find that compared to single star systems, planets in multiple-star systems occur $4.5\pm3.2$, $2.6\pm1.0$, and $1.7\pm0.5$ times less frequently when a stellar companion is present at a distance of 10, 100, and 1000 AU, respectively. This conclusion applies only to circumstellar planets; the planet occurrence rate for circumbinary planets requires further investigation.

preprint2014arXiv

MPISE: Symbolic Execution of MPI Programs

Message Passing Interfaces (MPI) plays an important role in parallel computing. Many parallel applications are implemented as MPI programs. The existing methods of bug detection for MPI programs have the shortage of providing both input and non-determinism coverage, leading to missed bugs. In this paper, we employ symbolic execution to ensure the input coverage, and propose an on-the-fly schedule algorithm to reduce the interleaving explorations for non-determinism coverage, while ensuring the soundness and completeness. We have implemented our approach as a tool, called MPISE, which can automatically detect the deadlock and runtime bugs in MPI programs. The results of the experiments on benchmark programs and real world MPI programs indicate that MPISE finds bugs effectively and efficiently. In addition, our tool also provides diagnostic information and replay mechanism to help understanding bugs.

preprint2014arXiv

On the Occurrence Rate of Hot Jupiters in Different Stellar Environments

Many Hot Jupiters (HJs) are detected by the Doppler and the transit techniques. From surveys using these two techniques, however, the measured HJ occurrence rates differ by a factor of two or more. Using the California Planet Survey sample and the Kepler sample, we investigate the causes for the difference of HJ occurrence rate. First, we find that $12.8\%\pm0.24\%$ of HJs are misidentified in the Kepler mission because of photometric dilution and subgiant contamination. Second, we explore the differences between the Doppler sample and the Kepler sample that can account for the different HJ occurrence rate. Third, we discuss how to measure the fundamental HJ occurrence rates by synthesizing the results from the Doppler and Kepler surveys. The fundamental HJ occurrence rates are a measure of HJ occurrence rate as a function of stellar multiplicity and evolutionary stage, e.g., the HJ occurrence rate for single and multiple stars or for main sequence and subgiant stars. While we find qualitative evidence that HJs occur less frequently in subgiants and multiple stellar systems, we conclude that our current knowledge of stellar properties and stellar multiplicity rate is too limited for us to reach any quantitative result for the fundamental HJ occurrence rates. This concern extends to $η_{\rm{Earth}}$, the occurrence rate of Earth-like planets.

preprint2014arXiv

Planet Hunters VII. Discovery of a New Low-Mass, Low-Density Planet (PH3 c) Orbiting Kepler-289 with Mass Measurements of Two Additional Planets (PH3 b and d)

We report the discovery of one newly confirmed planet ($P=66.06$ days, $R_{\rm{P}}=2.68\pm0.17R_\oplus$) and mass determinations of two previously validated Kepler planets, Kepler-289 b ($P=34.55$ days, $R_{\rm{P}}=2.15\pm0.10R_\oplus$) and Kepler-289-c ($P=125.85$ days, $R_{\rm{P}}=11.59\pm0.10R_\oplus$), through their transit timing variations (TTVs). We also exclude the possibility that these three planets reside in a $1:2:4$ Laplace resonance. The outer planet has very deep ($\sim1.3%$), high signal-to-noise transits, which puts extremely tight constraints on its host star's stellar properties via Kepler's Third Law. The star PH3 is a young ($\sim1$ Gyr as determined by isochrones and gyrochronology), Sun-like star with $M_*=1.08\pm0.02M_\odot$, $R_*=1.00\pm0.02R_\odot$, and $T_{\rm{eff}}=5990\pm38$ K. The middle planet's large TTV amplitude ($\sim5$ hours) resulted either in non-detections or inaccurate detections in previous searches. A strong chopping signal, a shorter period sinusoid in the TTVs, allows us to break the mass-eccentricity degeneracy and uniquely determine the masses of the inner, middle, and outer planets to be $M=7.3\pm6.8M_\oplus$, $4.0\pm0.9M_\oplus$, and $M=132\pm17M_\oplus$, which we designate PH3 b, c, and d, respectively. Furthermore, the middle planet, PH3 c, has a relatively low density, $ρ=1.2\pm0.3$ g/cm$^3$ for a planet of its mass, requiring a substantial H/He atmosphere of $2.1^{+0.8}_{-0.3}%$ by mass, and joins a growing population of low-mass, low-density planets.

preprint2014arXiv

Planet Hunters. VI: An Independent Characterization of KOI-351 and Several Long Period Planet Candidates from the Kepler Archival Data

We report the discovery of 14 new transiting planet candidates in the Kepler field from the Planet Hunters citizen science program. None of these candidates overlapped with Kepler Objects of Interest (KOIs) at the time of submission. We report the discovery of one more addition to the six planet candidate system around KOI-351, making it the only seven planet candidate system from Kepler. Additionally, KOI-351 bears some resemblance to our own solar system, with the inner five planets ranging from Earth to mini-Neptune radii and the outer planets being gas giants; however, this system is very compact, with all seven planet candidates orbiting $\lesssim 1$ AU from their host star. A Hill stability test and an orbital integration of the system shows that the system is stable. Furthermore, we significantly add to the population of long period transiting planets; periods range from 124-904 days, eight of them more than one Earth year long. Seven of these 14 candidates reside in their host star's habitable zone.

preprint2014arXiv

Revealing A Universal Planet-Metallicity Correlation For Planets of Different Sizes Around Solar-Type Stars

The metallicity of exoplanet systems serves as a critical diagnostic of planet formation mechanisms. Previous studies have demonstrated the planet-metallicity correlation for large planets ($R_P\ \geq\ 4\ R_E$); however, a correlation has not been found for smaller planets. With a sample of 406 $Kepler$ Objects of Interest whose stellar properties are determined spectroscopically, we reveal a universal planet-metallicity correlation: not only gas-giant planets ($3.9\ R_E\ < R_P\ \leq\ 22.0\ R_E$) but also gas-dwarf ($1.7\ R_E\ < R_P\ \leq\ 3.9\ R_E$) and terrestrial planets ($R_P\ \leq\ 1.7\ R_E$) occur more frequently in metal-rich stars. The planet occurrence rates of gas-giant planets, gas-dwarf planets, and terrestrial planets are $9.30^{+5.62}_{-3.04}$, $2.03^{+0.29}_{-0.26}$, and $1.72^{+0.19}_{-0.17}$ times higher for metal-rich stars than for metal-poor stars, respectively.

preprint2014arXiv

Thickness-shear Frequencies of an Infinite Quartz Plate with Material Property Variation Along the Thickness

Properties of the quartz crystal blank of a resonator is assumed homogeneous, uniform, and perfect in design, manufacturing, and applications. As end products, quartz crystal resonators are frequently exposed to gases and liquids which can cause surface damage and internal degradation of blanks under increasingly hostile conditions. The combination of service conditions and manufacturing process including chemical etching and polishing can inevitably modify the surface of quartz crystal blanks with changes of material properties, raising the question of what will happen to vibrations of quartz crystal resonators of thickness-shear type if such modifications to blanks are to be evaluated for sensitive applications. Such questions have been encountered in other materials and structures with property variations either on purpose or as the effect of environmental or natural processes commonly referred to as functionally graded materials, or FGMs. Analyses have been done in applications as part of studies on FGMs in structural as well as in acoustic wave device applications. A procedure based on series solutions has been developed in the evaluation of frequency changes and features in an infinite quartz crystal plate of AT-cut with the symmetric material variation pattern given in a cosine function with the findings that the vibration modes are now closely coupled. These results can be used in the evaluation of surface damage and corrosion of quartz crystal blanks of resonators in sensor applications or development of new structures of resonators.

preprint2013arXiv

An Analysis of Nonlinear Thickness-shear Vibrations of Quartz Crystal Plates by Two-Dimensional Finite Element Method

A nonlinear analysis of high-frequency thickness-shear vibrations of AT-cut quartz crystal plates is presented with the two-dimensional finite element method. We expanded both kinematic and constitutive nonlinear Mindlin plate equations and then truncated them to the first-order equations as an approximation, which is used later for the formulation of nonlinear finite element analysis with all zeroth- and first-order displacements and electric potentials. The matrix equation of motion is established with the first-order harmonic approximation and the generalized nonlinear eigensystem is solved by a direct iterative procedure. A backbone curve and corresponding mode shapes are obtained and analyzed. The nonlinear finite element program is developed based on earlier linear edition and can be utilized to predict nonlinear characteristics of miniaturized quartz crystal resonators in the design process.

preprint2013arXiv

Counterexample-Preserving Reduction for Symbolic Model Checking

The cost of LTL model checking is highly sensitive to the length of the formula under verification. We observe that, under some specific conditions, the input LTL formula can be reduced to an easier-to-handle one before model checking. In our reduction, these two formulae need not to be logically equivalent, but they share the same counterexample set w.r.t the model. In the case that the model is symbolically represented, the condition enabling such reduction can be detected with a lightweight effort (e.g., with SAT-solving). In this paper, we tentatively name such technique "Counterexample-Preserving Reduction" (CePRe for short), and finally the proposed technquie is experimentally evaluated by adapting NuSMV.

preprint2013arXiv

KELT-3b: A Hot Jupiter Transiting a V=9.8 Late-F Star

We report the discovery of KELT-3b, a moderately inflated transiting hot Jupiter with a mass of 1.477 (-0.067, +0.066) M_J, and radius of 1.345 +/- 0.072 R_J, with an orbital period of 2.7033904 +/- 0.000010 days. The host star, KELT-3, is a V=9.8 late F star with M_* = 1.278 (-0.061, +0.063) M_sun, R_* = 1.472 (-0.067, +0.065) R_sun, T_eff = 6306 (-49, +50) K, log(g) = 4.209 (-0.031, +0.033), and [Fe/H] = 0.044 (-0.082, +0.080), and has a likely proper motion companion. KELT-3b is the third transiting exoplanet discovered by the KELT survey, and is orbiting one of the 20 brightest known transiting planet host stars, making it a promising candidate for detailed characterization studies. Although we infer that KELT-3 is significantly evolved, a preliminary analysis of the stellar and orbital evolution of the system suggests that the planet has likely always received a level of incident flux above the empirically-identified threshold for radius inflation suggested by Demory & Seager (2011).

preprint2013arXiv

MARVELS-1: A face-on double-lined binary star masquerading as a resonant planetary system; and consideration of rare false positives in radial velocity planet searches

We have analyzed new and previously published radial velocity observations of MARVELS-1, known to have an ostensibly substellar companion in a ~6- day orbit. We find significant (~100 m/s) residuals to the best-fit model for the companion, and these residuals are naively consistent with an interior giant planet with a P = 1.965d in a nearly perfect 3:1 period commensuribility (|Pb/Pc - 3| < 10^{-4}). We have performed several tests for the reality of such a companion, including a dynamical analysis, a search for photometric variability, and a hunt for contaminating stellar spectra. We find many reasons to be critical of a planetary interpretation, including the fact that most of the three-body dynamical solutions are unstable. We find no evidence for transits, and no evidence of stellar photometric variability. We have discovered two apparent companions to MARVELS-1 with adaptive optics imaging at Keck; both are M dwarfs, one is likely bound, and the other is likely a foreground object. We explore false-alarm scenarios inspired by various curiosities in the data. Ultimately, a line profile and bisector analysis lead us to conclude that the ~100 m/s residuals are an artifact of spectral contamination from a stellar companion contributing ~15-30% of the optical light in the system. We conclude that origin of this contamination is the previously detected radial velocity companion to MARVELS-1, which is not, as previously reported, a brown dwarf, but in fact a G dwarf in a face-on orbit.

preprint2013arXiv

Planet Hunters. V. A Confirmed Jupiter-Size Planet in the Habitable Zone and 42 Planet Candidates from the Kepler Archive Data

We report the latest Planet Hunter results, including PH2 b, a Jupiter-size (R_PL = 10.12 \pm 0.56 R_E) planet orbiting in the habitable zone of a solar-type star. PH2 b was elevated from candidate status when a series of false positive tests yielded a 99.9% confidence level that transit events detected around the star KIC 12735740 had a planetary origin. Planet Hunter volunteers have also discovered 42 new planet candidates in the Kepler public archive data, of which 33 have at least three transits recorded. Most of these transit candidates have orbital periods longer than 100 days and 20 are potentially located in the habitable zones of their host stars. Nine candidates were detected with only two transit events and the prospective periods are longer than 400 days. The photometric models suggest that these objects have radii that range between Neptune to Jupiter. These detections nearly double the number of gas giant planet candidates orbiting at habitable zone distances. We conducted spectroscopic observations for nine of the brighter targets to improve the stellar parameters and we obtained adaptive optics imaging for four of the stars to search for blended background or foreground stars that could confuse our photometric modeling. We present an iterative analysis method to derive the stellar and planet properties and uncertainties by combining the available spectroscopic parameters, stellar evolution models, and transiting light curve parameters, weighted by the measurement errors. Planet Hunters is a citizen science project that crowd-sources the assessment of NASA Kepler light curves. The discovery of these 43 planet candidates demonstrates the success of citizen scientists at identifying planet candidates, even in longer period orbits with only two or three transit events.

preprint2013arXiv

Very Low Mass Stellar and Substellar Companions to Solar-like Stars From MARVELS IV: A Candidate Brown Dwarf or Low-Mass Stellar Companion to HIP 67526

We report the discovery of a candidate brown dwarf or a very low mass stellar companion (MARVELS-5b) to the star HIP 67526 from the Multi-object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The radial velocity curve for this object contains 31 epochs spread over 2.5 years. Our Keplerian fit using a Markov Chain Monte Carlo approach, reveals that the companion has an orbital period of $90.2695^{+0.0188}_{-0.0187}$ days, an eccentricity of $0.4375 \pm 0.0040$ and a semi-amplitude of $2948.14^{+16.65}_{-16.55}$ m s$^{-1}$. Using additional high-resolution spectroscopy, we find the host star has an effective temperature $T_{\rm{eff}}=6004 \pm 34$ K, a surface gravity $\log g$ [cgs] $=4.55 \pm 0.17$ and a metallicity [Fe/H] $=+0.04 \pm 0.06$. The stellar mass and radius determined through the empirical relationship of Torres et al. (2010), yields 1.10$\pm$0.09 $M_{\sun}$ and 0.92$\pm$0.19 $R_{\sun}$. The minimum mass of MARVELS-5b is $65.0 \pm 2.9 M_{Jup}$, indicating that it is likely to be either a brown dwarf or a very low mass star, thus occupying a relatively sparsely-populated region of the mass function of companions to solar-type stars. The distance to this system is 101$\pm$10 pc from the astrometric measurements of Hipparcos. No stellar tertiary is detected in the high-contrast images taken by either FastCam lucky imaging or Keck adaptive optics imaging, ruling out any star with mass greater than 0.2$M_{\sun}$ at a separation larger than 40 AU.

preprint2013arXiv

Very Low Mass Stellar and Substellar Companions to Solar-Like Stars From MARVELS V: A Low Eccentricity Brown Dwarf from the Driest Part of the Desert, MARVELS-6b

We describe the discovery of a likely brown dwarf (BD) companion with a minimum mass of 31.7 +/- 2.0 M_Jup to GSC 03546-01452 from the MARVELS radial velocity survey, which we designate as MARVELS-6b. For reasonable priors, our analysis gives a probability of 72% that MARVELS-6b has a mass below the hydrogen-burning limit of 0.072 M_Sun, and thus it is a high-confidence BD companion. It has a moderately long orbital period of 47.8929 +0.0063/-0.0062 days with a low eccentricty of 0.1442 +0.0078/-0.0073, and a semi-amplitude of 1644 +12/-13 m/s. Moderate resolution spectroscopy of the host star has determined the following parameters: T_eff = 5598 +/- 63, log g = 4.44 +/- 0.17, and [Fe/H] = +0.40 +/- 0.09. Based upon these measurements, GSC 03546-01452 has a probable mass and radius of M_star = 1.11 +/- 0.11 M_Sun and R_star = 1.06 +/- 0.23 R_Sun with an age consistent with less than ~6 Gyr at a distance of 219 +/- 21 pc from the Sun. Although MARVELS-6b is not observed to transit, we cannot definitively rule out a transiting configuration based on our observations. There is a visual companion detected with Lucky Imaging at 7.7 arcsec from the host star, but our analysis shows that it is not bound to this system. The minimum mass of MARVELS-6b exists at the minimum of the mass functions for both stars and planets, making this a rare object even compared to other BDs.

preprint2012arXiv

Accurate Group Delay Measurement for Radial Velocity Instruments Using the Dispersed Fixed Delay Interferometer Method

The dispersed fixed-delay Intereferometer (DFDI) method is attractive for its low cost, compact size, and multiobject capability in precision radial-velocity (RV) measurements. The phase shift of fringes of stellar absorption lines is measured and then converted to an RV shift via an important parameter, phase-to-velocity scale (PV scale), determined by the group delay (GD) of a fixed-delay interferometer. Two methods of GD measurement using a DFDI Doppler instrument are presented in this article: (1) GD measurement using white-light combs gen- erated by the fixed-delay interferometer and (2) GD calibration using an RV reference star. These two methods provide adequate precision of GD measurement and calibration, given the current RV precision achieved by a DFDI Doppler instrument. They can potentially be used to measure GD of an fixed-delay interferometer for submeter- precision Doppler measurement with a DFDI instrument. Advantages and limitations of each method are discussed in detail. The two methods can serve as standard procedures of PV-scale calibration for DFDI instruments and cross- checks for each other.

preprint2012arXiv

Speculative Symbolic Execution

Symbolic execution is an effective path oriented and constraint based program analysis technique. Recently, there is a significant development in the research and application of symbolic execution. However, symbolic execution still suffers from the scalability problem in practice, especially when applied to large-scale or very complex programs. In this paper, we propose a new fashion of symbolic execution, named Speculative Symbolic Execution (SSE), to speed up symbolic execution by reducing the invocation times of constraint solver. In SSE, when encountering a branch statement, the search procedure may speculatively explore the branch without regard to the feasibility. Constraint solver is invoked only when the speculated branches are accumulated to a specified number. In addition, we present a key optimization technique that enhances SSE greatly. We have implemented SSE and the optimization technique on Symbolic Pathfinder (SPF). Experimental results on six programs show that, our method can reduce the invocation times of constraint solver by 21% to 49% (with an average of 30%), and save the search time from 23.6% to 43.6% (with an average of 30%).

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

preprint2012arXiv

Very Low-Mass Stellar and Substellar Companions to Solar-Like Stars from MARVELS I: A Low Mass Ratio Stellar Companion to TYC 4110-01037-1 in a 79-day Orbit

TYC 4110-01037-1 has a low-mass stellar companion, whose small mass ratio and short orbital period are atypical amongst solar-like (Teff ~< 6000 K) binary systems. Our analysis of TYC 4110-01037-1 reveals it to be a moderately aged (~<5 Gyr) solar-like star having a mass of 1.07 +/- 0.08 MSun and radius of 0.99 +/- 0.18 RSun. We analyze 32 radial velocity measurements from the SDSS-III MARVELS survey as well as 6 supporting radial velocity measurements from the SARG spectrograph on the 3.6m TNG telescope obtained over a period of ~2 years. The best Keplerian orbital fit parameters were found to have a period of 78.994 +/- 0.012 days, an eccentricity of 0.1095 +/- 0.0023, and a semi-amplitude of 4199 +/- 11 m/s. We determine the minimum companion mass (if sin i = 1) to be 97.7 +/- 5.8 MJup. The system's companion to host star mass ratio, >0.087 +/- 0.003, places it at the lowest end of observed values for short period stellar companions to solar-like (Teff ~< 6000 K) stars. One possible way to create such a system would be if a triple-component stellar multiple broke up into a short period, low q binary during the cluster dispersal phase of its lifetime. A candidate tertiary body has been identified in the system via single-epoch, high contrast imagery. If this object is confirmed to be co-moving, we estimate it would be a dM4 star. We present these results in the context of our larger-scale effort to constrain the statistics of low mass stellar and brown dwarf companions to FGK-type stars via the MARVELS survey.

preprint2012arXiv

Very Low-mass Stellar and Substellar Companions to Solar-like Stars from MARVELS II: A Short-period Companion Orbiting an F Star with Evidence of a Stellar Tertiary And Significant Mutual Inclination

We report the discovery via radial velocity of a short-period (P = 2.430420 \pm 0.000006 days) companion to the F-type main sequence star TYC 2930-00872-1. A long-term trend in the radial velocities indicates the presence of a tertiary stellar companion with $P > 2000$ days. High-resolution spectroscopy of the host star yields T_eff = 6427 +/- 33 K, log(g) = 4.52 +/- 0.14, and [Fe/H]=-0.04 +/- 0.05. These parameters, combined with the broad-band spectral energy distribution and parallax, allow us to infer a mass and radius of the host star of M_1=1.21 +/- 0.08 M_\odot and R_1=1.09_{-0.13}^{+0.15} R_\odot. We are able to exclude transits of the inner companion with high confidence. The host star's spectrum exhibits clear Ca H and K core emission indicating stellar activity, but a lack of photometric variability and small v*sin(I) suggest the primary's spin axis is oriented in a pole-on configuration. The rotational period of the primary from an activity-rotation relation matches the orbital period of the inner companion to within 1.5 σ, suggesting they are tidally locked. If the inner companion's orbital angular momentum vector is aligned with the stellar spin axis, as expected through tidal evolution, then it has a stellar mass of M_2 ~ 0.3-0.4 M_\odot. Direct imaging limits the existence of stellar companions to projected separations < 30 AU. No set of spectral lines and no significant flux contribution to the spectral energy distribution from either companion are detected, which places individual upper mass limits of M < 1.0 M_\odot, provided they are not stellar remnants. If the tertiary is not a stellar remnant, then it likely has a mass of ~0.5-0.6 M_\odot, and its orbit is likely significantly inclined from that of the secondary, suggesting that the Kozai-Lidov mechanism may have driven the dynamical evolution of this system.

preprint2012arXiv

Very Low-mass Stellar and Substellar Companions to Solar-like Stars from Marvels III: A Short-Period Brown Dwarf Candidate Around An Active G0Iv Subgiant

We present an eccentric, short-period brown dwarf candidate orbiting the active, slightly evolved subgiant star TYC 2087-00255-1, which has effective temperature T_eff = 5903+/-42 K, surface gravity log (g) = 4.07+/-0.16 (cgs), and metallicity [Fe/H] = -0.23+/-0.07. This candidate was discovered using data from the first two years of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of Sloan Digital Sky Survey. From our 38 radial velocity measurements spread over a two-year time baseline, we derive a Keplerian orbital fit with semi-amplitude K=3.571+/-0.041 km/s, period P=9.0090+/-0.0004 days, and eccentricity e=0.226+/-0.011. Adopting a mass of 1.16+/-0.11 Msun for the subgiant host star, we infer that the companion has a minimum mass of 40.0+/-2.5 M_Jup. Assuming an edge-on orbit, the semimajor axis is 0.090+/-0.003 AU. The host star is photometrically variable at the \sim1% level with a period of \sim13.16+/-0.01 days, indicating that the host star spin and companion orbit are not synchronized. Through adaptive optics imaging we also found a point source 643+/-10 mas away from TYC 2087-00255-1, which would have a mass of 0.13 Msun if it is physically associated with TYC 2087-00255-1 and has the same age. Future proper motion observation should be able to resolve if this tertiary object is physically associated with TYC 2087-00255-1 and make TYC 2087-00255-1 a triple body system. Core Ca II H and K line emission indicate that the host is chromospherically active, at a level that is consistent with the inferred spin period and measured v_{rot}*sin i, but unusual for a subgiant of this T_eff. This activity could be explained by ongoing tidal spin-up of the host star by the companion.

preprint2011arXiv

Fundamental Performance of a Dispersed Fixed Delay Interferometer In Searching For Planets Around M Dwarfs

We present a new method to calculate fundamental Doppler measurement limits with a dispersed fixed-delay interferometer (DFDI) in the near infrared wavelength region for searching for exoplanets around M dwarfs in the coming decade. It is based on calculating the Q factor, a measure of flux-normalized Doppler sensitivity in the fringing spectra created with DFDI. We calculate the Q factor as a function of spectral resolution R, stellar projected rotational velocity V sini, stellar effective temperature T_eff and optical path difference (OPD) of the interferometer. We also compare the DFDI Q factor to that for the popular cross-dispersed echelle spectrograph method (the direct echelle (DE) method). Given the IR Doppler measurement is likely to be detector-limited for a while, we introduce new merit functions, which is directly related to photon-limited RV uncertainty, to evaluate Doppler performance with the DFDI and DE methods. We find that DFDI has strength in wavelength coverage and multi-object capability over the DE for a limited detector resource. We simulate the performance of the InfraRed Exoplanet Tracker (IRET) based on the DFDI design, being considered for the next generation IR Doppler measurements. The predicted photon-limited RV uncertainty suggests that IRET is capable of detecting Earth-like exoplanets in habitable zone around nearby bright M dwarfs if they exist. A new method is developed to quantitatively estimate the influence of telluric lines on RV uncertainty. Our study shows that photon-limited RV uncertainty can be reached if 99% of the strength of telluric lines can be removed from the measured stellar spectra. At low to moderate levels of telluric line strength removal (50% to 90%), the optimal RV uncertainty is typically a factor of 2-3 times larger than photon-limited RV uncertainty.

preprint2011arXiv

How Close Are We To Detecting Earth-like Planets in the Habitable Zone Using the Radial Velocity Technique?

Discovering an Earth-like exoplanet in habitable zone is an important milestone for astronomers in search of extra-terrestrial life. While the radial velocity (RV) technique remains one the most powerful tools in detecting and characterizing exo-planetary sys- tems, we calculate the uncertainties in precision RV measurements considering stellar spectral quality factors, RV calibration sources, stellar noise and telluric contamination in different observational bandpasses and for different spectral types. We predict the optimal observational bandpass for different spectral types using the RV technique un- der a variety of conditions. We compare the RV signal of an Earth-like planet in the habitable zone (HZ) to the near future state of the art RV precision and attempt to answer the question: How close are we to detecting Earth-like planet in the HZ using the RV technique?

preprint2011arXiv

On the Eccentricity Distribution of Short-Period Single-Planet Systems

We apply standard Markov chain Monte Carlo (MCMC) analysis techniques for 50 short- period, single-planet systems discovered with radial velocity technique. We develop a new method for accessing the significance of a non-zero orbital eccentricity, namely Γ analysis, which combines frequentist bootstrap approach with Bayesian analysis of each simulated data set. We find the eccentricity estimations from Γ analysis are generally consistent with results from both standard MCMC analysis and previous references. The Γ method is particular useful for assessing the significance of small eccentricities. Our results suggest that the current sample size is insufficient to draw robust conclusions about the roles of tidal interaction and perturbations in shaping the eccentricity distribution of short-period single-planet systems. We use a Bayesian population analysis to show that a mixture of analytical distributions is a good approximation of the underlying eccentricity distribution. For short-period planets, we find the most probable values of parameters in the analytical functions given the observed eccentricities. These analytical functions can be used in theoretical investigations or as priors for the eccentricity distribution when analyzing short-period planets. As the measurement precision improves and sample size increases, the method can be applied to more complex parametrizations for the underlying distribution of eccentricity for extrasolar planetary systems.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

preprint2010arXiv

MARVELS-1b: A Short-Period, Brown Dwarf Desert Candidate from the SDSS-III MARVELS Planet Search

We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.