Experimental Limits on Solar Reflected Dark Matter with a New Approach on Accelerated-Dark-Matter-Electron Analysis in Semiconductors

Recently a dark matter-electron (DM-electron) paradigm has drawn much attention. Models beyond the standard halo model describing DM accelerated by high energy celestial bodies are under intense examination as well. In this Letter, a velocity components analysis (VCA) method dedicated to swift analysis of accelerated DM-electron interactions via semiconductor detectors is proposed and the first HPGe detector-based accelerated DM-electron analysis is realized. Utilizing the method, the first germanium based constraint on sub-GeV solar reflected DM-electron interaction is presented with the 205.4  kg·day dataset from the CDEX-10 experiment. In the heavy mediator scenario, our result excels in the mass range of 5-15  keV/c^{2}, achieving a 3 orders of magnitude improvement comparing with previous semiconductor experiments. In the light mediator scenario, the strongest laboratory constraint for DM lighter than 0.1  MeV/c^{2} is presented. The result proves the feasibility and demonstrates the vast potential of the VCA technique in future accelerated DM-electron analyses with semiconductor detectors.

MARS Plus: An Improved Molecular Design Tool for Complex Compounds Involving Ionic, Stereo, and Cis-Trans Isomeric Structures

MARS (Molecular Assembling and Representation Suite) (Hsu et al. J. Chem. Inf. Model. 2019, 59, 3703-3713) is a toolbox for the molecular design of organic molecules. MARS uses integer arrays to represent the elements and connectivity between elements of a molecule. It provides a collection of operations to manipulate the elemental composition and connectivity of a molecule (or a pair of molecules), enabling the creation of novel chemical compounds. In this work, the original MARS is extended to handle complex molecular structures, including geometric (cis-trans) isomers, stereo isomers, cyclic compounds, and ionic species. The extended version of MARS, referred to as MARS+, has a more comprehensive coverage of the chemical space and therefore can explore molecules with a greater chemical and physical diversity. Compared to other molecular design tools, MARS+ is designed to perform all possible manipulations on a given molecule or a pair of molecules. Molecular structure manipulation can be conducted in either a controlled or a random fashion. Furthermore, every structure manipulation has a counterpart so that the operation can be reversed. Nearly any possible chemical structure can be generated with MARS+ via a combination of molecular operations. The capabilities of MARS+ are examined by the design of new ionic liquids (ILs). The results show that MARS+ is a useful tool for computer-aided molecular design (CAMD) and molecular structure enumeration.

Entropy of different phases formed by soft rods

The computation of entropy in liquids and liquid crystal (LC) phases is a big challenge in statistical physics. In this work, we extend the two-phase thermodynamic model (2PT) to shape anisotropic soft repulsive spherocylinders (SRSs) and report the absolute values of entropy for different LC phases for a range of aspect ratios L/D = 2 - 5. We calculate the density of states for different LC phases and decompose it into contributions arising from translational and rotational degrees of freedom. The translational and rotational modes are further partitioned into diffusive, gas-like, and non-diffusive, solid-like components using a fluidicity factor. In the dilute limit, the entropy values obtained from the 2PT method match exactly those of an ideal rigid rotor. We find that, for a given packing fraction, the magnitude of the total entropy is roughly equal regardless of the different LC phases associated with different aspect ratios. We also compute the excess entropy (for L/D = 5) and compare those with the values obtained using the standard integration approach of MD or Monte Carlo equation of state of SRS. The values obtained using both approaches match very well. The rotational and translational fluidicity factors are further used to determine the phase boundaries of different LC phases.

Mechanistic insight into the structure, thermodynamics and dynamics of equilibrium gels of multi-armed DNA nanostars

