Fabrication, Electrical Characterization and Simulation of Thin Film Solar Cells: CdTe and CIGS Materials

CdTe and Cu(In,Ga)Se2 (CIGS) thin film solar cells are fabricated, electrically characterized and modelled in this thesis. We start from the fabrication of CdTe thin film devices where the R.F. magnetron sputtering system is used to deposit the CdS/CdTe based solar cells. The chlorine post-growth treatment is modified in order to uniformly cover the cell surface and reduce the probability of pinholes and shunting pathways creation which, in turn, reduces the series resistance. The deionized water etching is proposed, for the first time, as the simplest solution to optimize the effect of shunt resistance, stability and metal-semiconductor inter-diffusion at the back contact. In continue, oxygen incorporation is proposed while CdTe layer deposition. This technique has been rarely examined through R.F sputtering deposition of such devices. The above experiments are characterized electrically and optically by current-voltage characterization, scanning electron microscopy, x-ray diffraction and optical spectroscopy. Furthermore, for the first time, the degradation rate of CdTe devices over time is numerically simulated through AMPS and SCAPS simulators. It is proposed that the instability of electrical parameters is coupled with the material properties and external stresses (bias, temperature and illumination). Then, CIGS materials are simulated and characterized by several techniques such as surface photovoltage spectroscopy is used (as a novel idea) to extract the band gap of graded band gap CIGS layers, surface or bulk defect states. The surface roughness is scanned by atomic force microscopy on nanometre scale to obtain the surface topography of the film. The modified equivalent circuits are proposed and the band gap graded profiles are simulated by AMPS simulator and several graded profiles are examined in order to optimize their thickness, grading strength and electrical parameters. Furthermore, the transport mechanisms and Auger generation phenomenon are modelled in CIGS devices.

Valorisation of organic waste: new developments from proton nuclear magnetic resonance characterization

The last half-century has seen a continuing population and consumption growth, increasing the competition for land, water and energy. The solution can be found in the new sustainability theories, such as the industrial symbiosis and the zero waste objective. Reducing, reusing and recycling are challenges that the whole world have to consider. This is especially important for organic waste, whose reusing gives interesting results in terms of energy release. Before reusing, organic waste needs a deeper characterization. The non-destructive and non-invasive features of both Nuclear Magnetic Resonance (NMR) relaxometry and imaging (MRI) make them optimal candidates to reach such characterization. In this research, NMR techniques demonstrated to be innovative technologies, but an important work on the hardware and software of the NMR LAGIRN laboratory was initially done, creating new experimental procedures to analyse organic waste samples. The first results came from soil-organic matter interactions. Remediated soils properties were described in function of the organic carbon content, proving the importance of limiting the addition of further organic matter to not inhibit soil processes as nutrients transport. Moreover NMR relaxation times and the signal amplitude of a compost sample, over time, showed that the organic matter degradation of compost is a complex process that involves a number of degradation kinetics, as a function of the mix of waste. Local degradation processes were studied with enhanced quantitative relaxation technique that combines NMR and MRI. The development of this research has finally led to the study of waste before it becomes waste. Since a lot of food is lost when it is still edible, new NMR experiments studied the efficiency of conservation and valorisation processes: apple dehydration, meat preservation and bio-oils production. All these results proved the readiness of NMR for quality controls on a huge kind of organic residues and waste.

Secondary Microseisms Characterization and Green’s Function Extraction at the Larderello-Travale Geothermal Field (Italy)

