Physics

Size-selected anionic silicon clusters, Sin- (n=14-20), have been investigated by photoelectron spectroscopy and density functional theory (DFT) calculations. Low-energy structures of the clusters are globally searched for by using a genetic algorithm based on DFT calculations. The electronic density of states and VDEs have been simulated by using ten DFT functionals and compared to the experimental results. We systematically evaluated the DFT functionals for the calculation of the energetics of silicon clusters. CCSD(T) single-point energies based on MP2 optimized geometries for selected isomers of Sin- are also used as benchmark for the energy sequence. The HSE06 functional with aug-cc-pVDZ basis set is found to show the best performance. Our global minimum search corroborates that most of the lowest-energy structures of Sin- (n=14-20) clusters can be derived from assembling tricapped trigonal prisms (TTP) in various ways. For most sizes previous structures are confirmed, whereas for Si20- a new structure has been found.

In molecular dynamics (MD) simulation, force field determines the capability of an individual model in capturing physical and chemistry properties. The method for generating proper parameters of the force field form is the key component for computational research in chemistry, biochemistry, and condensed-phase physics. Our study showed that the feasibility to predict experimental condensed phase properties (i.e., density and heat of vaporization) of methanol through problem specific force field from only quantum chemistry information. To acquire the satisfying parameter sets of the force field, the genetic algorithm (GA) is the main optimization method. For electrostatic potential energy, we optimized both the electrostatic parameters of methanol using the GA method, which leads to low deviations of between the quantum mechanics (QM) calculations and the GA optimized parameters. We optimized the van der Waals (vdW) parameters both using GA and guided GA methods by calibrating interaction energy of various methanol homo-clusters, such as nonamers, undecamers, or tridecamers. Excellent agreement between the training dataset from QM calculations (i.e., MP2) and GA optimized parameters can be achieved. However, only the guided GA method, which eliminates the overestimation of interaction energy from MP2 calculations in the optimization process, provides proper vdW parameters for MD simulation to get the condensed phase properties (i.e., density and heat of vaporization) of methanol. Throughout the whole optimization process, the experimental value were not involved in the objective functions, but were only used for the purpose of justifying models (i.e., nonamers, undecamers, or tridecamers) and validating methods (i.e., GA or guided GA). Our method shows the possibility of developing descriptive polarizable force field using only QM calculations.

Attosecond pulses in the soft-X-ray (SXR) to water-window energy region offer the tools for creating and studying target specific localized inner-shell electrons or holes in materials, enable monitoring or controlling charge and energy flows in a dynamic system on attosecond timescales. Recently, a number of laboratories have reported generation of continuum harmonics in the hundred-electron-volt to kilovolt region with few-cycle long-wavelength mid-infrared lasers. These harmonics have the bandwidth to support pulses with duration of few- to few-ten attoseconds. But harmonics generated in a gas medium have attochirps that cannot be fully compensated by materials over a broad spectral range; thus, realistically what are the typical shortest attosecond pulses that one can generate? To answer this question, it is essential that the temporal attosecond pulses be accurately characterized. By re-analyzing the soft X-ray harmonics reported in three recent experiments \cite{chang_natcom2017,Thomas_OE2017,Bieger_2017PRX} using a newly developed broadband phase retrieval algorithm, we show that their generated attosecond pulses are all longer than about 60 as. Since broadband pulses tend to have high-order chirps away from the spectral center of the pulse, the algorithm has to be able to retrieve accurately the phase over the whole bandwidth. Our re-evaluated pulse durations are found to be longer than those previously reported. We also introduce the autocorrelation (AC) of the streaking spectrogram. By comparing the ACs from the experiments and from the retrieved SXR pulses, the accuracy of the retrieved results can be directly visualized to ensure that correct phases have been obtained. Our retrieval method is fast and accurate, and it shall provide a powerful tool for the metrology of few-ten-attosecond pulses in the future.

The weakly-bound molecular clusters (D2O)2-CO and (D2O)3-CO are observed in the C-O stretch fundamental region (~2150 cm-1), and their rotationally-resolved infrared spectra yield precise rotational parameters. The corresponding H2O clusters are also observed, but their bands are broadened by predissociation, preventing detailed analysis. The rotational parameters are insufficient in themselves to determine cluster structures, so ab initio calculations are employed, and good agreement between the experimental and theory is found for the most stable cluster isomers, yielding the basic cluster geometries as well as confirming the assignments to (D2O)2-CO and (D2O)3-CO. The trimer, (D2O)2-CO, has a near-planar geometry with one D atom from each D2O slightly out of the plane. The tetramer, (D2O)3-CO, has the water molecules arranged in a cyclic quasi-planar ring similar to the water trimer, with the carbon monoxide located 'above' the ring and roughly parallel to its plane. The tunneling effects observed in the free water dimer and trimer are quenched by the presence of CO. The previously observed water-CO dimer together with the trimer and tetramer reported here represent the first three steps in the solvation of carbon monoxide.