The unique sequence specificity rule of DNA makes it an ideal molecular building block for constructing periodic arrays and devices with nanoscale accuracy and precision. Here, we present the self-assembly of DNA nanostars having three, four and five arms into a gel phase using a simplistic coarse-grained bead-spring model developed by Z. Xing, C. Ness, D. Frenkel and E. Eiser (Macromolecules, 2019, 52, 504-512). Our simulations show that the DNA nanostars form a thermodynamically stable fully bonded gel phase from an unstructured liquid phase with the lowering of temperature. We characterize the phase transition by calculating several structural features such as the radial distribution function and structure factor. The thermodynamics of gelation is quantified by the potential energy and translational pair-entropy of the system. The phase transition from an arrested gel phase to an unstructured liquid phase has been modelled using a two-state theoretical model. We find that this transition is enthalpy driven, and loss of configuration and translational entropy is counterpoised by enthalpic interaction of the DNA sticky-ends, which gives rise to a gel phase at low temperature. The absolute rotational and translational entropy of the systems, measured using a two-phase thermodynamic model, also substantiates the gel transition. The slowing down of the dynamics upon approaching the transition temperature from a high temperature demonstrates the phase transition to a gel phase. A detailed numerical simulation study of the morphology, dynamics and thermodynamics of DNA gelation can provide guidance for future experiments, is easily extensible to other polymeric systems, and is expected to help in understanding the physics of self-assembly.

Phase behavior of active and passive dumbbells

We report phase separation in a mixture of "hot" and "cold" three-dimensional dumbbells which interact by Lennard-Jones potential. We also have studied the effect of asymmetry of dumbbells and the variation of ratio of "hot" and "cold" dumbbells on their phase separation. The ratio of the temperature difference between hot and cold dumbbells to the temperature of cold dumbbells is a measure of the activity χ of the system. From constant density simulations of symmetric dumbbells, we observe that the "hot" and "cold" dumbbells phase separate at higher activity ratio (χ>5.80) compared to that of a mixture of hot and cold Lennard-Jones monomers (χ>3.44). We find that, in the phase-separated system, the hot dumbbells have high effective volume and hence high entropy which is calculated by two-phase thermodynamic method. The high kinetic pressure of hot dumbbells forces the cold dumbbells to form dense clusters such that at the interface the high kinetic pressure of hot dumbbells is balanced by the virial pressure of cold dumbbells. We find that phase separation pushes the cluster of cold dumbbells to have solidlike ordering. Bond orientation order parameters reveal that the cold dumbbells form solidlike ordering consisting of predominantly face-centered cubic and hexagonal-close packing packing, but the individual dumbbells have random orientations. The simulation of the nonequilibrium system of symmetric dumbbells at different ratios of number of hot dumbbells to cold dumbbells reveals that the critical activity of phase separation decreases with increase in fraction of hot dumbbells. The simulation of equal mixture of hot and cold asymmetric dumbbells revealed that the critical activity of phase separation was independent of the asymmetry of dumbbells. We also observed that the clusters of cold asymmetric dumbbells showed both crystalline and noncrystalline order depending on the asymmetry of dumbbells.

Constraints on Sub-GeV Dark Matter-Electron Scattering from the CDEX-10 Experiment

We present improved germanium-based constraints on sub-GeV dark matter via dark matter-electron (χ-e) scattering using the 205.4  kg·day dataset from the CDEX-10 experiment. Using a novel calculation technique, we attain predicted χ-e scattering spectra observable in high-purity germanium detectors. In the heavy mediator scenario, our results achieve 3 orders of magnitude of improvement for m_{χ} larger than 80  MeV/c^{2} compared to previous germanium-based χ-e results. We also present the most stringent χ-e cross-section limit to date among experiments using solid-state detectors for m_{χ} larger than 90  MeV/c^{2} with heavy mediators and m_{χ} larger than 100  MeV/c^{2} with electric dipole coupling. The result proves the feasibility and demonstrates the vast potential of a new χ-e detection method with high-purity germanium detectors in ultralow radioactive background.

Exotic Dark Matter Search with the CDEX-10 Experiment at China's Jinping Underground Laboratory

A search for exotic dark matter (DM) in the sub-GeV mass range has been conducted using 205 kg day data taken from a p-type point contact germanium detector of the CDEX-10 experiment at China's Jinping underground laboratory. New low-mass dark matter searching channels, neutral current fermionic DM absorption (χ+A→ν+A) and DM-nucleus 3→2 scattering (χ+χ+A→ϕ+A), have been analyzed with an energy threshold of 160 eVee. No significant signal was found; thus new limits on the DM-nucleon interaction cross section are set for both models at the sub-GeV DM mass region. A cross section limit for the fermionic DM absorption is set to be 2.5×10^{-46}  cm^{2} (90% C.L.) at DM mass of 10  MeV/c^{2}. For the DM-nucleus 3→2 scattering scenario, limits are extended to DM mass of 5 and 14  MeV/c^{2} for the massless dark photon and bound DM final state, respectively.