Over the past ten years, the cross-correlation of long-time series of ambient seismic noise (ASN) has been widely adopted to extract the surface-wave part of the Green’s Functions (GF). This stochastic procedure relies on the assumption that ASN wave-field is diffuse and stationary. At frequencies <1Hz, the ASN is mainly composed by surface-waves, whose origin is attributed to the sea-wave climate. Consequently, marked directional properties may be observed, which call for accurate investigation about location and temporal evolution of the ASN-sources before attempting any GF retrieval. Within this general context, this thesis is aimed at a thorough investigation about feasibility and robustness of the noise-based methods toward the imaging of complex geological structures at the local (∼10-50km) scale. The study focused on the analysis of an extended (11 months) seismological data set collected at the Larderello-Travale geothermal field (Italy), an area for which the underground geological structures are well-constrained thanks to decades of geothermal exploration. Focusing on the secondary microseism band (SM;f>0.1Hz), I first investigate the spectral features and the kinematic properties of the noise wavefield using beamforming analysis, highlighting a marked variability with time and frequency. For the 0.1-0.3Hz frequency band and during Spring- Summer-time, the SMs waves propagate with high apparent velocities and from well-defined directions, likely associated with ocean-storms in the south- ern hemisphere. Conversely, at frequencies >0.3Hz the distribution of back- azimuths is more scattered, thus indicating that this frequency-band is the most appropriate for the application of stochastic techniques. For this latter frequency interval, I tested two correlation-based methods, acting in the time (NCF) and frequency (modified-SPAC) domains, respectively yielding esti- mates of the group- and phase-velocity dispersions. Velocity data provided by the two methods are markedly discordant; comparison with independent geological and geophysical constraints suggests that NCF results are more robust and reliable.

Development of new analytical procedures aimet at the characterization of artistic samples

In this thesis, new advances in the development of spectroscopic based methods for the characterization of heritage materials have been achieved. As concern FTIR spectroscopy new approaches aimed at exploiting near and far IR region for the characterization of inorganic or organic materials have been tested. Paint cross-section have been analysed by FTIR spectroscopy in the NIR range and an “ad hoc” chemometric approach has been developed for the elaboration of hyperspectral maps. Moreover, a new method for the characterization of calcite based on the use of grinding curves has been set up both in MIR and in FAR region. Indeed, calcite is a material widely applied in cultural heritage, and this spectroscopic approach is an efficient and rapid tool to distinguish between different calcite samples. Different enhanced vibrational techniques for the characterisation of dyed fibres have been tested. First a SEIRA (Surface Enhanced Infra-Red Absorption) protocol has been optimised allowing the analysis of colorant micro-extracts thanks to the enhancement produced by the addition of gold nanoparticles. These preliminary studies permitted to identify a new enhanced FTIR method, named ATR/RAIRS, which allowed to reach lower detection limits. Regarding Raman microscopy, the research followed two lines, which have in common the aim of avoiding the use of colloidal solutions. AgI based supports obtained after deposition on a gold-coated glass slides have been developed and tested spotting colorant solutions. A SERS spectrum can be obtained thanks to the photoreduction, which the laser may induce on the silver salt. Moreover, these supports can be used for the TLC separation of a mixture of colorants and the analyses by means of both Raman/SERS and ATR-RAIRS can be successfully reached. Finally, a photoreduction method for the “on fiber” analysis of colorant without the need of any extraction have been optimised.

Functionalization and characterization of magnetic nanoparticles for biomedical applications

Die vorliegende Arbeit beschäftigt sich mit der Oberflächenfunktionalisierung von MnO Nanopartikeln (NP). Durch die Verwendung und Verbesserung verschiedener Polymere durch die Einbindung von Poly (Ethylen Glycol) (PEG), gelang es, die Löslichkeit dieser Nanopartikel in wässrigen Lösungen sowie in Körperflüssigkeiten zu erhöhen. Zusätzlich konnten diese Nanopartikel deutlich besser steril filtriert werden und zeigten eine erhöhte Aktivität alsrnKontrastmittel im MRT. Vorläufige Ergebnisse für die Verwendung von Silika als Schutzhülle für MnO NP werden ebenfalls kurz erläutert. Die verwendeten Polymere besaßen dabei zugängliche Aminogruppen, die eine weitere Funktionalisierung durch Bio-aktiver Gruppen ermöglichte. Der Nachweis einer erfolgreichen Bindung durch verschiedene Methoden wie SDS-PAGE, Western- und Northern Blot sowie die Verwendung unterschiedlicher FluoreszenzMessungen wird ebenfalls diskutiert. MnO NP und anderer magnetischer NP werden weiterhin auf ihr toxisches Verhalten gegenüber Caki1 und HeLa Zellen getestet. Dabei zeigte sich, dass MnO NP, im Gegensatz zu einigen Kupferoxiden, quasi nicht toxisch waren und das Proliferationsverhalten dieser Zellen quasi nicht beeinflussten. Weiterhin wurde ein Fluoreszenzfarbstoff, konkret Protoporphyrin IX, an die Oberfläche von MnO NP angebracht.Diese konnten dann erfolgreich als Kontrastmittel in der MRT verwendet werden und zeigten vielversprechende Ergebnisse für die Photodynamische Therapie. Desweiteren wird die Synthese des Antikörpers gegen p53 ausführlich erläutert. Dabei wurde genau darauf geachtet,dass dieser Antikörper dann an MnO NP gebunden werden kann.

