Numerical simulation of interdigitated back contact hetero-junction solar cells

The goal of this thesis is the application of an opto-electronic numerical simulation to heterojunction silicon solar cells featuring an all back contact architecture (Interdigitated Back Contact Hetero-Junction IBC-HJ). The studied structure exhibits both metal contacts, emitter and base, at the back surface of the cell with the objective to reduce the optical losses due to the shadowing by front contact of conventional photovoltaic devices. Overall, IBC-HJ are promising low-cost alternatives to monocrystalline wafer-based solar cells featuring front and back contact schemes, in fact, for IBC-HJ the high concentration doping diffusions are replaced by low-temperature deposition processes of thin amorphous silicon layers. Furthermore, another advantage of IBC solar cells with reference to conventional architectures is the possibility to enable a low-cost assembling of photovoltaic modules, being all contacts on the same side. A preliminary extensive literature survey has been helpful to highlight the specific critical aspects of IBC-HJ solar cells as well as the state-of-the-art of their modeling, processing and performance of practical devices. In order to perform the analysis of IBC-HJ devices, a two-dimensional (2-D) numerical simulation flow has been set up. A commercial device simulator based on finite-difference method to solve numerically the whole set of equations governing the electrical transport in semiconductor materials (Sentuarus Device by Synopsys) has been adopted. The first activity carried out during this work has been the definition of a 2-D geometry corresponding to the simulation domain and the specification of the electrical and optical properties of materials. In order to calculate the main figures of merit of the investigated solar cells, the spatially resolved photon absorption rate map has been calculated by means of an optical simulator. Optical simulations have been performed by using two different methods depending upon the geometrical features of the front interface of the solar cell: the transfer matrix method (TMM) and the raytracing (RT). The first method allows to model light prop-agation by plane waves within one-dimensional spatial domains under the assumption of devices exhibiting stacks of parallel layers with planar interfaces. In addition, TMM is suitable for the simulation of thin multi-layer anti reflection coating layers for the reduction of the amount of reflected light at the front interface. Raytracing is required for three-dimensional optical simulations of upright pyramidal textured surfaces which are widely adopted to significantly reduce the reflection at the front surface. The optical generation profiles are interpolated onto the electrical grid adopted by the device simulator which solves the carriers transport equations coupled with Poisson and continuity equations in a self-consistent way. The main figures of merit are calculated by means of a postprocessing of the output data from device simulation. After the validation of the simulation methodology by means of comparison of the simulation result with literature data, the ultimate efficiency of the IBC-HJ architecture has been calculated. By accounting for all optical losses, IBC-HJ solar cells result in a theoretical maximum efficiency above 23.5% (without texturing at front interface) higher than that of both standard homojunction crystalline silicon (Homogeneous Emitter HE) and front contact heterojuction (Heterojunction with Intrinsic Thin layer HIT) solar cells. However it is clear that the criticalities of this structure are mainly due to the defects density and to the poor carriers transport mobility in the amorphous silicon layers. Lastly, the influence of the most critical geometrical and physical parameters on the main figures of merit have been investigated by applying the numerical simulation tool set-up during the first part of the present thesis. Simulations have highlighted that carrier mobility and defects level in amorphous silicon may lead to a potentially significant reduction of the conversion efficiency.

Shear behaviour of reinforced cconcrete slab under concentrated load: an investigation through non-linear and sequentially linear analysis

