A QM/MM equation-of-motion coupled-cluster approach for predicting semiconductor color-center structure and emission frequencies

Valence excitation spectra are computed for all deep-center silicon-vacancy defect types in 3C, 4H, and 6H silicon carbide (SiC) and comparisons are made with literature photoluminescence measurements. Nuclear geometries surrounding the defect centers are optimized within a Gaussian basis-set framework using many-body perturbation theory or density functional theory (DFT) methods, with computational expenses minimized by a QM/MM technique called SIMOMM. Vertical excitation energies are subsequently obtained by applying excitation-energy, electron-attached, and ionized equation-of-motion coupled-cluster (EOMCC) methods, where appropriate, as well as time-dependent (TD) DFT, to small models including only a few atoms adjacent to the defect center. We consider the relative quality of various EOMCC and TD-DFT methods for (i) energy-ordering potential ground states differing incrementally in charge and multiplicity, (ii) accurately reproducing experimentally measured photoluminescence peaks, and (iii) energy-ordering defects of different types occurring within a given polytype. The extensibility of this approach to transition-metal defects is also tested by applying it to silicon-substitutional chromium defects in SiC and comparing with measurements. It is demonstrated that, when used in conjunction with SIMOMM-optimized geometries, EOMCC-based methods can provide a reliable prediction of the ground-state charge and multiplicity, while also giving a quantitative description of the photoluminescence spectra, accurate to within 0.1 eV of measurement in all cases considered.Comment: 13 pages, 4 figures, 6 tables, 5 equations, 100 reference

The closo-Si\u3csub\u3e12\u3c/sub\u3eC\u3csub\u3e12\u3c/sub\u3e Molecule from Cluster to Crystal: A Theoretical Prediction

The structure of closo-Si12C12 is unique among stable SinCm isomers (n, m \u3e 4) because of its high symmetry, π–π stacking of C6 rings and unsaturated silicon atoms at symmetrical peripheral positions. Dimerization potential surfaces reveal various dimerization reactions that form between two closo-Si12C12 molecules through Si–Si bonds at unsaturated Si atoms. As a result the closo-Si12C12 molecule is capable of polymerization to form stable 1D polymer chains, 2D crystal layers, and 3D crystals. 2D crystal structures formed by side-side polymerization satisfy eight Si valences on each monomer without large distortion of the monomer structure. 3D crystals are formed by stacking 2D structures in the Z direction, preserving registry of C6 rings in monomer moiety

Theoretical Investigation of Stabilities and Optical Properties of Si\u3csub\u3e12\u3c/sub\u3eC\u3csub\u3e12\u3c/sub\u3e Clusters

By sorting through hundreds of globally stable Si12C12 isomers using a potential surface search and using simulated annealing, we have identified low-energy structures. Unlike isomers knit together by Si–C bonds, the lowest energy isomers have segregated carbon and silicon regions that maximize stronger C–C bonding. Positing that charge separation between the carbon and silicon regions would produce interesting optical absorption in these cluster molecules, we used time-dependent density functional theory to compare the calculated optical properties of four isomers representing structural classes having different types of silicon and carbon segregation regions. Absorptions involving charge transfer between segregated carbon and silicon regions produce lower excitation energies than do structures having alternating Si–C bonding for which frontier orbital charge transfer is exclusively from separated carbon atoms to silicon atoms. The most stable Si12C12 isomer at temperatures below 1100 K is unique as regards its high symmetry and large optical oscillator strength in the visible blue. Its high-energy and low-energy visible transitions (1.15 eV and 2.56 eV) are nearly pure one-electron silicon-to-carbon transitions, while an intermediate energy transition (1.28 eV) is a nearly pure carbon-to-silicon one-electron charge transfer

Searching for Stable Si\u3csub\u3en\u3c/sub\u3eC\u3csub\u3en\u3c/sub\u3e Clusters: Combination of Stochastic Potential Surface Search and Pseudopotential Plane-Wave Car-Parinello Simulated Annealing Simulations

To find low energy SinCn structures out of hundreds to thousands of isomers we have developed a general method to search for stable isomeric structures that combines Stochastic Potential Surface Search and Pseudopotential Plane-Wave Density Functional Theory Car-Parinello Molecular Dynamics simulated annealing (PSPW-CPMD-SA). We enhanced the Sunders stochastic search method to generate random cluster structures used as seed structures for PSPW-CPMD-SA simulations. This method ensures that each SA simulation samples a different potential surface region to find the regional minimum structure. By iterations of this automated, parallel process on a high performance computer we located hundreds to more than a thousand stable isomers for each SinCn cluster. Among these, five to 10 of the lowest energy isomers were further optimized using B3LYP/cc-pVTZ method. We applied this method to SinCn (n = 4–12) clusters and found the lowest energy structures, most not previously reported. By analyzing the bonding patterns of low energy structures of each SinCn cluster, we observed that carbon segregations tend to form condensed conjugated rings while Si connects to unsaturated bonds at the periphery of the carbon segregation as single atoms or clusters when n is small and when n is large a silicon network spans over the carbon segregation region