Surface-enhanced vibrational techniques for the detection of synthetic organic colorants

The research project object of this thesis is focused on the development of an advanced analytical system based on the combination of an improved thin layer chromatography (TLC) plate coupled with infrared (FTIR) and Raman microscopies for the detection of synthetic dyes. Indeed, the characterization of organic colorants, which are commonly present in mixtures with other components and in a very limited amount, still represents a challenging task in scientific analyses of cultural heritage materials. The approach provides selective spectral fingerprints for each compound, foreseeing the complementary information obtained by micro ATR-RAIRS-FTIR and SERS-Raman analyses, which can be performed on the same separated spot. In particular, silver iodide (AgI) applied on a gold coated slide is proposed as an efficient stationary phase for the discrimination of complex analyte mixtures, such as dyes present in samples of art-historical interest. The gold-AgI-TLC plate shows high performances related both to the chromatographic separation of analytes and to the spectroscopic detection of components. The use of a mid-IR transparent inorganic salt as the stationary phase avoids interferences of the background absorption in FTIR investigations. Moreover, by ATR microscopy measurements performed on the gold-AgI surface, a considerable enhancement in the intensity of spectra is observed. Complementary information can be obtained by Raman analyses, foreseeing a SERS activity of the AgI substrate. The method has been tested for the characterization of a mixture of three synthetic organic colorants widely used in dyeing processes: Brilliant Green (BG1), Rhodamine B (BV10) and Methylene Blue (BB9).

Characterization of integrated Bragg gratings in silicon-on-insulator

Silicon-on-insulator (SOI) is rapidly emerging as a very promising material platform for integrated photonics. As it combines the potential for optoelectronic integration with the low-cost and large volume manufacturing capabilities and they are already accumulate a huge amount of applications in areas like sensing, quantum optics, optical telecommunications and metrology. One of the main limitations of current technology is that waveguide propagation losses are still much higher than in standard glass-based platform because of many reasons such as bends, surface roughness and the very strong optical confinement provided by SOI. Such high loss prevents the fabrication of efficient optical resonators and complex devices severely limiting the current potential of the SOI platform. The project in the first part deals with the simple waveguides loss problem and trying to link that with the polarization problem and the loss based on Fabry-Perot Technique. The second part of the thesis deals with the Bragg Grating characterization from again the point of view of the polarization effect which leads to a better stop-band use filters. To a better comprehension a brief review on the basics of the SOI and the integrated Bragg grating ends up with the fabrication techniques and some of its applications will be presented in both parts, until the end of both the third and the fourth chapters to some results which hopefully make its precedent explanations easier to deal with.

Structural characterization of intermediate and metastable phases by electron microscopy