English: The assessment of safety in existing bridges and viaducts led the Ministry of Public Works of the Netherlands to finance a specific campaing aimed at the study of the response of the elements of these infrastructures. Therefore, this activity is focused on the investigation of the behaviour of reinforced concrete slabs under concentrated loads, adopting finite element modeling and comparison with experimental results. These elements are characterized by shear behaviour and crisi, whose modeling is, from a computational point of view, a hard challeng, due to the brittle behavior combined with three-dimensional effects. The numerical modeling of the failure is studied through Sequentially Linear Analysis (SLA), an alternative Finite Element method, with respect to traditional incremental and iterative approaches. The comparison between the two different numerical techniques represents one of the first works and comparisons in a three-dimensional environment. It's carried out adopting one of the experimental test executed on reinforced concrete slabs as well. The advantage of the SLA is to avoid the well known problems of convergence of typical non-linear analysis, by directly specifying a damage increment, in terms of reduction of stiffness and resistance in particular finite element, instead of load or displacement increasing on the whole structure . For the first time, particular attention has been paid to specific aspects of the slabs, like an accurate constraints modeling and sensitivity of the solution with respect to the mesh density. This detailed analysis with respect to the main parameters proofed a strong influence of the tensile fracture energy, mesh density and chosen model on the solution in terms of force-displacement diagram, distribution of the crack patterns and shear failure mode. The SLA showed a great potential, but it requires a further developments for what regards two aspects of modeling: load conditions (constant and proportional loads) and softening behaviour of brittle materials (like concrete) in the three-dimensional field, in order to widen its horizons in these new contexts of study.

Comparison of three optical tracking systems in a complex navigation scenario

Three-dimensional rotational X-ray imaging with the SIREMOBIL Iso-C3D (Siemens AG, Medical Solutions, Erlangen, Germany) has become a well-established intra-operative imaging modality. In combination with a tracking system, the Iso-C3D provides inherently registered image volumes ready for direct navigation. This is achieved by means of a pre-calibration procedure. The aim of this study was to investigate the influence of the tracking system used on the overall navigation accuracy of direct Iso-C3D navigation. Three models of tracking system were used in the study: Two Optotrak 3020s, a Polaris P4 and a Polaris Spectra system, with both Polaris systems being in the passive operation mode. The evaluation was carried out at two different sites using two Iso-C3D devices. To measure the navigation accuracy, a number of phantom experiments were conducted using an acrylic phantom equipped with titanium spheres. After scanning, a special pointer was used to pinpoint these markers. The difference between the digitized and navigated positions served as the accuracy measure. Up to 20 phantom scans were performed for each tracking system. The average accuracy measured was 0.86 mm and 0.96 mm for the two Optotrak 3020 systems, 1.15 mm for the Polaris P4, and 1.04 mm for the Polaris Spectra system. For the Polaris systems a higher maximal error was found, but all three systems yielded similar minimal errors. On average, all tracking systems used in this study could deliver similar navigation accuracy. The passive Polaris system showed ? as expected ? higher maximal errors; however, depending on the application constraints, this might be negligible.

Clinical and biochemical description of a novel CYP21A2 gene mutation 962_963insA using a new 3D model for the P450c21 protein

OBJECTIVE: A severely virilized 46, XX newborn girl was referred to our center for evaluation and treatment of congenital adrenal hyperplasia (CAH) because of highly elevated 17alpha-hydroxyprogesterone levels at newborn screening; biochemical tests confirmed the diagnosis of salt-wasting CAH. Genetic analysis revealed that the girl was compound heterozygote for a previously reported Q318X mutation in exon 8 and a novel insertion of an adenine between nucleotides 962 and 963 in exon 4 of the CYP21A2 gene. This 962_963insA mutation created a frameshift leading to a stop codon at amino acid 161 of the P450c21 protein. AIM AND METHODS: To better understand structure-function relationships of mutant P450c21 proteins, we performed multiple sequence alignments of P450c21 with three mammalian P450s (P450 2C8, 2C9 and 2B4) with known structures as well as with human P450c17. Comparative molecular modeling of human P450c21 was then performed by MODELLER using the X-ray crystal structure of rabbit P450 2B4 as a template. RESULTS: The new three dimensional model of human P450c21 and the sequence alignment were found to be helpful in predicting the role of various amino acids in P450c21, especially those involved in heme binding and interaction with P450 oxidoreductase, the obligate electron donor. CONCLUSION: Our model will help in analyzing the genotype-phenotype relationship of P450c21 mutations which have not been tested for their functional activity in an in vitro assay.

Different amino acid substitutions at the same position in rhodopsin lead to distinct phenotypes