Promotion mechanism for the growth of CO2 hydrate with urea using molecular dynamics simulations

The role of urea as a kinetic promoter for the growth of CO2 hydrates is revealed for the first time using molecular dynamics simulations. Analysis of simulation trajectories shows that urea plays two important roles in the growth process: increasing mass transport of CO2 and catalyzing cage formation at the solid-liquid interface.

Effect of Small Cage Guests on Dissociation Properties of Tetrahydrofuran Hydrates

It is well understood that tetrahydrofuran (THF) molecules are able to stabilize the large cages (5¹²6⁴) of structure II to form the THF hydrate with empty small cages even at atmospheric pressure. This leaves the small cages to store gas molecules at relatively lower pressures and higher temperatures. The dissociation enthalpy and temperature strongly depend on the size of gas molecules enclathrated in the small cages of structure II THF hydrate. A high-pressure microdifferential scanning calorimeter was applied to measure the dissociation enthalpies and temperatures of THF hydrates pressurized by helium and methane under a constant pressure ranging from 0.10 to 35.00 MPa and a wide THF concentration ranging from 0.25 to 8.11 mol %. The dissociation temperature of binary He + THF and methane + THF hydrates increases along with an increase in the THF concentration in the liquid phase at a fixed pressure (e.g., 30 MPa), reaching a maximum value of 280.8 and 312.8 K, respectively, at stoichiometric concentration (5.56 mol % THF), and then remains nearly constant for even higher THF concentrations (>5.56 mol %). The effect of gas occupancy in the small cages on the dissociation enthalpy of He + THF and methane + THF mixed hydrates was further examined by using molecular dynamics (MD) simulations. The dissociation enthalpy of the He-THF mixed hydrates is independent of pressure with an average of 5.68 kJ/mol H₂O over the pressure ranging from 0.10 to 30.0 MPa, consistent with the MD results of the He-THF mixed hydrates with low single occupancy (<23%) of helium molecules in the small cages. Consequently, the heat of adsorption of helium molecules in the small cages of the He-THF mixed hydrates is rather too weak to be identified. On the other hand, the dissociation enthalpy of the methane-THF hydrates increases from 9.11 to 10.01 kJ/mol H₂O along with an increase in methane pressure over the pressure ranging from 5.0 to 30.0 MPa, consistent with the MD results of the methane-THF mixed hydrates with full occupancy of methane molecules in the small cages. These findings provide important information for the design of a potential medium of gas storage and transportation.

A Benchmark Open-Source Implementation of COSMO-SAC

The COSMO-SAC modeling approach has found wide application in science as well as in a range of industries due to its good predictive capabilities. While other models for liquid phases, as for example UNIFAC, are in general more accurate than COSMO-SAC, these models typically contain many adjustable parameters and can be limited in their applicability. In contrast, the COSMO-SAC model only contains a few universal parameters and subdivides the molecular surface area into charged segments that interact with each other. In recent years, additional improvements to the construction of the sigma profiles and evaluation of activity coefficients have been made. In this work, we present a comprehensive description of how to postprocess the results of a COSMO calculation through to the evaluation of thermodynamic properties. We also assembled a large database of COSMO files, consisting of 2261 compounds, freely available to academic and noncommercial users. We especially focus on the documentation of the implementation and provide the optimized source code in C++, wrappers in Python, and sample sigma profiles calculated from each approach, as well as tests and validation results. The misunderstandings in the literature relating to COSMO-SAC are described and corrected. The computational efficiency of the implementation is demonstrated.

Direct Detection Constraints on Dark Photons with the CDEX-10 Experiment at the China Jinping Underground Laboratory

