Clinical and pharmacokinetic evaluation of S-ketamine for intravenous general anaesthesia in horses undergoing field castration.

BACKGROUND Intravenous anaesthetic drugs are the primary means for producing general anaesthesia in equine practice. The ideal drug for intravenous anaesthesia has high reliability and pharmacokinetic properties indicating short elimination and lack of accumulation when administered for prolonged periods. Induction of general anaesthesia with racemic ketamine preceded by profound sedation has already an established place in the equine field anaesthesia. Due to potential advantages over racemic ketamine, S-ketamine has been employed in horses to induce general anaesthesia, but its optimal dose remains under investigation. The objective of this study was to evaluate whether 2.5 mg/kg S-ketamine could be used as a single intravenous bolus to provide short-term surgical anaesthesia in colts undergoing surgical castration, and to report its pharmacokinetic profile. RESULTS After premedication with romifidine and L-methadone, the combination of S-ketamine and diazepam allowed reaching surgical anaesthesia in the 28 colts. Induction of anaesthesia as well as recovery were good to excellent in the majority (n = 22 and 24, respectively) of the colts. Seven horses required additional administration of S-ketamine to prolong the duration of surgical anaesthesia. Redosing did not compromise recovery quality. Plasma concentration of S-ketamine decreased rapidly after administration, following a two-compartmental model, leading to the hypothesis of a consistent unchanged elimination of the parent compound into the urine beside its conversion to S-norketamine. The observed plasma concentrations of S-ketamine at the time of first movement were various and did not support the definition of a clear cut-off value to predict the termination of the drug effect. CONCLUSIONS The administration of 2.5 mg/kg IV S-ketamine after adequate premedication provided good quality of induction and recovery and a duration of action similar to what has been reported for racemic ketamine at the dose of 2.2 mg/kg. Until further investigations will be provided, close monitoring to adapt drug delivery is mandatory, particularly once the first 10 minutes after injection are elapsed. Taking into account rapid elimination of S-ketamine, significant inter-individual variability and rapid loss of effect over a narrow range of concentrations a sudden return of consciousness has to be foreseen.

Cultural Legacies and Political Preferences: The Failure of Separatism in the Swiss Jura

The study of secession generally stresses the causal influence of cultural identities, political preferences, or ecological factors. Whereas these different views are often considered to be mutually exclusive, this paper proposes a two-stage model in which they are complementary. We posit that cultural identities matter for explaining secessionism, but not because of primordial attachments. Rather, religious and linguistic groups matter because their members are imbued with cultural legacies that lead to distinct political preferences – in this case preferences over welfare statism. Further, ecological constraints such as geography and topography affect social interaction with like-minded individuals. On the basis of both these political preferences and ecological constraints, individuals then make rational choices about the desirability of secession. Instrumental considerations are therefore crucial in explaining the decision to secede, but not in a conventional pocketbook manner. To examine this theory, we analyze the 2013 referendum on the secession of the Jura Bernois region from the Canton of Berne in Switzerland, using municipal level census and referendum data. The results lend support to the theory and suggest one way in which the politics of identity, based on factors like language and religion, can be fused with the politics of interest (preferences for more or less state intervention into the polity and economy) to better understand group behavior.

Lipid topogenesis - 35years on

Glycerophospholipids are the principal fabric of cellular membranes. The pathways by which these lipids are synthesized were elucidated mainly through the work of Kennedy and colleagues in the late 1950s and early 1960s. Subsequently, attention turned to cell biological aspects of lipids: Where in the cell are lipids synthesized? How are lipids integrated into membranes to form a bilayer? How are they sorted and transported from their site of synthesis to other cellular destinations? These topics, collectively termed 'lipid topogenesis', were the subject of a review article in 1981 by Bell, Ballas and Coleman. We now assess what has been learned about early events of lipid topogenesis, i.e. "lipid synthesis, the integration of lipids into membranes, and lipid translocation across membranes", in the 35years since the publication of this important review. We highlight the recent elucidation of the X-ray structures of key membrane enzymes of glycerophospholipid synthesis, progress on identifying lipid scramblase proteins needed to equilibrate lipids across membranes, and new complexities in the subcellular location and membrane topology of phosphatidylinositol synthesis revealed through a comparison of two unicellular model eukaryotes.

Insolation, erosion, and morphology of comet 67P/Churyumov-Gerasimenko

