Bounded model checking is among the most efficient techniques for the automatic verification of concurrent programs. However, encoding all possible interleavings often requires a huge and complex formula, which significantly limits the salability. This paper proposes a novel and efficient abstraction refinement method for multi-threaded program verification. Observing that the huge formula is usually dominated by the exact encoding of the scheduling constraint, this paper proposes a \tsc based abstraction refinement method, which avoids the huge and complex encoding of BMC. In addition, to obtain an effective refinement, we have devised two graph-based algorithms over event order graph for counterexample validation and refinement generation, which can always obtain a small yet effective refinement constraint. Enhanced by two constraint-based algorithms for counterexample validation and refinement generation, we have proved that our method is sound and complete w.r.t. the given loop unwinding depth. Experimental results on \svcompc benchmarks indicate that our method is promising and significantly outperforms the existing state-of-the-art tools.

Atomically engineered oxide heterostructures provide a fertile ground for creating novel states. For example, a two-dimensional electron gas at the interface between two oxide insulators, giant thermoelectric Seebeck coefficient, emergent ferromagnetism from otherwise nonmagnetic components, and colossal ionic conductivity. Extensive research efforts reveal that oxygen deficiency or lattice strain play an important role in determining these unexpected properties. Herein, by studying the abrupt presence of robust ferromagnetism (up to 1.5 uB/Mn) in LaMnO3-based heterostructures, we find the multivalence states of Mn that play a decisive role in the emergence of ferromagnetism in the otherwise antiferromagnetic LaMnO3 thin films. Combining spatially resolved electron energy-loss spectroscopy, X-ray absorption spectroscopy and X-ray magnetic circular dichroism techniques, we determine unambiguously that the ferromagnetism results from a conventional Mn3+-O-Mn4+ double-exchange mechanism rather than an interfacial effect. In contrast, the magnetic dead layer of 5 unit cell in proximity to the interface is found to be accompanied with the accumulation of Mn2+ induced by electronic reconstruction. These findings provide a hitherto-unexplored multivalence state of Mn on the emergent magnetism in undoped manganite epitaxial thin films, such as LaMnO3 and BiMnO3, and shed new light on all-oxide spintronic devices.

Message Passing Interface (MPI) is the standard paradigm of programming in high performance computing. MPI programming takes significant effort, and is error-prone. Thus, effective tools for analyzing MPI programs are much needed. On the other hand, analyzing MPI programs itself is challenging because of non-determinism caused by program inputs and non-deterministic operations. Existing approaches for analyzing MPI programs either do not handle inputs or fail to support programs with mixed blocking and non-blocking operations. This paper presents MPI symbolic verifier (MPI-SV), the first symbolic execution based tool for verifying MPI programs having both blocking and non-blocking operations. To ensure soundness, we propose a blockingdriven matching algorithm to safely handle non-deterministic operations, and a method to soundly and completely model the equivalent behavior of a program execution path. The models of MPI program paths are generated on-the-fly during symbolic execution, and verified w.r.t. the expected properties by model checking. To improve scalability, MPI-SV uses the results of model checking to prune redundant paths. We have implemented MPI-SV and evaluated it on the verification of deadlock freedom for 108 real-world MPI tasks. The pure symbolic execution based technique can successfully verify 61 out of the 108 tasks (56%) within one hour, while in comparison, MPI-SV can verify 94 tasks (87%), a 31% improvement. On average, MPI-SV also achieves 7.25X speedup on verifying deadlock freedom and 2.64X speedup on finding deadlocks. These experimental results are promising, and demonstrate MPI-SV's effectiveness and efficiency.

EarthFinder is a Probe Mission concept selected for study by NASA for input to the 2020 astronomy decadal survey. This study is currently active and a final white paper report is due to NASA at the end of calendar 2018. We are tasked with evaluating the scientific rationale for obtaining precise radial velocity (PRV) measurements in space, which is a two-part inquiry: What can be gained from going to space? What can't be done form the ground? These two questions flow down to these specific tasks for our study - Identify the velocity limit, if any, introduced from micro- and macro-telluric absorption in the Earth's atmosphere; Evaluate the unique advantages that a space-based platform provides to emable the identification and mitigation of stellar acitivity for multi-planet signal recovery.