We report constraints on the dark photon effective kinetic mixing parameter (κ) with data taken from two p-type point-contact germanium detectors of the CDEX-10 experiment at the China Jinping Underground Laboratory. The 90% confidence level upper limits on κ of solar dark photon from 205.4 kg-day exposure are derived, probing new parameter space with masses (m_{V}) from 10 to 300  eV/c^{2} in direct detection experiments. Considering dark photon as the cosmological dark matter, limits at 90% confidence level with m_{V} from 0.1 to 4.0  keV/c^{2} are set from 449.6 kg-day data, with a minimum of κ=1.3×10^{-15} at m_{V}=200  eV/c^{2}.

Search for Light Weakly-Interacting-Massive-Particle Dark Matter by Annual Modulation Analysis with a Point-Contact Germanium Detector at the China Jinping Underground Laboratory

We present results on light weakly interacting massive particle (WIMP) searches with annual modulation (AM) analysis on data from a 1-kg mass p-type point-contact germanium detector of the CDEX-1B experiment at the China Jinping Underground Laboratory. Datasets with a total live time of 3.2 yr within a 4.2-yr span are analyzed with analysis threshold of 250 eVee. Limits on WIMP-nucleus (χ-N) spin-independent cross sections as function of WIMP mass (m_{χ}) at 90% confidence level (C.L.) are derived using the dark matter halo model. Within the context of the standard halo model, the 90% C.L. allowed regions implied by the DAMA/LIBRA and CoGeNT AM-based analysis are excluded at >99.99% and 98% C.L., respectively. These results correspond to the best sensitivity at m_{χ}<6  GeV/c^{2} among WIMP AM measurements to date.

Constraints on Spin-Independent Nucleus Scattering with sub-GeV Weakly Interacting Massive Particle Dark Matter from the CDEX-1B Experiment at the China Jinping Underground Laboratory

We report results on the searches of weakly interacting massive particles (WIMPs) with sub-GeV masses (m_{χ}) via WIMP-nucleus spin-independent scattering with Migdal effect incorporated. Analysis on time-integrated (TI) and annual modulation (AM) effects on CDEX-1B data are performed, with 737.1 kg day exposure and 160 eVee threshold for TI analysis, and 1107.5 kg day exposure and 250 eVee threshold for AM analysis. The sensitive windows in m_{χ} are expanded by an order of magnitude to lower DM masses with Migdal effect incorporated. New limits on σ_{χN}^{SI} at 90% confidence level are derived as 2×10^{-32}∼7×10^{-35}  cm^{2} for TI analysis at m_{χ}∼50-180  MeV/c^{2}, and 3×10^{-32}∼9×10^{-38}  cm^{2} for AM analysis at m_{χ}∼75  MeV/c^{2}-3.0  GeV/c^{2}.

New Data Structure for Computational Molecular Design with Atomic or Fragment Resolution

A new molecular data structure and molecular structure operation algorithms are proposed for general purpose molecular design. The data structure allows for a variety of molecular operations for creating new molecules. Two types of molecular operations were developed, unimolecular and bimolecular operations. In unimolecular operations, a child molecule can be created from a parent via addition of a functional group, deletion of a fragment, mutation of an atom, etc. In bimolecular operations, children molecules are generated from two parent molecules through combination or crossover (hybridization). These molecular operations are essential for the creation and modification of molecules for the purpose of molecular design. The data structure is capable of representing linear, branched, multifunctional, and multivalent compounds. Algorithms are developed for deriving the molecular data structure of a molecule from its atomic coordinates and vice versa. We show that this new molecular data structure and the developed algorithms, referred to as Molecular Assembling and Representation Suite, allow one to generate a comprehensive library of new molecules via performing every possible molecular structure modification.

Stable, color-tunable 2D SCN-based perovskites: revealing the critical influence of an asymmetric pseudo-halide on constituent ions