Context. The complex shape of comet 67P and its oblique rotation axis cause pronounced seasonal effects. Irradiation and hence activity vary strongly. Aims. We investigate the insolation of the cometary surface in order to predict the sublimation of water ice. The strongly varying erosion levels are correlated with the topography and morphology of the present cometary surface and its evolution. Methods. The insolation as a function of heliocentric distance and diurnal (spin dependent) variation is calculated using >10(5) facets of a detailed digital terrain model. Shading, but also illumination and thermal radiation by facets in the field of view of a specific facet are iteratively taken into account. We use a two-layer model of a thin porous dust cover above an icy surface to calculate the water sublimation, presuming steady state and a uniform surface. Our second model, which includes the history of warming and cooling due to thermal inertia, is restricted to a much simpler shape model but allows us to test various distributions of active areas. Results. Sublimation from a dirty ice surface yields maximum erosion. A thin dust cover of 50 pm yields similar rates at perihelion. Only about 6% of the surface needs to be active to match the observed water production rates at perihelion. A dust layer of 1 mm thickness suppresses the activity by a factor of 4 to 5. Erosion on the south side can reach more than 10 m per orbit at active spots. The energy input to the concave neck area (Hapi) during northern summer is enhanced by about 50% owing to self-illumination. Here surface temperatures reach maximum values along the foot of the Hathor wall. Integrated over the whole orbit this area receives the least energy input. Based on the detailed shape model, the simulations identify "hot spots" in depressions and larger pits in good correlation with observed dust activity. Three-quarters of the total sublimation is produced while the sub-solar latitude is south, resulting in a distinct dichotomy in activity and morphology. Conclusions. The northern areas display a much rougher morphology than what is seen on Imhotep, an area at the equator that will be fully illuminated when 67P is closer to the Sun. Self-illumination in concave regions enhance the energy input and hence erosion. This explains the early activity observed at Hapi. Cliffs are more prone to erosion than horizontal, often dust covered, areas, which leads to surface planation. Local activity can only persist if the forming cliff walls are eroding. Comet 67P has two lobes and also two distinct sides. Transport of material from the south to the north is probable. The morphology of the Imhotep plain should be typical for the terrains of the yet unseen southern hemisphere.

The differential staurosporine mediated G1 arrest in normal versus tumor cells is dependent on the retinoblastoma and Chk1 proteins

The ability to regulate cell cycle progression is one of the differences that separates normal from tumor cells. A protein, which is frequently mutated or deleted in a majority of tumor cells, is the retinoblastoma protein (pRb). Previously, we reported that normal cells, which have a wild-type Rb pathway, can be reversibly arrested in the G1 phase of the cell cycle by staurosporine (ST), while tumor cells were unaffected by this treatment. As a result, ST may be used to protect normal cells against the toxic affects of chemotherapy. Here we set out to determine the mechanism(s) by which ST can mediate a reversible G1 arrest in pRb positive cells. To this end, we used an isogenic cell model system of normal human mammary epithelial cells (HMEC) with either intact pRb+ (p53-) or p53+ (pRb-) treated with ST. Our results show that pRb+ cells treated with low concentrations of ST, arrested in the G1 phase of the cell cycle; however, in pRb - cells there was no response. This was verified as a true G 1 arrest in pRb+ cells by two different methods for monitoring cell cycle kinetics and in two additional model systems for Rb (i.e. pRb -/- mouse embryo fibroblasts, and downregulation of RB with siRNA). Our results indicated that ST-mediated G1 arrest required pRb, which in turn initiated a cascade of events leading to inhibition of CDK4 and CDK2 activities and up-regulation of p21 protein. Further assessment of this pathway revealed the novel finding that Chk1 expression and activity were required for the Rb-dependent, ST-mediated G1 arrest. In fact, overexpression of Chk1 facilitated recovery from ST-mediated G1 arrest, an effect only observed in RB+ cells. Collectively, our data suggest pRb is able to cooperate with Chk1 to mediate a G1 arrest in pRb+ cells, but not in pRb- cells. The elucidation of this pathway can help identify novel agents that can be used to protect cancer patients against the debilitating affects of chemotherapy, by targeting only the normal proliferating cells in the body that are otherwise destroyed. ^

Rb, Sr, Sm and Nd concentrations and isotopic composition of basalts from DSDP Leg 74 (Table 1)