PURPOSE: Identification of a novel rhodopsin mutation in a family with retinitis pigmentosa and comparison of the clinical phenotype to a known mutation at the same amino acid position. METHODS: Screening for mutations in rhodopsin was performed in 78 patients with retinitis pigmentosa. All exons and flanking intronic regions were amplified by PCR, sequenced, and compared to the reference sequence derived from the National Center for Biotechnology Information (NCBI, Bethesda, MD) database. Patients were characterized clinically according to the results of best corrected visual acuity testing (BCVA), slit lamp examination (SLE), funduscopy, Goldmann perimetry (GP), dark adaptometry (DA), and electroretinography (ERG). Structural analyses of the rhodopsin protein were performed with the Swiss-Pdb Viewer program available on-line (http://www.expasy.org.spdvbv/ provided in the public domain by Swiss Institute of Bioinformatics, Geneva, Switzerland). RESULTS: A novel rhodopsin mutation (Gly90Val) was identified in a Swiss family of three generations. The pedigree indicated autosomal dominant inheritance. No additional mutation was found in this family in other autosomal dominant genes. The BCVA of affected family members ranged from 20/25 to 20/20. Fundus examination showed fine pigment mottling in patients of the third generation and well-defined bone spicules in patients of the second generation. GP showed concentric constriction. DA demonstrated monophasic cone adaptation only. ERG revealed severely reduced rod and cone signals. The clinical picture is compatible with retinitis pigmentosa. A previously reported amino acid substitution at the same position in rhodopsin leads to a phenotype resembling night blindness in mutation carriers, whereas patients reported in the current study showed the classic retinitis pigmentosa phenotype. The effect of different amino acid substitutions on the three-dimensional structure of rhodopsin was analyzed by homology modeling. Distinct distortions of position 90 (shifts in amino acids 112 and 113) and additional hydrogen bonds were found. CONCLUSIONS: Different amino acid substitutions at position 90 of rhodopsin can lead to night blindness or retinitis pigmentosa. The data suggest that the property of the substituted amino acid distinguishes between the phenotypes.

Do women have impaired regional systolic function in hypertensive heart disease? A 3-dimensional echocardiography study

AIMS: In pressure overload left ventricular (LV) hypertrophy, gender-related differences in global LV systolic function have been previously reported. The goal of this study was to determine regional systolic function of the left ventricle in male and female patients with hypertensive heart disease. METHODS AND RESULTS: Regional LV function was analyzed from multiplane transesophageal echocardiography with three-dimensional (3D) reconstruction of the left ventricle. In 24 patients (13 males and 11 females), four parallel (2 basal and 2 apical) equidistant short axis cross-sections from base to apex were obtained from the reconstructed LV. In each short axis 24 wall-thickness measurements were carried out at 15 degrees intervals at end-diastole and end-systole. Thus, a total of 192 measurements were obtained in each patient. Wall thickening was calculated as difference of end-diastolic and end-systolic wall thickness, and fractional thickening as thickening divided by end-diastolic thickness. Fractional thickening and wall stress were inversely related to end-diastolic wall thickness in both, males and females. Females showed less LV systolic function when compared to males (p<0.001). However, when corrected for wall stress, which was higher in females, there was no gender difference in systolic function. CONCLUSION: There are regional differences in LV systolic function in females and males which are directly related to differences in wall stress. Thus, gender-related differences in LV regional function are load-dependent and not due to structural differences.

Pseudodominant Inheritance of Goitrous Congenital Hypothyroidism Caused by TPO Mutations: Molecular and in Silico Studies

Context and Objective: Most cases of goitrous congenital hypothyroidism (CH) from thyroid dyshormonogenesis 1) follow a recessive mode of inheritance and 2) are due to mutations in the thyroid peroxidase gene (TPO). We report the genetic mechanism underlying the apparently dominant inheritance of goitrous CH in a nonconsanguineous family of French Canadian origin. Design, Setting, and Participants: Two brothers identified by newborn TSH screening had severe hypothyroidism and a goiter with increased (99m)Tc uptake. The mother was euthyroid, but the father and two paternal uncles had also been diagnosed with goitrous CH. After having excluded PAX8 gene mutations, we hypothesized that the underlying defect could be TPO mutations. Results: Both compound heterozygous siblings had inherited a mutant TPO allele carried by their mother (c.1496delC; p.Pro499Argfs2X), and from their father, one brother had inherited a missense mutation (c.1978C-->G; p.Gln660Glu) and the other an insertion (c.1955insT; p.Phe653Valfs15X). The thyroid gland of one uncle who is a compound heterozygote for TPO mutations (p.Phe653Valfs15X/p.Gln660Glu) was removed because of concurrent multiple endocrine neoplasia type 2A. Immunohistochemistry revealed normal TPO staining, implying that Gln660Glu TPO is expressed properly. Modeling of this mutant in silico suggests that its three-dimensional structure is conserved, whereas the electrostatic binding energy between the Gln660Glu TPO and its heme group becomes repulsive. Conclusion: We report a pedigree presenting with pseudodominant goitrous CH due to segregation of three different TPO mutations. Although goitrous CH generally follows a recessive mode of inheritance, the high frequency of TPO mutations carriers may lead to pseudodominant inheritance.