Solving the atomic structure of metallic clusters is fundamental to understanding their optical, electronic, and chemical properties. We report the structure of Au 146 (p-MBA) 57 at subatomic resolution (0.85 Å) using electron diffraction (MicroED) and atomic resolution by X-ray diffraction. The 146 gold atoms may be decomposed into two constituent sets consisting of 119 core and 27 peripheral atoms. The core atoms are organized in a twinned FCC structure whereas the surface gold atoms follow a C 2 rotational symmetry about an axis bisecting the twinning plane. The protective layer of 57 p-MBAs fully encloses the cluster and comprises bridging, monomeric, and dimeric staple motifs. Au 146 (p-MBA) 57 is the largest cluster observed exhibiting a bulk-like FCC structure as well as the smallest gold particle exhibiting a stacking fault.

We demonstrate an efficient control of 29 Si nuclear spin orientation for specific lattice sites near 31 P donors in silicon crystals at temperatures below 1 K and in high magnetic field of 4.6 T. Excitation of the forbidden electron-nuclear transitions leads to a pattern of narrow holes and peaks in the ESR lines of 31 P. The pattern originates from dynamic polarization the 29 Si nuclear spins near the donors via the solid effect. This method can be used for initialization of qubits based on 29 Si nuclear spins in the all-silicon quantum computer. In comparison, polarization of 29 Si performed by pumping the allowed ESR transitions, did not create any patterns. Instead, a single narrow spectral hole was burnt in the ESR line. The difference is explained by a rapid spin diffusion during the microwave pumping of the allowed transitions.

Since its discovery in 1938, superfluidity has been the subject of much investigation because it provides a unique example of a macroscopic manifestation of quantum mechanics. About 60 years later, scientists successfully observed this phenomenon in the microscopic world though the spectroscopic Andronikashvili experiment in helium nano-droplets. This reduction of scale suggests that not only helium but also para-H2 (pH2) can be a candidate for superfluidity. This expectation is based on the fact that the smaller number of neighbours and surface effects of a finite-size cluster may hinder solidification and promote a liquid-like phase. The first prediction of superfluidity in pH2 clusters was reported in 1991 based on quantum Monte Carlo simulations. The possible superfluidity of pH2 was later indirectly observed in a spectroscopic Andronikashvili experiment in 2000. Since then, a growing number of studies have appeared, and theoretical simulations have been playing a special role because they help guide and interpret experiments. In this review, we go over the theoretical studies of pH2 superfluid clusters since the experiment of 2000. We provide a historical perspective and introduce the basic theoretical formalism along with key experimental advances. We then present illustrative results of the theoretical studies and comment on the possible future developments in the field. We include sufficient theoretical details such that the review can serve as a guide for newcomers to the field.

We find that negative charges on an armchair single-walled carbon nanotube (SWCNT) can significantly enhance the migration of a carbon adatom on the external surfaces of SWCNTs, along the direction of the tube axis. Nanotube charging results in stronger binding of adatoms to SWCNTs and consequent longer lifetimes of adatoms before desorption, which in turn increases their migration distance several orders of magnitude. These results support the hypothesis of diffusion enhanced SWCNT growth in the volume of arc plasma. This process could enhance effective carbon flux to the metal catalyst.

The addition of small impurities, such as a single proton charge carrier, in noble gas clusters has recently been shown to have considerable effects on their geometries and stabilities. Here we report on a mass spectrometric study of cationic clusters of N 2 molecules and the effects that adding hydrogen, in the form of D 2 , has on the systems. Protonated nitrogen clusters formed by the breakup of D 2 are shown to have similar behaviors as protonated rare gas clusters. For larger systems consisting of different mixtures of intact N 2 and D 2 , the different molecular species are found to sometimes be interchangeable with regards to magic numbers. This is especially true for the (N 2 ) n (D 2 ) m D + systems with n+m=17 , which is particularly abundant for all measured combinations of n and m .

Gold nanoparticles (AuNPs) have been shown to possess properties beneficial for the treatment of cancerous tumors by acting as radiosensitizers for both photon and ion radiation. Blood circulation time is usually increased by coating the AuNPs with poly(ethylene glycol) (PEG) ligands. The effectiveness of the PEG coating, however, depends on both the ligand surface density and length of the PEG molecules, making it important to understand the structure of the coating. In this paper the thickness, ligand surface density, and density of the PEG coating is studied with classical molecular dynamics using the software package MBN Explorer. AuNPs consisting of 135 atoms (approximately 1.4 nm diameter) in a water medium have been studied with the number of PEG ligands varying between 32 and 60. We find that the thickness of the coating is only weakly dependent on the surface ligand density and that the degree of water penetration increased when there is a smaller number of attached ligands.