Blockchain technologies are gaining massive momentum in the last few years. Blockchains are distributed ledgers that enable parties who do not fully trust each other to maintain a set of global states. The parties agree on the existence, values and histories of the states. As the technology landscape is expanding rapidly, it is both important and challenging to have a firm grasp of what the core technologies have to offer, especially with respect to their data processing capabilities. In this paper, we first survey the state of the art, focusing on private blockchains (in which parties are authenticated). We analyze both in-production and research systems in four dimensions: distributed ledger, cryptography, consensus protocol and smart contract. We then present BLOCKBENCH, a benchmarking framework for understanding performance of private blockchains against data processing workloads. We conduct a comprehensive evaluation of three major blockchain systems based on BLOCKBENCH, namely Ethereum, Parity and Hyperledger Fabric. The results demonstrate several trade-offs in the design space, as well as big performance gaps between blockchain and database systems. Drawing from design principles of database systems, we discuss several research directions for bringing blockchain performance closer to the realm of databases.

In this paper, we construct a special kind of breather solution of the nonlinear Schr\"{o}dinger (NLS) equation, the so-called breather-positon ({\it b-positon} for short), which can be obtained by taking the limit $\lambda_{j}$ $\rightarrow$ $\lambda_{1}$ of the Lax pair eigenvalues in the order-$n$ periodic solution which is generated by the $n$-fold Darboux transformation from a special "seed" solution--plane wave. Further, an order-$n$ {\it b-positon} gives an order-$n$ rogue wave under a limit $\lambda_1\rightarrow \lambda_0$. Here $\lambda_0$ is a special eigenvalue in a breather of the NLS equation such that its period goes to infinity. Several analytical plots of order-2 breather confirm visually this double degeneration. The last limit in this double degeneration can be realized approximately in an optical fiber governed by the NLS equation, in which an injected initial ideal pulse is created by a frequency comb system and a programable optical filter (wave shaper) according to the profile of an analytical form of the {\it b-positon} at a certain position $z_0$. We also suggest a new way to observe higher-order rogue waves generation in an optical fiber, namely, measure the patterns at the central region of the higher-order {\it b-positon} generated by above ideal initial pulses when $\lambda_1$ is very close to the $\lambda_0$. The excellent agreement between the numerical solutions generated from initial ideal inputs with a low signal noise ratio and analytical solutions of order-2 {\it b-positon}, supports strongly this way in a realistic optical fiber system. Our results also show the validity of the generating mechanism of a higher-order rogue waves from a multi-breathers through the double degeneration.

Direct imaging of exoplanets presents a formidable technical challenge owing to the small angular separation and high contrast between exoplanets and their host stars. High Dispersion Coronagraphy (HDC) is a pathway to achieve unprecedented sensitivity to Earth-like planets in the habitable zone. Here, we present a framework to simulate HDC observations and data analyses. The goal of these simulations is to perform a detailed analysis of the trade-off between raw star light suppression and spectral resolution for various instrument configurations, target types, and science cases. We predict the performance of an HDC instrument at Keck observatory for characterizing directly imaged gas-giant planets in near infrared bands. We also simulate HDC observations of an Earth-like planet using next-generation ground-based (TMT) and spaced-base telescopes (HabEx and LUVOIR). We conclude that ground-based ELTs are more suitable for HDC observations of an Earth-like planet than future space-based missions owing to the considerable difference in collecting area. For ground-based telescopes, HDC observations can detect an Earth-like planet in the habitable zone around an M dwarf star at 10$^{-4}$ starlight suppression level. Compared to the 10$^{-7}$ planet/star contrast, HDC relaxes the starlight suppression requirement by a factor of 10$^3$. For space-based telescopes, detector noise will be a major limitation at spectral resolutions higher than 10$^4$. Considering detector noise and speckle chromatic noise, R=400 (1600) is the optimal spectral resolutions for HabEx(LUVOIR). The corresponding starlight suppression requirement to detect a planet with planet/star contrast=$6.1\times10^{-11}$ is relaxed by a factor of 10 (100) for HabEx (LUVOIR).

The height of an $n$th-order fundamental rogue wave $q_{\rm rw}^{[n]}$ for the nonlinear Schr\"odinger equation, namely $(2n+1)c$, is proved directly by a series of row operations on matrices appeared in the $n$-fold Darboux transformation. Here the positive constant $c$ denotes the height of the asymptotical plane of the rogue wave.