Prediction of spatial orientation and morphology of calcaneonavicular coalitions

BACKGROUND: Calcaneonavicular coalitions (CNC) have been reported to be associated with anatomical aberrations of either the calcaneus and/or navicular bones. These morphological abnormalities may complicate accurate surgical resection. Three-dimensional analysis of spatial orientation and morphological characteristics may help in preoperative planning of resection. MATERIALS AND METHODS: Sixteen feet with a diagnosis of CNC were evaluated by means of 3-D CT modeling. Three angles were defined that were expressed in relation to one reproducible landmark (lateral border of the calcaneus): the dorsoplantar inclination, anteroposterior inclination, and socket angle. The depth and width of the coalitions were measured and calculated to obtain the estimated contact surface. Three-dimensional reconstructions of the calcanei served to evaluate the presence, distortion or absence of the anterior calcaneal facet and presence of a navicular beak. The interrater correlations were assessed in order to obtain values for the accuracy of the measurement methods. Sixteen normal feet were used as controls for comparison of the socket angle; anatomy of the anterior calcaneal facet and navicular beak as well. RESULTS: The dorsoplantar inclination angle averaged 50 degrees (+/-17), the anteroposterior inclination angle 64 degrees (+/-15), and the pathologic socket angle 98 degrees (+/-11). The average contact area was 156 mm(2). Ninety-four percent of all patients in the CNC group revealed a plantar navicular beak. In 50% of those patients the anterior calcaneal facet was replaced by the navicular portion and in 44% the facet was totally missing. In contrast, the socket angle in the control group averaged 77 degrees (+/-18), which was found to be statistically different than the CNC group (p = 0.0004). Only 25% of the patients in the control group had a plantar navicular beak. High, statistically significant interrater correlations were found for all measured angles. CONCLUSION: Computer-aided CT analysis and reconstructions help to determine the spatial orientations of CNC in space and provide useful information in order to anticipate morphological abnormalities of the calcaneus and navicular.

Integration of computational models and experimental characterization to study internal frost damage in cementitious materials