Two-dimensional (2D) layered perovskites (An+1BnX3n+1, n = 1, 2, …) have recently attracted significant research interest because of their enhanced ambient stability compared to their conventional 3D counterparts. In addition to the common A-site cation engineering, using an asymmetric pseudo-halide anion, SCN-, in the anion X-site has been recently proven to be another effective approach to constitute 2D perovskites. Among these, 2D (MA)2Pb(SCN)2I2 was the most widely investigated and was considered to be a promising material owing to its good optoelectronic properties; however, its poor stability has aroused concerns in recent researches. In this study, systematical composition engineering of A2Pb(SCN)2X2 (A = FA+, MA+, Cs+ and X = Br-, I-) was conducted. Our results revealed that the linear SCN- anion dictates critical restrictions on the constituent ions of its derived 2D framework (PbX4(SCN)2), which has not yet been extensively discussed. We demonstrated that using a smaller Cs+ cation can afford a more favorable 2D structure compared with the MA+ cation. Cs2Pb(SCN)2I2 was revealed to possess improved stability and photo-response compared to (MA)2Pb(SCN)2I2. Interestingly, Cs2Pb(SCN)2I2 and (MA)2Pb(SCN)2I2 appear to possess distinct electronic band structures. This is indicated by their discrepant photoluminescence spectra, in which the former exhibits a rather intense singlet emission at room temperature in contrast with the latter, which shows a dominant emission associated with triplet or defective states. Furthermore, using a smaller Cs+ cation enables facile replacement of a smaller halide anion. A series of mixed-halide 2D Cs2Pb(SCN)2(I1-xBrx)2 (x = 0, 1/3, 1/2, 2/3, 1) with varying vivid colors was explored by both calculation and experimental efforts to corroborate the enhanced stability when the x value increases. The results revealed in this study might represent a novel discovery of an inherent trait of the 2D SCN-based perovskites and also suggest that the all-inorganic 2D Cs2Pb(SCN)2X2 perovskite system is a promising class of materials with good stability and color-tunability that deserves further exploration.

Limits on Light Weakly Interacting Massive Particles from the First 102.8 kg×day Data of the CDEX-10 Experiment

We report the first results of a light weakly interacting massive particles (WIMPs) search from the CDEX-10 experiment with a 10 kg germanium detector array immersed in liquid nitrogen at the China Jinping Underground Laboratory with a physics data size of 102.8 kg day. At an analysis threshold of 160 eVee, improved limits of 8×10^{-42} and 3×10^{-36}  cm^{2} at a 90% confidence level on spin-independent and spin-dependent WIMP-nucleon cross sections, respectively, at a WIMP mass (m_{χ}) of 5  GeV/c^{2} are achieved. The lower reach of m_{χ} is extended to 2  GeV/c^{2}.

[A prospective pilot study of combined intra-operative radiotherapy for centrally located hepatocellular carcinomas]

Objective: To carry out a prospective cohort study of combined intra-operative radiotherapy for centrally located hepatocellular carcinomas (HCC) and to observe the safety and postoperative complications. Methods: A total of 79 patients with centrally located HCC who underwent hepatectomy were divided into two groups: experimental group (combined with targeted intra-operative radiotherapy, 32 cases) and control group (single surgical operation, 47 cases). Patients in the experimental group received intra-operative electron radiotherapy after tumor resection, while patients in the control group received to intra-operative electron radiotherapy.The haemorrhagia amount and operation time during the operation, intra-operative liver function and the recovery of liver and gastrointestinal tract of patients in these two groups were compared. Results: No postoperative 30-day mortality was observed in all of the patients. The average total operation time of patients in the experimental group was (319±76) min, significantly longer than (233±76) min of the control group (P<0.001). The average aspartate transaminase (AST) level of patients in the experimental group at postoperative day 1 was 562.5 U/L, significantly higher than 347.0 U/L of control group (P=0.031). However, the average prothrombin activity levels of patients in the experimental group at postoperative day 3 and day 7 were (68.3±17.9)% and (73.4±10.2)%, respectively, significantly lower than (78.9±15.9)% and (80.0±10.6)% of control group (both P<0.05). There were no significant differences of tumor volume, differentiation degree, satellite lesion, dorsal membrane invasion, microvascular invasion between these two groups (all P>0.05). There were no significant differences of hospital stay, ventilation time, the incidence of hepatic insufficiency, ascites, pleural effusion, infection, biliary fistula between these two groups (all P>0.05). There were no significant differences of alanine aminotransferase (ALT), albumin, total bilirubin between these two groups at postoperative day 1, 3, 5 and 7 (all of P>0.05). Conclusion: The resection of centrally located HCC combined with intra-operative radiotherapy may increase the total operation time, delay the early postoperative recovery of liver function, but it is still safe and feasible. Trial registration: National Cancer Centre /Cancer Hospital, Chinese Academy of Medical Sciences, ChiCTR-TRC-12002802.