Basement intersected in DSDP holes 525A, 528 and 527 on the Walvis Ridge consists of submarine basalt flows and pillows with minor intercalated sediments. These holes are situated on the crest and mid and lower northwest flank of a NNW-SSE-trending ridge block which would have closely paralleled the paleo mid-ocean ridge (Rabinowitz and LaBrecque, 1979 doi:10.1029/JB084iB11p05973, Moore et al. (1983 doi:10.1130/0016-7606(1983)94<907:TWRTDS>2.0.CO;2). The basalts were erupted approximately 70 m.y. ago, an age equivalent to that of immediately adjacent oceanic crust in the Angola Basin and coraistent with formation at the paleo mid-ocean ridge (Moore et al., 1983). The basalt types vary from aphyric quartz tholeiites on the ridge crest to highly plagioclase phyric olivine tholeiites on the ridge flank. These show systematic differences in incompatible trace element and isotopic composition. Many element and isotope ratio pairs form systematic trends with the ridge crest basalts at one end and the highly phyric ridge flank basalts at the other. The low 143Nd/144Nd (0.51238), 206Pb/204Pb (17.54), 207Pb/204Pb (15.47), 208Pb/204Pb (38.14) and high 87Sr/86Sr (0.70512) ratios of the ridge crest basalts suggest derivation from an old Nd/Sm-, Rb/Sr- and Pb/U-enriched mantle source. This isotopic signature is similar to that of alkaline basalts on Tristan da Cunha but offset to significantly lower Nd and Pb isotopic ratios. The isotopic ratio trends may be extrapolated beyond the ridge flank basalts with higher 143Nd/144Nd (0.51270), 206Pb/204Pb (18.32), 207Pb/204Pb (15.52), 208Pb/204Pb (38.77) and lower 87Sr/86Sr (0.70417) ratios in the direction of increasingly Nd/Sm-, Rb/Sr- and Pb/U-depleted source compositions. These isotopic correlations are equally consistent with mixing of depleted and enriched end member melts or partial melting of an inhomogeneous, variably enriched mantle source. However, observed Zr-Ba-Nb-Y interelement relationships are inconsistent with any simple two-component model of magma mixing, as might result from the rise of a lower mantle plume through the upper mantle. Incompatible element and Pb isotopic systematics also preclude extensive involvement of depleted (N-type) MORB material or its mantle sources. In our preferred petrogenetic model the Walvis Ridge basalts were derived by partial melting of mantle similar to an enriched (E-type) MORB source which had become heterogeneous on a small scale due to the introduction of small-volume melts and metasomatic fluids.

Concentrations of natural and anthropogenic organohalogens in water and air samples in the Canadian Arctic during IPY (2007-2008)

Shipboard measurements of organohalogen compounds in air and surface seawater were conducted in the Canadian Arctic in 2007-2008. Study areas included the Labrador Sea, Hudson Bay, and the southern Beaufort Sea. High volume air samples were collected at deck level (6 m), while low volume samples were taken at 1 and 15 m above the water or ice surface. Water samples were taken within 7 m. Water concentration ranges (pg/L) were as follows: alpha-hexachlorocyclohexane (alpha-HCH) 465-1013, gamma-HCH 150-254, hexachlorobenzene (HCB) 4.0-6.4, 2,4-dibromoanisole (DBA) 8.5-38, and 2,4,6-tribromoanisole (TBA) 4.7-163. Air concentration ranges (pg/m**3) were as follows: alpha-HCH 7.5-48, gamma-HCH 2.1-7.7, HCB 48-71, DBA 4.8-25, and TBA 6.4-39. Fugacity gradients predicted net deposition of HCB in all areas, while exchange directions varied for the other chemicals by season and locations. Net evasion of alpha-HCH from Hudson Bay and the Beaufort Sea during open water conditions was shown by air concentrations that averaged 14% higher at 1 m than 15 m. No significant difference between the two heights was found over ice cover. The alpha-HCH in air over the Beaufort Sea was racemic in winter (mean enantiomer fraction, EF = 0.504 ± 0.008) and nonracemic in late spring-early summer (mean EF = 0.476 ± 0.010). This decrease in EF was accompanied by a rise in air concentrations due to volatilization of nonracemic alpha-HCH from surface water (EF = 0.457 ± 0.019). Fluxes of chemicals during the southern Beaufort Sea open water season (i.e., Leg 9) were estimated using the Whitman two-film model, where volatilization fluxes are positive and deposition fluxes are negative. The means ± SD (and ranges) of net fluxes (ng/m**2/d) were as follows: alpha-HCH 6.8 ± 3.2 (2.7-13), gamma-HCH 0.76 ± 0.40 (0.26-1.4), HCB -9.6 ± 2.7 (-6.1 to -15), DBA 1.2 ± 0.69 (0.04-2.0), and TBA 0.46 ± 1.1 ng/m**2/d (-1.6 to 2.0).

Nd and Sr isotopic composition of Cretaceous MORB

Chemical and isotopic (Nd and Sr) compositions have been determined for 12 Cretaceous basaltic samples (108 Ma old) from Holes 417D and 418A of Legs 51,52 and 53. We have found that: (1) The chemical compositions are typical of MORB. They do not vary systematically with the stratigraphic positions of the analyzed samples; thus, the chemical evolution is independent of the eruption sequence that occurred at this Cretaceous ridge. (2) REE patterns for all rocks are characterized by a strong LREE depletion with (La/Sm)N = 0.38-0.50; no significant Eu anomalies are found; HREE are nearly flat or slightly depleted towards Yb-Lu and have 12-18 * chondritic abundances. Combining the results of previous studies, it suggests that no significant temporal and spatial variation in magma chemistry (especially for LIL elements) has occurred in the 'normal' ridge segments over the last 150 Ma. (3) lsotopically, 143Nd/144Nd ratios vary from 0.513026 to 0.513154, corresponding to epsilon-Nd(0) = +7.5 to +10, and they fall in the typical range of MORB. However, these rocks have unexpectedly high 87Sr/86Sr ratios (0.70355-0.70470) which are attributed to the result of seawater-rock interaction. (4) The Nd model ages (Tin), ranging from 1.53 to 2.47 (average 2.06) AE, suggest that the upper mantle source(s) underwent a large scale chemical differentiation leading to LREE and other LIL element depletion about 2 AE ago, assuming a simple two-stage model. More realistically, the variation in Tm(Nd) or epsilon-Nd could be derived from mixing of heterogeneous mantle sources that were a consequence of continuous mantle differentiation and continental formation. (5) Because of the low mg values (0.52-0.63), the analyzed basaltic rocks do not represent primary liquids of mantle melting. The variation in La/Sm ratios and TiO2 are not compatible with a model in which all rocks are genetically related by a simple fractional crystallization. Rather, it is proposed that the basaltic rocks might have been derived from some heterogeneous upper mantle source with or without later magmatic mixing, and followed by some shallow-level fraetionations.

Barium concentration in sediments of the central Equatorial Pacific

High resolution pore-water dissolved Ba concentration-depth profiles were determined at seven sites across an Equatorial Pacific productivity gradient from 12°S to 9°N, at 140°W. These data are important for understanding the physical, chemical, and biological controls on Ba recycling in the ocean, and for evaluating the paleo-oceanographic significance of Ba content in central Equatorial Pacific sediments. Pore-water Ba concentrations at all sites are higher than in the overlying bottom water, leading to a diffusive flux of Ba into the ocean. A pronounced subsurface concentration maximum exceeding barite solubility characterizes the dissolved Ba pore-water profiles, suggesting that the Ba regenerated in the upper few millimeters of sediment is not controlled by barite solubility. A few centimeters down-core Ba concentrations reach a relatively constant value of approximately barite saturation. The benthic Ba flux shows a clear zonal trend, with a maximum between 2°S and 2°N, most probably due to higher productivity at the equatorial divergence zone, and with lowest values at the southern and northern extremes of the transect. The dissolved Ba flux between 2°S and 2°N is ~30 nmol/cm**2 yr and drops to 6 nmol/cm**2 yr at 12°S. Even the lowest fluxes are significantly higher than those previously reported for the open ocean. In the Equatorial Pacific the calculated Ba recycling efficiency is about 70%. Thus, ~30% of the particulate Ba flux to the deep ocean is preserved in the sediments, compared with less than 1% for organic carbon and ~5% for biogenic silica. Mass balance calculation of the oceanic Ba cycle, using a two-box model, implies benthic Ba fluxes similar to those reported here for a steady-state ocean.

Noble gases sampled from three cruises in the eastern tropical Atlantic in September 2005 and June/July 2006

Upwelling velocities w in the equatorial band are too small to be directly observed. Here, we apply a recently proposed indirect method, using the observed helium isotope (3He or 4He) disequilibria in the mixed layer. The helium data were sampled from three cruises in the eastern tropical Atlantic in September 2005 and June/July 2006. A one-dimensional two-box model was applied, where the helium air-sea gas exchange is balanced by upwelling from 3He-rich water below the mixed layer and by vertical mixing. The mixing coefficients Kv were estimated from microstructure measurements, and on two of the cruises, Kv exceeded 1 x 10**-4 m**2/s, making the vertical mixing term of the same order of magnitude as the gas exchange and the upwelling term. In total, helium disequilibrium was observed on 54 stations. Of the calculated upwelling velocities, 48% were smaller than 1.0 x 10**-5 m/s, 19% were between 1.0 and 2.0 x 10**-5 m/s, 22% were between 2.0 and 4.0 x 10**-5 m/s, and on 11% of upwelling velocities exceeded this limit. The highest upwelling velocities were found in late June 2006. Meridional upwelling distribution indicated an equatorial asymmetry with higher vertical velocities between the equator and 1° to 2° south compared to north of the equator, particularly at 10°W. Associated heat flux into the mixed layer could be as high as 138 W/m**2, but this depends strongly on the chosen depths where the upwelled water comes from. By combining upwelling velocities with sea surface temperature and productivity distributions, a mean monthly equatorial upwelling rate of 19 Sv was estimated for June 2006 and a biweekly mean of 24 Sv was estimated for September 2005.

C and H geochemistry of altered oceanic crus of ODP/IODP Hole 1256D

Carbon and hydrogen concentrations and isotopic compositions were measured in 19 samples from altered oceanic crust cored in ODP/IODP Hole 1256D through lavas, dikes down to the gabbroic rocks. Bulk water content varies from 0.32 to 2.14 wt% with dD values from -64per mil to -25per mil. All samples are enriched in water relative to fresh basalts. The dD values are interpreted in terms of mixing between magmatic water and another source that can be either secondary hydrous minerals and/or H contained in organic compounds such as hydrocarbons. Total CO2, extracted by step-heating technique, ranges between 564 and 2823 ppm with d13C values from -14.9per mil to -26.6per mil. As for water, these altered samples are enriched in carbon relative to fresh basalts. The carbon isotope compositions are interpreted in terms of a mixing between two components: (1) a carbonate with d13C = -4.5per mil and (2) an organic compound with d13C = -26.6per mil. A mixing model calculation indicates that, for most samples (17 of 19), more than 75% of the total C occurs as organic compounds while carbonates represent less than 25%. This result is also supported by independent estimates of carbonate content from CO2 yield after H3PO4 attack. A comparison between the carbon concentration in our samples, seawater DIC (Dissolved Inorganic Carbon) and DOC (Dissolved Organic Carbon), and hydrothermal fluids suggests that CO2 degassed from magmatic reservoirs is the main source of organic C addition to the crust during the alteration process. A reduction step of dissolved CO2 is thus required, and can be either biologically mediated or not. Abiotic processes are necessary for the deeper part of the crust (>1000 mbsf) because alteration temperatures are greater than any hyperthermophilic living organism (i.e. T > 110 °C). Even if not required, we cannot rule out the contribution of microbial activity in the low-temperature alteration zones. We propose a two-step model for carbon cycling during crustal alteration: (1) when "fresh" oceanic crust forms at or close to ridge axis, alteration starts with hot hydrothermal fluids enriched in magmatic CO2, leading to the formation of organic compounds during Fischer-Tropsch-type reactions; (2) when the crust moves away from the ridge axis, these interactions with hot hydrothermal fluids decrease and are replaced by seawater interactions with carbonate precipitation in fractures. Taking into account this organic carbon, we estimate C isotope composition of mean altered oceanic crust at ? -4.7per mil, similar to the d13C of the C degassed from the mantle at ridge axis, and discuss the global carbon budget. The total flux of C stored in the altered oceanic crust, as carbonate and organic compound, is 2.9 ± 0.4 * 10**12 molC/yr.

Volcanic rocks from the submarine Vityaz Ridge: their age, chemical and isotopic compositions

Results of geological studies at the submarine Vityaz Ridge carried out during cruises 37 and 41 of R/V Akademik Lavrent'ev in 2005 and 2006 are reported. The studied area is located at an near-island trench of the slope in the central part of the Kuril Island arc. Morphologically it consists of two parts: an inner volcanic arc represented by the Great Kuril Range and an outer arc corresponding to the submarine Vityaz Ridge. Diverse rocks composing the basement and the sedimentary cover of the ridge were recovered by dredging. Based on K-Ar dating and geochemistry, volcanics were divided into Paleocene, Eocene, late Oligocene, and Pliocene-Pleistocene complexes. Each of the complexes reflects a tectonomagmatic stage in the ridge evolution. Geochemical and isotope data on the volcanics indicate contribution of ancient crustal material in the magma source and, correspondingly, formation of this structure on the continental basement. Two-stage model ages (TDM2) vary in a wide range from zero values in mafic rocks to 0.77 Ga in felsic varieties, pointing to presence of Precambrian protolith in the source of the felsic rocks of the Vityaz Ridge. The Pliocene-Pleistocene volcanics are classed with tholeiitic, calc-alkaline, and subalkaline series, which differ in alkali contents and REE fractionation. Values of (La/Sm)_n and (La/Yb)_n ratios vary from 0.74 and 0.84 in the tholeiitic varieties to 1.19 and 1.44 in the calc-alkaline and 2.32 and 3.73 in the subalkaline rocks. All three varieties occur within the same volcanic edifices and formed during differentiation of magmatic melts that were channeled along fault zones from the mantle source slightly enriched in crustal component.

Contraste de las formulaciones existentes para la estimación de la tasa de rebase en diques de abrigo portuario

La investigación que se resume en este artículo, resultante de diversos trabajos realizados en el CEDEX [10, 12], tiene por objeto contrastar las formulaciones desarrolladas por diferentes autores en materia de rebases en distintas tipologías de diques de abrigo. Para ello, se realizó en primer lugar la recopilación y el análisis de las formulaciones existentes para estimar la tasa de rebase sobre diques en talud y verticales. Posteriormente, se llevó a cabo el contraste de dichas formulaciones con los resultados obtenidos en dos ensayos realizados en el Centro de Estudios de Puertos y Costas del CEDEX, en Cabo Prioriño, Ferrol y en la alineación tercera de Melilla. The research that is summarized in this article, resultant of several studies carried out in the CEDEX [10, 12], is focused in the existing methods to estimate overtopping rates developed by different authors. For that, the summary was carry out firstly we collected and analyzed the existing formulae, particularly for rubble mound and vertical breakwaters. These formulations was compared with the results obtained in two hydraulic model tests of the Hydraulic Research Laboratory (the Center of Studies of Ports and Coasts of the CEDEX, Madrid, Spain). Keywords: overtopping rates, rubble mound breakwater, vertical breakwater, crown wall, run-up.

Diabatic quantum gates

At present, several models for quantum computation have been proposed. Adiabatic quantum computation scheme particularly offers this possibility and is based on a slow enough time evolution of the system, where no transitions take place. In this work, a new strategy for quantum computation is provided from the opposite point of view. The objective is to control the non-adiabatic transitions between some states in order to produce the desired exit states after the evolution. The model is introduced by means of an analogy between the adiabatic quantum computation and an inelastic atomic collision. By means of a simple two-state model, several quantum gates are reproduced, concluding the possibility of diabatic universal faulttolerant quantum computation. Going a step further, a new quantum diabatic computation model is glimpsed, where a carefully chosen Hamiltonian could carry out a non-adiabatic transition between the initial and the sought final state.

Laser ablation modelling of aluminium, silver and crystalline silicon for applications in photovoltaic technologies

Laser material processing is being extensively used in photovoltaic applications for both the fabrication of thin film modules and the enhancement of the crystalline silicon solar cells. The two temperature model for thermal diffusion was numerically solved in this paper. Laser pulses of 1064, 532 or 248 nm with duration of 35, 26 or 10 ns were considered as the thermal source leading to the material ablation. Considering high irradiance levels (108–109 W cm−2), a total absorption of the energy during the ablation process was assumed in the model. The materials analysed in the simulation were aluminium (Al) and silver (Ag), which are commonly used as metallic electrodes in photovoltaic devices. Moreover, thermal diffusion was also simulated for crystalline silicon (c-Si). A similar trend of temperature as a function of depth and time was found for both metals and c-Si regardless of the employed wavelength. For each material, the ablation depth dependence on laser pulse parameters was determined by means of an ablation criterion. Thus, after the laser pulse, the maximum depth for which the total energy stored in the material is equal to the vaporisation enthalpy was considered as the ablation depth. For all cases, the ablation depth increased with the laser pulse fluence and did not exhibit a clear correlation with the radiation wavelength. Finally, the experimental validation of the simulation results was carried out and the ability of the model with the initial hypothesis of total energy absorption to closely fit experimental results was confirmed.