The objective of this doctoral research is to investigate the internal frost damage due to crystallization pore pressure in porous cement-based materials by developing computational and experimental characterization tools. As an essential component of the U.S. infrastructure system, the durability of concrete has significant impact on maintenance costs. In cold climates, freeze-thaw damage is a major issue affecting the durability of concrete. The deleterious effects of the freeze-thaw cycle depend on the microscale characteristics of concrete such as the pore sizes and the pore distribution, as well as the environmental conditions. Recent theories attribute internal frost damage of concrete is caused by crystallization pore pressure in the cold environment. The pore structures have significant impact on freeze-thaw durability of cement/concrete samples. The scanning electron microscope (SEM) and transmission X-ray microscopy (TXM) techniques were applied to characterize freeze-thaw damage within pore structure. In the microscale pore system, the crystallization pressures at sub-cooling temperatures were calculated using interface energy balance with thermodynamic analysis. The multi-phase Extended Finite Element Modeling (XFEM) and bilinear Cohesive Zone Modeling (CZM) were developed to simulate the internal frost damage of heterogeneous cement-based material samples. The fracture simulation with these two techniques were validated by comparing the predicted fracture behavior with the captured damage from compact tension (CT) and single-edge notched beam (SEB) bending tests. The study applied the developed computational tools to simulate the internal frost damage caused by ice crystallization with the two dimensional (2-D) SEM and three dimensional (3-D) reconstructed SEM and TXM digital samples. The pore pressure calculated from thermodynamic analysis was input for model simulation. The 2-D and 3-D bilinear CZM predicted the crack initiation and propagation within cement paste microstructure. The favorably predicted crack paths in concrete/cement samples indicate the developed bilinear CZM techniques have the ability to capture crack nucleation and propagation in cement-based material samples with multiphase and associated interface. By comparing the computational prediction with the actual damaged samples, it also indicates that the ice crystallization pressure is the main mechanism for the internal frost damage in cementitious materials.

Structural model of the CopA copper ATPase of Enterococcus hirae based on chemical cross-linking

The CopA copper ATPase of Enterococcus hirae belongs to the family of heavy metal pumping CPx-type ATPases and shares 43% sequence similarity with the human Menkes and Wilson copper ATPases. Due to a lack of suitable protein crystals, only partial three-dimensional structures have so far been obtained for this family of ion pumps. We present a structural model of CopA derived by combining topological information obtained by intramolecular cross-linking with molecular modeling. Purified CopA was cross-linked with different bivalent reagents, followed by tryptic digestion and identification of cross-linked peptides by mass spectrometry. The structural proximity of tryptic fragments provided information about the structural arrangement of the hydrophilic protein domains, which was integrated into a three-dimensional model of CopA. Comparative modeling of CopA was guided by the sequence similarity to the calcium ATPase of the sarcoplasmic reticulum, Serca1, for which detailed structures are available. In addition, known partial structures of CPx-ATPase homologous to CopA were used as modeling templates. A docking approach was used to predict the orientation of the heavy metal binding domain of CopA relative to the core structure, which was verified by distance constraints derived from cross-links. The overall structural model of CopA resembles the Serca1 structure, but reveals distinctive features of CPx-type ATPases. A prominent feature is the positioning of the heavy metal binding domain. It features an orientation of the Cu binding ligands which is appropriate for the interaction with Cu-loaded metallochaperones in solution. Moreover, a novel model of the architecture of the intramembranous Cu binding sites could be derived.

Three members of the connective tissue growth factor family CCN are differentially regulated by mechanical stress

Expression of connective tissue growth factor (CTGF), a member of the CCN gene family, is known to be significantly induced by mechanical stress. We have therefore investigated whether other members of the CCN gene family, including Cyr61 and Nov, might reveal a similar stress-dependent regulation. Fibroblasts growing under stressed conditions within a three-dimensional collagen gel showed at least a 15 times higher level of Cyr61 mRNA than cells growing under relaxed conditions. Upon relaxation, the decline of the Cyr61 mRNA to a lower level occurred within 2 h, and was thus quicker than the response of CTGF. The regulation was fully reversible when stress was reapplied. Thus, Cyr61 represents another typical example of a stress-responsive gene. The level of the Nov mRNA was low in the stressed state, but increased in the relaxed state. This CCN gene therefore shows an inverted regulation relative to that of Cyr61 and CTGF. Inhibition of protein kinases by means of staurosporine suppressed the stress-induced expression of Cyr61 and CTGF. Elevated levels of cAMP induced by forskolin mimicked the effects of relaxation on the regulation of Cyr61, CTGF and Nov. Thus, adenylate cyclase as well as one or several protein kinases might be involved in the mechanoregulation of these CCN genes.

Molecular structure of antibodies against small ligands