Semiconductor Color-center Structure and Excitation Spectra: Equation-of-motion Coupled-cluster Description of Vacancy and Transition-metal Defect Photoluminescence

Valence excitation spectra are computed for deep-center silicon-vacancy defects in 3C, 4H, and 6H silicon carbide (SiC), and comparisons are made with literature photoluminescence measurements. Optimizations of nuclear geometries surrounding the defect centers are performed within a Gaussian basis-set framework using many-body perturbation theory or density functional theory (DFT) methods, with computational expenses minimized by a QM/MM technique called SIMOMM. Vertical excitation energies are subsequently obtained by applying excitation-energy, electron-attached, and ionized equation-of-motion coupled-cluster (EOMCC) methods, where appropriate, as well as time-dependent (TD) DFT, to small models including only a few atoms adjacent to the defect center. We consider the relative quality of various EOMCC and TD-DFT methods for (i) energy-ordering potential ground states differing incrementally in charge and multiplicity, (ii) accurately reproducing experimentally measured photoluminescence peaks, and (iii) energy-ordering defects of different types occurring within a given polytype. The extensibility of this approach to transition-metal defects is also tested by applying it to silicon-substituted chromium defects in SiC and comparing with measurements. It is demonstrated that, when used in conjunction with SIMOMM-optimized geometries, EOMCC-based methods can provide a reliable prediction of the ground-state charge and multiplicity, while also giving a quantitative description of the photoluminescence spectra, accurate to within 0.1 eV of measurement for all cases considered. Abstract ©2018 American Physical Societ

Predictive coupled-cluster isomer orderings for some Sin{}_nCm{}_m (m,n≤12m, n\le 12) clusters; A pragmatic comparison between DFT and complete basis limit coupled-cluster benchmarks

The accurate determination of the preferred ${\rm Si}_{12}{\rm C}_{12}$ isomer is important to guide experimental efforts directed towards synthesizing SiC nano-wires and related polymer structures which are anticipated to be highly efficient exciton materials for opto-electronic devices. In order to definitively identify preferred isomeric structures for silicon carbon nano-clusters, highly accurate geometries, energies and harmonic zero point energies have been computed using coupled-cluster theory with systematic extrapolation to the complete basis limit for set of silicon carbon clusters ranging in size from SiC$_3$ to ${\rm Si}_{12}{\rm C}_{12}$. It is found that post-MBPT(2) correlation energy plays a significant role in obtaining converged relative isomer energies, suggesting that predictions using low rung density functional methods will not have adequate accuracy. Utilizing the best composite coupled-cluster energy that is still computationally feasible, entailing a 3-4 SCF and CCSD extrapolation with triple-$\zeta$ (T) correlation, the {\it closo} ${\rm Si}_{12}{\rm C}_{12}$ isomer is identified to be the preferred isomer in support of previous calculations [J. Chem. Phys. 2015, 142, 034303]. Additionally we have investigated more pragmatic approaches to obtaining accurate silicon carbide isomer energies, including the use of frozen natural orbital coupled-cluster theory and several rungs of standard and double-hybrid density functional theory. Frozen natural orbitals as a way to compute post MBPT(2) correlation energy is found to be an excellent balance between efficiency and accuracy

Valence and Charge-transfer Optical Properties for Some Si\u3csub\u3en\u3c/sub\u3eC\u3csub\u3em\u3c/sub\u3e (m, n ≤ 12) Clusters: Comparing TD-DFT, Complete-basis-limit EOMCC, and Benchmarks from Spectroscopy

Accurate optical characterization of the closo-Si12C12 molecule is important to guide experimental efforts toward the synthesis of nano-wires, cyclic nano-arrays, and related array structures, which are anticipated to be robust and efficient exciton materials for opto-electronic devices. Working toward calibrated methods for the description of closo-Si12C12 oligomers, various electronic structure approaches are evaluated for their ability to reproduce measured optical transitions of the SiC2, Si2Cn (n = 1–3), and Si3Cn (n = 1, 2) clusters reported earlier by Steglich and Maier [Astrophys. J. 801, 119 (2015)]. Complete-basis-limit equation-of-motion coupled-cluster (EOMCC) results are presented and a comparison is made between perturbative and renormalized non-iterative triples corrections. The effect of adding a renormalized correction for quadruples is also tested. Benchmark test sets derived from both measurement and high-level EOMCC calculations are then used to evaluate the performance of a variety of density functionals within the time-dependent density functional theory (TD-DFT) framework. The best-performing functionals are subsequently applied to predict valence TD-DFT excitation energies for the lowest-energy isomers of SinC and Sin−1C7−n (n = 4–6). TD-DFT approaches are then applied to the SinCn (n = 4–12) clusters and unique spectroscopic signatures of closo-Si12C12 are discussed. Finally, various long-range corrected density functionals, including those from the CAM-QTP family, are applied to a charge-transfer excitation in a cyclic (Si4C4)4 oligomer. Approaches for gauging the extent of charge-transfer character are also tested and EOMCC results are used to benchmark functionals and make recommendations