Blockchain technologies are taking the world by storm. Public blockchains, such as Bitcoin and Ethereum, enable secure peer-to-peer applications like crypto-currency or smart contracts. Their security and performance are well studied. This paper concerns recent private blockchain systems designed with stronger security (trust) assumption and performance requirement. These systems target and aim to disrupt applications which have so far been implemented on top of database systems, for example banking, finance applications. Multiple platforms for private blockchains are being actively developed and fine tuned. However, there is a clear lack of a systematic framework with which different systems can be analyzed and compared against each other. Such a framework can be used to assess blockchains' viability as another distributed data processing platform, while helping developers to identify bottlenecks and accordingly improve their platforms. In this paper, we first describe BlockBench, the first evaluation framework for analyzing private blockchains. It serves as a fair means of comparison for different platforms and enables deeper understanding of different system design choices. Any private blockchain can be integrated to BlockBench via simple APIs and benchmarked against workloads that are based on real and synthetic smart contracts. BlockBench measures overall and component-wise performance in terms of throughput, latency, scalability and fault-tolerance. Next, we use BlockBench to conduct comprehensive evaluation of three major private blockchains: Ethereum, Parity and Hyperledger Fabric. The results demonstrate that these systems are still far from displacing current database systems in traditional data processing workloads. Furthermore, there are gaps in performance among the three systems which are attributed to the design choices at different layers of the software stack.

High-dispersion coronagraphy (HDC) optimally combines high contrast imaging techniques such as adaptive optics/wavefront control plus coronagraphy to high spectral resolution spectroscopy. HDC is a critical pathway towards fully characterizing exoplanet atmospheres across a broad range of masses from giant gaseous planets down to Earth-like planets. In addition to determining the molecular composition of exoplanet atmospheres, HDC also enables Doppler mapping of atmosphere inhomogeneities (temperature, clouds, wind), as well as precise measurements of exoplanet rotational velocities. Here, we demonstrate an innovative concept for injecting the directly-imaged planet light into a single-mode fiber, linking a high-contrast adaptively-corrected coronagraph to a high-resolution spectrograph (diffraction-limited or not). Our laboratory demonstration includes three key milestones: close-to-theoretical injection efficiency, accurate pointing and tracking, on-fiber coherent modulation and speckle nulling of spurious starlight signal coupling into the fiber. Using the extreme modal selectivity of single-mode fibers, we also demonstrated speckle suppression gains that outperform conventional image-based speckle nulling by at least two orders of magnitude.

Many brown dwarfs exhibit photometric variability at levels from tenths to tens of percents. The photometric variability is related to magnetic activity or patchy cloud coverage, characteristic of brown dwarfs near the L-T transition. Time-resolved spectral monitoring of brown dwarfs provides diagnostics of cloud distribution and condensate properties. However, current time-resolved spectral studies of brown dwarfs are limited to low spectral resolution (R$\sim$100) with the exception of the study of Luhman 16 AB at resolution of 100,000 using the VLT$+$CRIRES. This work yielded the first map of brown dwarf surface inhomogeneity, highlighting the importance and unique contribution of high spectral resolution observations. Here, we report on the time-resolved high spectral resolution observations of a nearby brown dwarf binary, 2MASSW J0746425+200032AB. We find no coherent spectral variability that is modulated with rotation. Based on simulations we conclude that the coverage of a single spot on 2MASSW J0746425+200032AB is smaller than 1\% or 6.25\% if spot contrast is 50\% or 80\% of its surrounding flux, respectively. Future high spectral resolution observations aided by adaptive optics systems can put tighter constraints on the spectral variability of 2MASSW J0746425+200032AB and other nearby brown dwarfs.

Today's storage systems expose abstractions which are either too low-level (e.g., key-value store, raw-block store) that they require developers to re-invent the wheels, or too high-level (e.g., relational databases, Git) that they lack generality to support many classes of applications. In this work, we propose and implement a general distributed data storage system, called UStore, which has rich semantics. UStore delivers three key properties, namely immutability, sharing and security, which unify and add values to many classes of today's applications, and which also open the door for new applications. By keeping the core properties within the storage, UStore helps reduce application development efforts while offering high performance at hand. The storage embraces current hardware trends as key enablers. It is built around a data-structure similar to that of Git, a popular source code versioning system, but it also synthesizes many designs from distributed systems and databases. Our current implementation of UStore has better performance than general in-memory key-value storage systems, especially for version scan operations. We port and evaluate four applications on top of UStore: a Git-like application, a collaborative data science application, a transaction management application, and a blockchain application. We demonstrate that UStore enables faster development and the UStore-backed applications can have better performance than the existing implementations.