Study of inactive layer uniformity and charge collection efficiency of a p-type point-contact germanium detector

The characteristics of the surface inactive layer of a 1-kg-mass p-type point-contact germanium detector were studied. The thickness of the inactive layer and its uniformity on the top and lateral surfaces were measured. A charge collection efficiency function was developed according to the Monte Carlo simulation to describe the charge collection capacity along the depth within this inactive layer. In the energy range below 18keV, the surface, bulk, and total spectra of ⁵⁷Co, ¹³³Ba, ¹³⁷Cs, and ⁶⁰Co from simulations based on the charge collection efficiency function were well consistent with those from experiments.

[Analysis of clinicopathological features and prognosis of 98 cases of small hepatocellular carcinoma]

Objective: To analyze the clinicopathological features and prognosis of patients with small hepatocellular carcinoma. Methods: The clinicopathological and follow-up data of 98 patients with small hepatocellular carcinoma who underwent R0 resection from January 2009 to December 2013 were analyzed retrospectively. Results: All of the patients were followed up. Their postoperative 1-year, 3-year and 5-year overall survival rates were 99.0%, 91.7%, and 76.3%, respectively. Their postoperative median overall survival (OS) period was 52 months. The postoperative progression-free survival rates were 86.7%, 66.2% and 55.0%, respectively, and the median progression-free survival (PFS) period was 43.5 months. The univariate analysis showed that satellite nodules, liver capsule invasion and postoperative recurrence time were associated with OS (P<0.05), and long-term heavy drinking, satellite nodules and liver capsule invasion with PFS (P<0.05). The multivariate analysis indicated that long-term heavy drinking was an independent factor influencing the progression-free survival period of patients with small hepatocellular carcinoma (P=0.003) and postoperative recurrence time and liver capsule invasion were independent factors affecting their overall survival period (P<0.05). Conclusions: The treatment of small hepatocellular carcinoma still concentrates on the active treatment of surgery. It is beneficial to patients to minimize the resection scope of normal liver under the premise of R0 removal of tumor. Postoperative recurrence time of ≤2 years suggests poor prognosis of small hepatocellular carcinoma. Long-term heavy drinking can accelerate the recurrence of small hepatocellular carcinoma.

Explicit consideration of spatial hydrogen bonding direction for activity coefficient prediction based on implicit solvation calculations

Directional hydrogen bonding is introduced to implicit solvation calculations for improved prediction of solvation properties and phase equilibria of associating fluids.

Molecular dynamics study on the nucleation of methane + tetrahydrofuran mixed guest hydrate

The nucleation of methane (CH4), tetrahydrofuran (THF), and CH4 + THF hydrates are investigated by microsecond MD simulations. These three systems exhibit distinct structural developments in the aqueous phase quantified by the formation of cage structures of hydrogen bonded water molecules. The development of a cluster of cages in the CH4 system is limited by the scarce CH4 molecules in the solution, while in the THF system it is limited by the short lifetime of cages. In the CH4 + THF mixed guest system, a small cluster of caged CH4 molecules can be rapidly stabilized by abundant neighboring cages of THF molecules. Therefore, the induction time of the CH4 + THF mixed guest system is found to be significantly shorter than that of the pure CH4 and pure THF systems. Furthermore, the structure of cages found in the initially formed cage clusters are often different from the typical 5(12)6(n) (n = 0, 2, 3, 4) cages observed in clathrate hydrate systems. The cluster of cages may grow or transform into structure I or II clathrate hydrate in the later stages.

Assessment of density functional methods for exciton binding energies and related optoelectronic properties

Mean absolute error (MAE) in exciton binding energy (E_b) from 9 DFT methods against benchmark CCSD and EOM-CCSD.

Theoretical study on the torsional potential of alkyl, donor, and acceptor substituted bithiophene: the hidden role of noncovalent interaction and backbone conjugation

The torsional potential of substituted bithiophene is influenced not only by steric repulsion, but also by backbone conjugation and noncovalent interactions.