Structure characterization of nanocrystalline intermediates and metastable phases is of primary importance for a deep understanding of synthetic processes undergoing solid-to-solid state phase transitions. Understanding the evolution from the first nucleation stage to the final synthetic product supports not only the optimization of existing processes, but might assist in tailoring new synthetic paths. A systematic investigation of intermediates and metastable phases is hampered because it is impossible to produce large crystals and only in few cases a pure synthetic product can be obtained. Structure investigation by X-ray powder diffraction methods is still challenging on nanoscale, especially when the sample is polyphasic. Electron diffraction has the advantage to collect data from single nanoscopic crystals, but is limited by data incompleteness, dynamical effects and fast deterioration of the sample under the electron beam. Automated diffraction tomography (ADT), a recently developed technique, making possible to collect more complete three-dimensional electron diffraction data and to reduce at the same time dynamical scattering and beam damage, thus allowing to investigate even beam sensitive materials (f.e. hydrated phases and organics). At present, ADT is the only technique able to deliver complete three-dimensional structural information from single nanoscopic grains, independently from other surrounding phases. Thus, ADT is an ideal technique for the study of on-going processes where different phases exist at the same time and undergo several structural transitions. In this study ADT was used as the main technique for structural characterization for three different systems and combined subsequently with other techniques, among which high-resolution transmission electron microscopy (HRTEM), cryo-TEM imaging, X-ray powder diffraction (XRPD) and energy disperse X-ray spectroscopy (EDX).rnAs possible laser host materials, i.e. materials with a broad band emission in the near-infrared region, two unknown phases were investigated in the ternary oxide system M2O-Al2O3-WO3 (M = K, Na). Both phases exhibit low purity as well as non-homogeneous size distribution and particle morphology. The structures solved by ADT are also affected by pseudo-symmetry. rnSodium titanate nanotubes and nanowires are both intermediate products in the synthesis of TiO2 nanorods which are used as additives to colloidal TiO2 film for improving efficiency of dye-sensitized solar cells (DSSC). The structural transition from nantubes to nanowires was investigated in a step by step time-resolved study. Nanowires were discovered to consist of a hitherto unknown phase of sodium titanate. This new phase, typically affected by pervasive defects like mutual layer shift, was structurally determined ab-initio on the basis of ADT data. rnThe third system is related with calcium carbonate nucleation and early crystallization. The first part of this study is dedicated to the extensive investigations of calcium carbonate formation in a step by step analysis, up to the appearance of crystalline individua. The second part is dedicated to the structure determination by ADT of the first-to-form anhydrated phase of CaCO3: vaterite. An exhaustive structure analysis of vaterite had previously been hampered by diffuse scattering, extra periodicities and fast deterioration of the material under electron irradiation. rn

Development, characterization, and application of flowing atmospheric-pressure afterglow ionization for mass spectrometric analysis of ambient organic aerosols