Broad emission-line outflows of active galactic nuclei (AGNs) have been proposed for many years but are very difficult to quantitatively study because of the coexistence of the gravitationally-bound and outflow emission. We present detailed analysis of a heavily reddened quasar, SDSS J000610.67+121501.2, whose normal ultraviolet (UV) broad emission lines (BELs) are heavily suppressed by the Dusty Torus as a natural "Coronagraph", thus the blueshifted BELs (BBELs) can be reliably measured. The physical properties of the emission-line outflows are derived as follows: ionization parameter $U \sim 10^{-0.5}$, column density $N_{\rm H}\sim 10^{22.0}$ cm$^{-2}$, covering fraction of $\sim 0.1$ and upper limit density of $n_{\rm H}\sim 10^{5.8}$ cm$^{-3}$. The outflow gases are located at least 41 pc away from the central engine, which suggests that they have expanded to the scale of the dust torus or beyond. Besides, Lya shows a narrow symmetric component, to our surprise, which is undetected in any other lines. After inspecting the narrow emission-line region and the starforming region as the origin of the Lya narrow line, we propose the end-result of outflows, diffusing gases in the larger region, acts as the screen of Lya photons. Future high spatial resolution spectrometry and/or spectropolarimetric observation are needed to make a final clarification.

We find transient, transit-like dimming events within the K2 time series photometry of the young star RIK-210 in the Upper Scorpius OB association. These dimming events are variable in depth, duration, and morphology. High spatial resolution imaging revealed the star is single, and radial velocity monitoring indicated that the dimming events can not be due to an eclipsing stellar or brown dwarf companion. Archival and follow-up photometry suggest the dimming events are transient in nature. The variable morphology of the dimming events suggests they are not due to a single, spherical body. The ingress of each dimming event is always shallower than egress, as one would expect for an orbiting body with a leading tail. The dimming events are periodic and synchronous with the stellar rotation. However, we argue it is unlikely the dimming events could be attributed to anything on the stellar surface based on the observed depths and durations. Variable obscuration by a protoplanetary disk is unlikely on the basis that the star is not actively accreting and lacks the infrared excess associated with an inner disk. Rather, we explore the possibilities that the dimming events are due to magnetospheric clouds, a transiting protoplanet surrounded by circumplanetary dust and debris, eccentric orbiting bodies undergoing periodic tidal disruption, or an extended field of dust or debris near the corotation radius.

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.

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.

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\sigma$. 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\sigma$ 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.

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 $\Delta m_{I}=4.4$~mag and $\Delta 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.

$\mu$ 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.

The census of exoplanets is incomplete for orbital distances larger than 1 AU. Here, we present 41 long-period planet candidates in 38 systems identified by Planet Hunters based on Kepler archival data (Q0-Q17). Among them, 17 exhibit only one transit, 14 have two visible transits and 10 have more than three visible transits. For planet candidates with only one visible transit, we estimate their orbital periods based on transit duration and host star properties. The majority of the planet candidates in this work (75%) have orbital periods that correspond to distances of 1-3 AU from their host stars. We conduct follow-up imaging and spectroscopic observations to validate and characterize planet host stars. In total, we obtain adaptive optics images for 33 stars to search for possible blending sources. Six stars have stellar companions within 4". We obtain high-resolution spectra for 6 stars to determine their physical properties. Stellar properties for other stars are obtained from the NASA Exoplanet Archive and the Kepler Stellar Catalog by Huber et al. (2014). We validate 7 planet candidates that have planet confidence over 0.997 (3-{\sigma} level). These validated planets include 3 single-transit planets (KIC-3558849b, KIC-5951458b, and KIC-8540376c), 3 planets with double transits (KIC-8540376b, KIC-9663113b, and KIC-10525077b), and 1 planet with 4 transits (KIC-5437945b). This work provides assessment regarding the existence of planets at wide separations and the associated false positive rate for transiting observation (17%-33%). More than half of the long-period planets with at least three transits in this paper exhibit transit timing variations up to 41 hours, which suggest additional components that dynamically interact with the transiting planet candidates. The nature of these components can be determined by follow-up radial velocity and transit observations.

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\'eel 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.

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.

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{\sigma} 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%).

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.

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.