Different Effects of Farrerol on an OVA-Induced Allergic Asthma and LPS-induced Acute Lung Injury

BACKGROUND: Farrerol, isolated from rhododendron, has been shown to have the anti-bacterial activity, but no details on the anti-inflammatory activity. We further evaluated the effects of this compound in two experimental models of lung diseases. METHODOLOGY/PRINCIPAL FINDINGS: For the asthma model, female BALB/c mice were challenged with ovalbumin (OVA), and then treated daily with farrerol (20 and 40 mg/kg, i.p.) as a therapeutic treatment from day 22 to day 26 post immunization. To induce acute lung injury, female BALB/c mice were injected intranasally with LPS and treated with farrerol (20 and 40 mg/kg, i.p.) 1 h prior to LPS stimulation. Inflammation in the two different models was determined using ELISA, histology, real-time PCR and western blot. Farrerol significantly regulated the phenotype challenged by OVA, like cell number, Th1 and Th2 cytokines levels in the BALF, the OVA-specific IgE level in the serum, goblet cell hyperplasia in the airway, airway hyperresponsiveness to inhaled methacholine and mRNA expression of chemokines and their receptors. Furthermore, farrerol markedly attenuated the activation of phosphorylation of Akt and nuclear factor-κB (NF-κB) subunit p65 both in vivo and in vitro. However, farrerol has no effect on the acute lung injury model. CONCLUSION/SIGNIFICANCE: Our finding demonstrates that the distinct anti-inflammatory effect of farrerol in the treatment of asthma acts by inhibiting the PI3K and NF-κB pathway

Association of PPARγ2 polymorphisms with carcass and meat quality traits in a Pietrain x Jinhua F2 population

The PPARγ2 gene is a key regulator of both proliferation and preadipocyte differentiation in mammals. Herein its genotype and allele frequencies were analyzed using PCR-SSCP in eight pig breeds (N = 416). Two kinds of polymorphisms of the PPARγ2 gene were detected, including a previously reported shift SNP A177G (Met59Val) in exon 1 and a novel silent mutation G876A in exon 5. The results revealed that European pig breeds carry a higher allele A frequency at the A177G locus and a fixed GG genotype at the G876A locus. Allele A at the G876A locus was only found in Jinhua pigs. The association between haplotype (A177G/G876A) and carcass and meat quality traits was analyzed in a Pietrain x Jinhua F2 population (N = 248). The PPARγ2 gene was found to be significantly associated with backfat thickness at the shoulder (p < 0.05), 6–7th ribs (p < 0.01), last rib (p < 0.01), gluteus medius (p <0.05) and ham weight (p < 0.01). Significant effects of different haplotypes on ham weight and backfat thickness at the 6–7th ribs, last rib, and gluteus medius were also observed

A Human-Specific De Novo Protein-Coding Gene Associated with Human Brain Functions

To understand whether any human-specific new genes may be associated with human brain functions, we computationally screened the genetic vulnerable factors identified through Genome-Wide Association Studies and linkage analyses of nicotine addiction and found one human-specific de novo protein-coding gene, FLJ33706 (alternative gene symbol C20orf203). Cross-species analysis revealed interesting evolutionary paths of how this gene had originated from noncoding DNA sequences: insertion of repeat elements especially Alu contributed to the formation of the first coding exon and six standard splice junctions on the branch leading to humans and chimpanzees, and two subsequent substitutions in the human lineage escaped two stop codons and created an open reading frame of 194 amino acids. We experimentally verified FLJ33706's mRNA and protein expression in the brain. Real-Time PCR in multiple tissues demonstrated that FLJ33706 was most abundantly expressed in brain. Human polymorphism data suggested that FLJ33706 encodes a protein under purifying selection. A specifically designed antibody detected its protein expression across human cortex, cerebellum and midbrain. Immunohistochemistry study in normal human brain cortex revealed the localization of FLJ33706 protein in neurons. Elevated expressions of FLJ33706 were detected in Alzheimer's brain samples, suggesting the role of this novel gene in human-specific pathogenesis of Alzheimer's disease. FLJ33706 provided the strongest evidence so far that human-specific de novo genes can have protein-coding potential and differential protein expression, and be involved in human brain functions