Addressing current limitations of state-of-the-art instrumentation in aerosol research, the aim of this work was to explore and assess the applicability of a novel soft ionization technique, namely flowing atmospheric-pressure afterglow (FAPA), for the mass spectrometric analysis of airborne particulate organic matter. Among other soft ionization methods, the FAPA ionization technique was developed in the last decade during the advent of ambient desorption/ionization mass spectrometry (ADI–MS). Based on a helium glow discharge plasma at atmospheric-pressure, excited helium species and primary reagent ions are generated which exit the discharge region through a capillary electrode, forming the so-called afterglow region where desorption and ionization of the analytes occurs. Commonly, fragmentation of the analytes during ionization is reported to occur only to a minimum extent, predominantly resulting in the formation of quasimolecular ions, i.e. [M+H]+ and [M–H]– in the positive and the negative ion mode, respectively. Thus, identification and detection of signals and their corresponding compounds is facilitated in the acquired mass spectra. The focus of the first part of this study lies on the application, characterization and assessment of FAPA–MS in the offline mode, i.e. desorption and ionization of the analytes from surfaces. Experiments in both positive and negative ion mode revealed ionization patterns for a variety of compound classes comprising alkanes, alcohols, aldehydes, ketones, carboxylic acids, organic peroxides, and alkaloids. Besides the always emphasized detection of quasimolecular ions, a broad range of signals for adducts and losses was found. Additionally, the capabilities and limitations of the technique were studied in three proof-of-principle applications. In general, the method showed to be best suited for polar analytes with high volatilities and low molecular weights, ideally containing nitrogen- and/or oxygen functionalities. However, for compounds with low vapor pressures, containing long carbon chains and/or high molecular weights, desorption and ionization is in direct competition with oxidation of the analytes, leading to the formation of adducts and oxidation products which impede a clear signal assignment in the acquired mass spectra. Nonetheless, FAPA–MS showed to be capable of detecting and identifying common limonene oxidation products in secondary OA (SOA) particles on a filter sample and, thus, is considered a suitable method for offline analysis of OA particles. In the second as well as the subsequent parts, FAPA–MS was applied online, i.e. for real time analysis of OA particles suspended in air. Therefore, the acronym AeroFAPA–MS (i.e. Aerosol FAPA–MS) was chosen to refer to this method. After optimization and characterization, the method was used to measure a range of model compounds and to evaluate typical ionization patterns in the positive and the negative ion mode. In addition, results from laboratory studies as well as from a field campaign in Central Europe (F–BEACh 2014) are presented and discussed. During the F–BEACh campaign AeroFAPA–MS was used in combination with complementary MS techniques, giving a comprehensive characterization of the sampled OA particles. For example, several common SOA marker compounds were identified in real time by MSn experiments, indicating that photochemically aged SOA particles were present during the campaign period. Moreover, AeroFAPA–MS was capable of detecting highly oxidized sulfur-containing compounds in the particle phase, presenting the first real-time measurements of this compound class. Further comparisons with data from other aerosol and gas-phase measurements suggest that both particulate sulfate as well as highly oxidized peroxyradicals in the gas phase might play a role during formation of these species. Besides applying AeroFAPA–MS for the analysis of aerosol particles, desorption processes of particles in the afterglow region were investigated in order to gain a more detailed understanding of the method. While during the previous measurements aerosol particles were pre-evaporated prior to AeroFAPA–MS analysis, in this part no external heat source was applied. Particle size distribution measurements before and after the AeroFAPA source revealed that only an interfacial layer of OA particles is desorbed and, thus, chemically characterized. For particles with initial diameters of 112 nm, desorption radii of 2.5–36.6 nm were found at discharge currents of 15–55 mA from these measurements. In addition, the method was applied for the analysis of laboratory-generated core-shell particles in a proof-of-principle study. As expected, predominantly compounds residing in the shell of the particles were desorbed and ionized with increasing probing depths, suggesting that AeroFAPA–MS might represent a promising technique for depth profiling of OA particles in future studies.

Synthesis, characterization and formulation of a cathode active material: Copper nitroprusside

Nowadays, rechargeable Li-ion batteries play an important role in portable consumer devices. Formulation of such batteries is improvable by researching new cathodic materials that present higher performances of cyclability and negligible efficiency loss over cycles. Goal of this work was to investigate a new cathodic material, copper nitroprusside, which presents a porous 3D framework. Synthesis was carried out by a low-cost and scalable co-precipitation method. Subsequently, the product was characterized by means of different techniques, such as TGA, XRF, CHN elemental analysis, XRD, Mössbauer spectroscopy and cyclic voltammetry. Electrochemical tests were finally performed both in coin cells and by using in situ cells: on one hand, coin cells allowed different formulations to be easily tested, on the other operando cycling led a deeper insight to insertion process and both chemical and physical changes. Results of several tests highlighted a non-reversible electrochemical behavior of the material and a rapid capacity fading over time. Moreover, operando techniques report that amorphisation occurs during the discharge.

Structure determination of channel and transport proteins by high-resolution microscopy techniques

High-resolution microscopy techniques provide a plethora of information on biological structures from the cellular level down to the molecular level. In this review, we present the unique capabilities of transmission electron and atomic force microscopy to assess the structure, oligomeric state, function and dynamics of channel and transport proteins in their native environment, the lipid bilayer. Most importantly, membrane proteins can be visualized in the frozen-hydrated state and in buffer solution by cryo-transmission electron and atomic force microscopy, respectively. We also illustrate the potential of the scintillation proximity assay to study substrate binding of detergent-solubilized transporters prior to crystallization and structural characterization.

Hydrodynamic Characterization of Bile Salt Host-Guest Complexes by Pulsed Field Gradient NMR Spectroscopy