Antibodies which bind bioactive ligands can serve as a template for the generation of a second antibody which may react with the physiological receptor. This phenomenon of molecular mimicry by antibodies has been described in a variety of systems. In order to understand the chemical and molecular mechanisms involved in these interactions, monoclonal antibodies directed against two pharmacologically active alkaloids, morphine and nicotine, were carefully studied using experimental and theoretical molecular modeling techniques. The molecular characterization of these antibodies involved binding studies with ligand analogs and determination of the variable region amino acid sequence. A three-dimensional model of the anti-morphine binding site was constructed using computational and graphics display techniques. The antibody response in BALB/c mice to morphine appears relatively restricted, in that all of the antibodies examined in this study contained a $\lambda$ light chain, which is normally found in only 5% of mouse immunoglobulins. This study represents the first use of theoretical and experimental modeling techniques to describe the antigen binding site of a mouse Fv region containing a $\lambda$ light chain. The binding site model indicates that a charged glutamic acid residue and aromatic side chains are key features in ionic and hydrophobic interactions with the ligand morphine. A glutamic acid residue is found in the identical position in the anti-nicotine antibody and may play a role in binding nicotine. ^

Simulating Cortical Development as a Self Constructing Process: A Novel Multi-Scale Approach Combining Molecular and Physical Aspects

Current models of embryological development focus on intracellular processes such as gene expression and protein networks, rather than on the complex relationship between subcellular processes and the collective cellular organization these processes support. We have explored this collective behavior in the context of neocortical development, by modeling the expansion of a small number of progenitor cells into a laminated cortex with layer and cell type specific projections. The developmental process is steered by a formal language analogous to genomic instructions, and takes place in a physically realistic three-dimensional environment. A common genome inserted into individual cells control their individual behaviors, and thereby gives rise to collective developmental sequences in a biologically plausible manner. The simulation begins with a single progenitor cell containing the artificial genome. This progenitor then gives rise through a lineage of offspring to distinct populations of neuronal precursors that migrate to form the cortical laminae. The precursors differentiate by extending dendrites and axons, which reproduce the experimentally determined branching patterns of a number of different neuronal cell types observed in the cat visual cortex. This result is the first comprehensive demonstration of the principles of self-construction whereby the cortical architecture develops. In addition, our model makes several testable predictions concerning cell migration and branching mechanisms.

ON THE EXISTENCE OF SHOCKS IN IRRADIATED EXOPLANETARY ATMOSPHERES

Supersonic flows are expected to exist in the atmospheres of irradiated exoplanets, but the question of whether shocks develop lingers. Specifically, it reduces to whether continuous flow in a closed loop may become supersonic and if some portions of the supersonic flow steepen into shocks. We first demonstrate that continuous, supersonic flow may exist in two flavors: isentropic and non-isentropic, with shocks being included in the latter class of solutions. Supersonic flow is a necessary but insufficient condition for shocks to develop. The development of a shock requires the characteristics of neighboring points in a flow to intersect. We demonstrate that the intersection of characteristics may be quantified via the knowledge of the Mach number. Finally, we examine three-dimensional simulations of hot Jovian atmospheres and demonstrate that shock formation is expected to occur mostly on the dayside hemisphere, upstream of the substellar point, because the enhanced temperatures near the substellar point provide a natural pressure barrier for the returning flow. Understanding the role of shocks in irradiated exoplanetary atmospheres is relevant to correctly modeling observables such as the peak offsets of infrared phase curves.

Two-Dimensional Crystallisation of Membrane Proteins and Structural Assessment

Two-dimensional (2D) crystallisation of Membrane proteins reconstitutes them into their native environment, the lipid bilayer. Electron crystallography allows the structural analysis of these regular protein–lipid arrays up to atomic resolution. The crystal quality depends on the protein purity, ist stability and on the crystallisation conditions. The basics of 2D crystallisation and different recent advances are reviewed and electron crystallography approaches summarised. Progress in 2D crystallisation, sample preparation, image detectors and automation of the data acquisition and processing pipeline makes 2D electron crystallography particularly attractive for the structural analysis of membrane proteins that are too small for single-particle analyses and too unstable to form three-dimensional (3D) crystals.