The work described herein is aimed at understanding primary and secondary aggregation of bile salt micelles and how micelles can perform chiral recognition of binapthyl analytes. Previous work with cholate and deoxycholate using micellar electrokinetic chromatography (MEKC) and nuclear magnetic resonance (NMR) has provided insightinto cholate and deoxycholate micelle formation, especially with respect to the critical micelle concentration (CMC). Chiral separations of the model analyte, 1,1â??-binaphthyl-2,2â??-diyl hydrogen phosphate (BNDHP), via cholate (C) and deoxycholate (DC) mediated MEKC separataions previously have shown the DC CMC to be 7-10 mM andthe cholate CMC at 14 mM at ph 12. A second model analyte,1,1â??-binaphthol (BN), was also previously investigated to probe micellar structure, but the MEKC data for this analyte implied a higher CMC, which may be interpreted as secondary aggregation. Thiswork extends the investigation of bile salts to include pulsed field gradient spin echo (PFGSE) NMR experiments being used to gain information about the size and degree of polydispersity of cholate and deoxycholate micelles. Concentrations of cholate below 10mM show a large variation in effective radius likely due to the existence of transient preliminary aggregates. The onset of the primary micelle shows a dramatic increase in effective radius of the micelle in cholate and deoxycholate. In the region of expectedsecondary aggregation a gradual increase of effective radius was observed with cholate; deoxycholate showed a persistent aggregate size in the secondary micelle region that is modulated by the presence of an analyte molecule. Effective radii of cholate anddeoxycholate (individually) were compared with and without R- and S-BNDHP in order to observe the effective radius difference of micelles with and without analyte present. The presence of S-BNDHP consistently resulted in a larger effective aggregate radius incholate and deoxycholate, confirming previous data of the S-BNDHP interacting more with the micelle than R-BNDHP. In total, various NMR techniques, like diffusion NMR can be used to gain a greater understanding of the bile salt micellization process and chiral resolution.

Characterization of harmonic tremor at Santiaguito volcano and its implications for eruption mechanisms

We observed Santiaguito volcano in southwestern Guatemala from March 2008 - March 2010. Seismic and infrasound data collected between January and March of 2009 contain records of many diverse processes occurring at the dacitic dome complex, including the recurrence of short lived (30-200 seconds in duration) harmonic tremor concurrent with ash poor gas emissions from the volcano. We employ several different analytical techniques to examine different portions of the tremor and source mechanisms. We use the parameters derived by this analysis to compare the feasibility of several suggested models of eruption mechanisms, and determine that this type of harmonic tremor is most justifiably generated by the flow of gas through crack networks generated by shear fracture along the magma conduit margin.

Fabrication and characterization of room temperature operating single electron transistors using focused ion beam technologies

The single electron transistor (SET) is a Coulomb blockade device, whose operation is based on the controlled manipulation of individual electrons. Single electron transistors show immense potential to be used in future ultra lowpower devices, high density memory and also in high precision electrometry. Most SET devices operate at cryogenic temperatures, because the charging energy is much smaller than the thermal oscillations. The room temperature operation of these devices is possible with sub- 10nm nano-islands due to the inverse dependance of charging energy on the radius of the conducting nano-island. The fabrication of sub-10nm features with existing lithographic techniques is a technological challenge. Here we present the results for the first room temperature operating SET device fabricated using Focused Ion Beam deposition technology. The SET device, incorporates an array of tungsten nano-islands with an average diameter of 8nm. The SET devices shows clear Coulomb blockade for different gate voltages at room temperature. The charging energy of the device was calculated to be 160.0 meV; the capacitance per junction was found to be 0.94 atto F; and the tunnel resistance per junction was calculated to be 1.26 G Ω. The tunnel resistance is five orders of magnitude larger than the quantum of resistance (26 k Ω) and allows for the localization of electrons on the tungsten nano-island. The lower capacitance of the device combined with the high tunnel resistance, allows for the Coulomb blockade effects observed at room temperature. Different device configurations, minimizing the total capacitance of the device have been explored. The effect of the geometry of the nano electrodes on the device characteristics has been presented. Simulated device characteristics, based on the soliton model have been discussed. The first application of SET device as a gas sensor has been demonstrated.

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.