Fabrication of optochemical and electrochemical sensors using thin films of porphyrin and phthalocyanine derivatives

This paper describes the fabrication of thin films of porphyrin and metallophthalocyanine derivatives on different substrates for the optochemical detection of HCl gas and electrochemical determination of L-cysteine (CySH). Solid state gas sensor for HCl gas was fabricated by coating meso-substituted porphyrin derivatives on glass slide and examined optochemical sensing of HCl gas. The concentration of gaseous HCl was monitored from the changes in the absorbance of Soret band. Among the different porphyrin derivatives, meso- tetramesitylporphyrin (MTMP) coated film showed excellent sensitivity towards HCl and achieved a detection limit of 0.03ppm HCl. Further, we have studied the self-assembly of 1,8,15,22-tetraaminometallophthalocyanine (4α-MTAPc; M = Co and Ni) from DMF on GC electrode. The CVs for the self-assembled monolayers (SAMs) of 4α-CoIITAPc and 4α-NiIITAPc show two pairs of well-defined redox couple corresponding to metal and ring. Using the 4α-CoIITAPc SAM modified electrode, sensitive and selective detection of L-cysteine was demonstrated. Further, the SAM modified electrode also successfully separates the oxidation potentials of AA and CySH with a peak separation of 320mV.

SERRS OF PHTHALOCYANINE DERIVATIVES LB MONOLAYER

SERRS of copper tetrakis(2,4-ditert-penthyl phenoxy) Phthalocyanine LB monolayer on Ag island film and RRS of the compound in solid state have been obtained. On comparison between the spectra, it is proposed that the macrocycle of the molecule on the Ag substrate is oriented parallel to the surface.

Adsorption kinetic investigations of phthalocyanine derivatives self assembled monolayers (SAMs) on gold: Temperature influence on the SAM formation process and quality

The ordered nano-structured surfaces, like self-assembled monolayers (SAMs) are of a great scientific interest, due to the low cost, simplicity, and versatility of this method. SAMs found numerous of applications in molecular electronics, biochemistry and optical devices. Phthalocyanine (Pc) complexes are of particular interest for the SAM preparation. These molecules exhibit fascinating physical properties and are chemically and thermally stable. Moreover their complex structure is advantageous for the fabrication of switchable surfaces. In this work the adsorption process of Pcs derivatives, namely, subphthalocyanines (SubPcB) and terbium (2TbPc) sandwich complexes on gold has been investigated. The influence of the molecular concentration, chain length of peripheral groups, and temperature on the film formation process has been examined using a number of techniques. The SAMs formation process has been followed in situ and in real time by means of second harmonic generation (SHG) and surface plasmon resonance (SPR) spectroscopy. To investigate the quality of the SAMs prepared at different temperatures atomic force microscopy (AFM) and X-Ray photoelectron spectroscopy (XPS)measurements were performed. Valuable information about SubPcB and 2TbPc adsorbtion process has been obtained in the frame of this work. The kinetic data, obtained with SHG and SPR, shows the best conformance with the first order Langmuir kinetic model. Comparing SHG and SPR results, it has been found, that the film formation occurs faster than the formation of chemical bonds. Such, the maximum amount of molecules on the surface is reached after 6 min for SubPcB and 30 min for 2TbPc. However, at this time the amount of formed chemicals bonds is only 10% and 40% for SubPcB and 2TbPc, respectively. The most intriguing result, among others, was obtained at T = 2 °C, where the formation of the less dense SAMs have been detected with SHG.However, analyzing XPS and AFM data, it has been revealed, that there is the same amount of molecules on the surface at both temperature T = 2 °C, and T = 21 °C, but the amount of formed chemicals bond is different. At T = 2 °C molecules form aggregates, therefore many of available anchor groups stay unattached.

Síntese de derivados oligoméricos de porfirinas e ftalocianinas

A presente dissertação contempla estudos de funcionalização de 5,10,15,20- tetrafenilporfirina via grupos nitro e amino e a preparação de sistemas porfirina-ftalocianina. Este trabalho encontra-se dividido em quatro partes. Na primeira parte descrevem-se as características gerais de porfirinas e ftalocianinas bem como algumas metodologias de síntese utilizadas na sua preparação e suas potenciais aplicações. Na segunda parte desta dissertação descreve-se a funcionalização de 5,10,15,20-tetrafenilporfirina com arilaminas recorrendo a duas rotas sintéticas diferentes. A reacção de 2-nitro-5,10,15,20-tetrafenilporfirina com anilina ou aminas aromáticas substituídas com grupos dadores de electrões permitiu, através do ataque do nucleófilo ao carbono beta-pirrólico onde está ligado o grupo nitro, ataque ipso, a obtenção de derivados do tipo 2-arilaminoporfirinas e derivados porfirínicos de anéis fundidos, sendo estes últimos resultantes da ciclização oxidativa de 2-arilaminoporfirinas. A reacção entre (2-amino- 5,10,15,20-tetrafenilporfirinato)níquel(II) e brometos de arilo na presença de paládio, reacção de aminação de Buchwald-Hartwig, permitiu, após descomplexação, a preparação de novos derivados do tipo 2- arilaminoporfirinas com grupos substituintes dadores e aceitadores de electrões. Um dos derivados porfirínicos de anéis fundidos foi submetido a reacção de complexação com diferentes iões metálicos e foram estudadas as respectivas propriedades fotoquímicas e electroquímicas. Esses estudos revelaram que estes compostos são bons geradores de oxigénio singuleto e que sofrem processos de oxidação-redução electroquimicamente reversíveis. Esta metodologia foi estendida ainda a brometos de hetarilo (derivados de piridina e tiofeno). Recorrendo ao acoplamento, em condições de Buchwald-Hartwig, de complexos metálicos da 5,10,15,20-tetrafenilporfirina, funcionalizados com grupos amino e bromo, preparam-se dímeros porfirina-amino-porfirina, cujos espectros electrónicos revelam a existência de uma boa “comunicação electrónica” entre as duas subunidades. A terceira parte descreve a síntese de sistemas porfirina-ftalocianina. Recorrendo à condensação estatística entre a 5,10,15,20-tetrafenilporfirina substituída com um grupo ftalonitrilo na posição beta-pirrólica com ftalonitrilo ou ftalonitrilo substituído foram obtidas díades porfirina-ftalocianina onde as duas subunidades se encontram directamente ligadas ou fundidas. Os porfirinilftalonitrilos necessários para a síntese das diferentes díades foram preparados através da reacção de adição do fumaronitrilo à 5,10,15,20- tetrafenilporfirina funcionalizada com o grupo 1,3-butadienilo ou vinilo, seguida de oxidação do aducto resultante. O acoplamento catalisado por paládio entre (2-bromo-5,10,15,20- tetrafenilporfirinato)zinco(II) e [9(10),16(17),23(24)-tri-terc-butil-2- etinilftalocianinato]zinco(II) permitiu a síntese de uma díade porfirinaftalocianina com as duas unidades ligadas por um grupo etinilo. Uma pentíade porfirina-ftalocianina foi obtida através da ciclotetramerização de um dos porfirinilftalonitrilos. A comparação dos espectros electrónicos das diferentes classes de sistemas revela que as correspondentes propriedades electrónicas são altamente afectadas pela distância entre as subunidades e também pelo número de cromóforos presentes no sistema. Os estudos fotofísicos de alguns dos novos compostos acima referidos permitiram verificar a ocorrência eficiente de transferência de energia da subunidade porfirínica para a da ftalocianina, capacidade essa que permitirá a estes sistemas serem usados para modelar o processo fotossintético. Na última parte descrevem-se, pormenorizadamente, todas as experiências efectuadas e as caracterizações espectroscópicas, nomeadamente de espectroscopia de ressonância magnética nuclear (RMN), espectrometria de massa e espectrofotometria de UV-vis, dos compostos sintetizados. Nalguns casos recorreu-se ainda a técnicas de RMN bidimensionais como COSY, HSQC, HMBC, NOESY e ROESY.

Selbstassemblierte Monolagen und magnetisch strukturierte Submonolagen mit dia- und paramagnetischen Molekülen sowie Einzelmolekülmagneten auf Phthalocyanin- und Subphthalocyanin-Basis

Mit dieser Arbeit wurde die Selbstassemblierung von dia- und paramagnetischen Molekülen sowie Einzelmolekülmagneten auf Goldsubstraten und magnetisch strukturierten Substraten untersucht. Dazu wurden drei verschiedene Klassen an Phthalocyaninderivaten verwendet: Diamagnetische Subphthalocyanine, paramagnetische Phthalocyaninatometalle und Diphthalocyaninatolanthanidkomplexe. Alle synthetisierten Verbindungen sind peripher thioethersubstituiert. Die Alkylketten (a: n-C8H17, b: n-C12H25) vermitteln die Löslichkeit in vielen organischen Solventien und sorgen für eine geordnete Assemblierung auf einer Oberfläche, wobei die Bindung auf Gold hauptsächlich über die Schwefelatome stattfindet. Die aus Lösung abgeschiedenen selbstassemblierten Monolagen wurden mit XPS, NEXAFS-Spektroskopie und ToF-SIMS untersucht. Bei der Selbstassemblierung auf magnetisch strukturierten Substraten stehen die Moleküle unter dem Einfluss magnetischer Streufelder und binden bevorzugt nur in bestimmten Bereichen. Die gebildeten Submonolagen wurden zusätzlich mit X-PEEM untersucht. Die erstmals dargestellten Manganphthalocyanine [MnClPc(SR)8] 1 wurden ausgehend von MnCl2 erhalten. Hier fand bei der Aufarbeitung an Luft eine Oxidation zu Mangan(III) statt; +III ist die stabilste Oxidationsstufe von Mangan in Phthalocyaninen. Der Nachweis des axialen Chloridoliganden erfolgte mit Massenspektrometrie und FIR- sowie Raman-Spektroskopie. SQUID-Messungen haben gezeigt, dass die Komplexe 1 vier ungepaarte Elektronen haben. Bei den Subphthalocyaninen [BClSubpc(SR)6] 2 wurde der axiale Chloridoligand mit dem stäbchenförmigen Phenolderivat 29-H substituiert und die erfolgreiche Ligandensubstitution durch NMR- und IR-Spektroskopie sowie Massenspektrometrie an den Produkten [BSubpc(SR)6(29)] 30 belegt. Der Radikalcharakter der synthetisierten Terbiumkomplexe [Tb{Pc(SR)8}2] 3 wurde spektroskopisch nachgewiesen; SQUID-Messungen ergaben, dass es sich um Einzelmolekülmagnete mit einer Energiebarriere U des Doppelpotentialtopfs von 880 K oder 610 cm-1 bei 3a handelt. Zunächst wurden die SAMs der Komplexverbindungen 1, 2, 30 und 3 auf nicht magnetisch strukturierten Goldsubstraten untersucht. Die Manganphthalocyanine 1 bilden geordnete SAMs mit größtenteils flach liegenden Molekülen, wie die XPS-, NEXAFS- und ToF-SIMS-Analyse zeigte. Die Mehrzahl der Thioether-Einheiten ist auf Gold koordiniert und die Alkylketten zeigen ungeordnet von der Oberfläche weg. Bei der Adsorption findet eine Reduktion zu Mangan(II) statt und der axiale Chloridoligand wird abgespalten. Das beruht auf dem sog. Oberflächen-trans-Effekt. Im vorliegenden Fall übt die Metalloberfläche einen stärkeren trans-Effekt als der axiale Ligand aus, was bisher experimentell noch nicht beobachtet wurde. Die thioethersubstituierten Subphthalocyanine 2 und 30 sowie die Diphthalocyaninatoterbium-Komplexe 3 sind ebenfalls für SAMs geeignet. Ihre Monolagen wurden mit XPS und NEXAFS-Spektroskopie untersucht, und trotz einer gewissen Unordnung in den Filmen liegen die Moleküle jeweils im Wesentlichen flach auf der Goldoberfläche. Vermutlich sind bei diesen Systemen auch die Alkylketten größtenteils parallel zur Oberfläche orientiert. Im Gegensatz zu den Manganphthalocyaninen 1 tritt bei 2b, 30a, 30b und 3b neben der koordinativen Bindung der Schwefelatome auf Gold auch eine für Thioether nicht erwartete kovalente Au–S-Bindung auf, die durch C–S-Bindungsbruch unter Abspaltung der Alkylketten ermöglicht wird. Der Anteil, zu dem dieser Prozess stattfindet, scheint nicht mit der Molekülstruktur zu korrelieren. Selbstassemblierte Submonolagen auf magnetisch strukturierten Substraten wurden mit dem diamagnetischen Subphthalocyanin 2b hergestellt. Der Nachweis der Submonolagen war schwierig und gelang schließlich durch eine Kombination von ToF-SIMS, NEXAFS Imaging und X-PEEM. Die Analyse der ToF-SIMS-Daten zeigte, dass tatsächlich eine Modulation der Verteilung der Moleküle auf einem unterwärts magnetisch strukturierten Substrat eintritt. Mit X-PEEM konnte die magnetische Struktur der ferromagnetischen Schicht des Substrats direkt der Verteilung der adsorbierten Moleküle zugeordnet werden. Die Subphthalocyanine 2b adsorbieren nicht an den Domänengrenzen, sondern vermehrt dazwischen. Auf Substraten mit abwechselnd 6.5 und 3.5 µm breiten magnetischen Domänen binden die Moleküle bevorzugt in den Bereichen geringster magnetischer Streufeldgradienten, also den größeren Domänen. Solche Substrate wurden für die ToF-SIMS- und X-PEEM-Messungen verwendet. Bei größeren magnetischen Strukturen mit ca. 400 µm breiten Domänen, wie sie aufgrund der geringeren Ortsauflösung dieser Methode für NEXAFS Imaging eingesetzt wurden, binden die Moleküle dann in allen Domänen. Die diamagnetischen Moleküle werden nach dieser Interpretation aus dem inhomogenen Magnetfeld über der Probenoberfläche heraus gedrängt und verhalten sich analog makroskopischer Diamagnete. Die eindeutige Detektion der Moleküle auf den magnetisch strukturierten Substraten konnte bisher nur für die diamagnetischen Subphthalocyanine 2b erfolgen. Um die Interpretation ihres Verhaltens bei der Selbstassemblierung in einem inhomogenen Magnetfeld weiter voranzutreiben, wurde das Subphthalocyanin 37b dargestellt, welches ein stabiles organisches TEMPO-Radikal in seinem axialen Liganden enthält. Das paramagnetische Subphthalocyanin 37b sollte auf den magnetisch strukturierten Substraten in Regionen starker magnetischer Streufelder binden und damit das entgegengesetzte Verhalten zu den diamagnetischen Subphthalocyaninen 2b zeigen. Aus Zeitgründen konnte dieser Nachweis im Rahmen dieser Arbeit noch nicht erbracht werden.

Electroreduction of oxygen on some novel cobalt phthalocyanine complexes

The catalytic activity of cobalt phthalocyanine monomer and some of its polymeric derivatives towards the electroreduction of molecular oxygen in salt and alkaline solutions is examined. It is found that most of these complexes exhibit a higher catalytic activity than the cobalt phthalocyanine monomer.

Phthalocyanine-based systems: synthesis, properties and applications

A presente dissertação descreve a síntese e o estudo de propriedades fotofísicas e biológicas de sistemas baseados em ftalocianinas. Os sistemas sintetizados consistem em ftalocianinas contendo unidades de D-galactose, ciclodextrinas, [60]fulereno ou porfirinas. Na primeira parte da dissertação aborda-se a síntese de ftalocianinas substituídas com unidades de D-galactose, com o intuito de serem utilizadas como fotossensibilizadores em terapia fotodinâmica. A eficiência de geração de oxigénio singuleto induzida pelas glicoftalocianinas, bem como a sua actividade fotodinâmica in vitro em células HeLa são avaliadas. Adicionalmente são descritos os estudos de MALDI-MS/MS de dois pares de glicoftalocianinas isoméricas. Numa segunda parte, descrevem-se a síntese e a caracterização de sistemas supramoleculares de ftalocianinas ligadas a ciclodextrinas ou [60]fulereno. Foram desenvolvidas metodologias sintéticas para obter complexos ftalocianínicos de ruténio(II) substituídos axialmente com derivados de ciclodextrinas, e complexos ftalocianínicos de zinco(II) ligados covalentemente a ciclodextrinas. Descreve-se também a síntese de uma díade ftalocianina- [60]fulereno. A estratégia sintética envolve um peculiar primeiro passo, uma reacção de aminação redutiva, no qual se obtém uma ftalocianina funcionalizada com glicina, seguida da habitual reacção de cicloadição 1,3- dipolar. Na terceira parte, abordam-se sistemas dador-aceitador de díades porfirinaftalocianina, como modelos promissores para mimetizar sistemas fotossintéticos naturais. Diferentes tipos de díades porfirina-ftalocianina, ligadas através do grupo fenilo da posição meso ou da posição β-pirrólica da porfirina, são sintetizadas e caracterizadas. As abordagens sintéticas envolvem reacções de aminação catalisadas por paládio e reacções de condensação cruzada de dois ftalonitrilos adequadamente substituídos. Também se descrevem processos de transferência de energia e electrónica fotoinduzidos que se formam a partir das díades porfirina-ftalocianina.

Vibrational spectroscopy of phthalocyanine and naphthalocyanine in sandwich-type (na)phthalocyaninato and porphyrinato rare earth complexes : (Part 10) The infrared and Raman characteristics of phthalocyanine in heteroleptic bis(phthalocyaninato) rare earth complexes with decreased molecular symmetry

The infrared (IR) spectroscopic data for a series of eleven heteroleptic bis(phthalocyaninato) rare earth complexes MIII(Pc)[Pc(α-OC5H11)4] (M = Sm–Lu, Y) [H2Pc = unsubstituted phthalocyanine, H2Pc(α-OC5H11)4 = 1,8,15,22-tetrakis(3-pentyloxy)phthalocyanine] have been collected with 2 cm−1 resolution. Raman spectroscopic properties in the range of 500–1800 cm−1 for these double-decker molecules have also been comparatively studied using laser excitation sources emitting at 632.8 and 785 nm. Both the IR and Raman spectra for M(Pc)[Pc(α-OC5H11)4] are more complicated than those of homoleptic bis(phthalocyaninato) rare earth analogues due to the decreased molecular symmetry of these double-decker compounds, namely C4. For this series, the IR Pc√− marker band appears as an intense absorption at 1309–1317 cm−1, attributed to the pyrrole stretching. With laser excitation at 632.8 nm, Raman vibrations derived from isoindole ring and aza stretchings in the range of 1300–1600 cm−1 are selectively intensified. In contrast, when excited with laser radiation of 785 nm, the ring radial vibrations of isoindole moieties and dihedral plane deformations between 500 and 1000 cm−1 for M(Pc)[Pc(α-OC5H11)4] intensify to become the strongest scatterings. Both techniques reveal that the frequencies of pyrrole stretching, isoindole breathing, isoindole stretchings, aza stretchings and coupling of pyrrole and aza stretchings depend on the rare earth ionic size, shifting to higher energy along with the lanthanide contraction due to the increased ring-ring interaction across the series. The assignments of the vibrational bands for these compounds have been made and discussed in relation to other unsubstituted and substituted bis(phthalocyaninato) rare earth analogues, such as M(Pc)2 and M(OOPc)2 [H2OOPc = 2,3,9,10,16,17,23,24-octakis(octyloxy)phthalocyanine].

Vibrational spectroscopy of phthalocyanine and naphthalocyanine in sandwich-type (na)phthalocyaninato and porphyrinato rare earth complexes: Part 14. The infrared and Raman characteristics of phthalocyanine in “pinwheel-like” homoleptic bis[1,8,15,22-tetrakis(3-pentyloxy)phthalocyaninato] rare earth(III) double-decker complexes

The infrared (IR) spectroscopic data and Raman spectroscopic properties for a series of 13 “pinwheel-like” homoleptic bis(phthalocyaninato) rare earth complexes M[Pc(α-OC5H11)4]2 [M = Y and Pr–Lu except Pm; H2Pc(α-OC5H11)4 = 1,8,15,22-tetrakis(3-pentyloxy)phthalocyanine] have been collected and comparatively studied. Both the IR and Raman spectra for M[Pc(α-OC5H11)4]2 are more complicated than those of homoleptic bis(phthalocyaninato) rare earth analogues, namely M(Pc)2 and M[Pc(OC8H17)8]2, but resemble (for IR) or are a bit more complicated (for Raman) than those of heteroleptic counterparts M(Pc)[Pc(α-OC5H11)4], revealing the decreased molecular symmetry of these double-decker compounds, namely S8. Except for the obvious splitting of the isoindole breathing band at 1110–1123 cm−1, the IR spectra of M[Pc(α-OC5H11)4]2 are quite similar to those of corresponding M(Pc)[Pc(α-OC5H11)4] and therefore are similarly assigned. With laser excitation at 633 nm, Raman bands derived from isoindole ring and aza stretchings in the range of 1300–1600 cm−1 are selectively intensified. The IR spectra reveal that the frequencies of pyrrole stretching and pyrrole stretching coupled with the symmetrical CH bending of –CH3 groups are sensitive to the rare earth ionic size, while the Raman technique shows that the bands due to the isoindole stretchings and the coupled pyrrole and aza stretchings are similarly affected. Nevertheless, the phthalocyanine monoanion radical Pc′− IR marker band of bis(phthalocyaninato) complexes involving the same rare earth ion is found to shift to lower energy in the order M(Pc)2 > M(Pc)[Pc(α-OC5H11)4] > M[Pc(α-OC5H11)4]2, revealing the weakened π–π interaction between the two phthalocyanine rings in the same order.

Raman spectroscopic characteristics of phthalocyanine and naphthalocyanine in sandwich-type phthalocyaninato and porphyrinato rare earth complexes Part 5. – Raman spectroscopic characteristics of naphthalocyanine in mixed [tetrakis(4-tert-butylphenyl)porphyrinato]-(naphthalocyaninato) rare earth double-deckers

Raman spectra were recorded in the range 400–1800 cm−1 for a series of 15 mixed \[tetrakis(4-tert-butylphenyl)porphyrinato](2,3-naphthalocyaninato) rare earth double-deckers M(TBPP)(Nc) (M = Y; La–Lu except Pm) using laser excitation at 632.8 and 785 nm. Comparisons with bis(naphthalocyaninato) rare earth counterparts reveal that the vibrations of the metallonaphthalocyanine M(Nc) fragment dominate the Raman features of M(TBPP)(Nc). When excited with radiation of 632.8 nm, the most intense vibration appears at about 1595 cm−1, due to the naphthalene stretching. These complexes exhibit the marker Raman band for Nc•− as a medium-intense band in the range 1496–1507 cm−1, attributed to the coupling of pyrrole and aza stretching, while the marker Raman band of Nc2− in intermediate-valence Ce(TBPP)(Nc) appears as a strong band at 1493 cm−1 and is due to the isoindole stretchings. By contrast, when excited with radiation of 785 nm that is in close resonance with the main Q absorption band of the naphthalocyanine ligand, the ring radial vibrations at ca 680 and 735 cm−1 for MIII(TBPP)(Nc) are selectively intensified and are the most intense bands. For the cerium double-decker, the most intense vibration also acting as the marker Raman band of Nc2− appears at 1497 cm−1 with contributions from both pyrrole CC and aza CN stretches. The same vibrational modes show weak to medium intensity scattering at 1506–1509 cm−1 for MIII(TBPP)(Nc) and this is the marker Raman band of Nc•− when thus excited. The scatterings due to the Nc breathings, ring radial vibration, aza group stretchings, naphthalene stretchings, benzoisoindole stretchings and the coupling of pyrrole CC and aza CN stretchings in MIII(TBPP)(Nc) are all slightly blue shifted along with the decrease in rare earth ionic radius, confirming the effects of increased ring–ring interactions on the Raman characteristics of naphthalocyanine in the mixed ring double-deckers.

Numerical simulation for the 3D seepage flow with fractional derivatives in porous media

In this paper, the numerical simulation of the 3D seepage flow with fractional derivatives in porous media is considered under two special cases: non-continued seepage flow in uniform media (NCSFUM) and continued seepage flow in non-uniform media (CSF-NUM). A fractional alternating direction implicit scheme (FADIS) for the NCSF-UM and a modified Douglas scheme (MDS) for the CSF-NUM are proposed. The stability, consistency and convergence of both FADIS and MDS in a bounded domain are discussed. A method for improving the speed of convergence by Richardson extrapolation for the MDS is also presented. Finally, numerical results are presented to support our theoretical analysis.

Modified alternating direction methods for solving a two-dimensional non-continuous seepage flow with fractional derivatives

In this paper, a two-dimensional non-continuous seepage flow with fractional derivatives (2D-NCSF-FD) in uniform media is considered, which has modified the well known Darcy law. Using the relationship between Riemann-Liouville and Grunwald-Letnikov fractional derivatives, two modified alternating direction methods: a modified alternating direction implicit Euler method and a modified Peaceman-Rachford method, are proposed for solving the 2D-NCSF-FD in uniform media. The stability and consistency, thus convergence of the two methods in a bounded domain are discussed. Finally, numerical results are given.

Numerical methods for fractional partial differential equations with Riesz space fractional derivatives

In this paper, we consider the numerical solution of a fractional partial differential equation with Riesz space fractional derivatives (FPDE-RSFD) on a finite domain. Two types of FPDE-RSFD are considered: the Riesz fractional diffusion equation (RFDE) and the Riesz fractional advection–dispersion equation (RFADE). The RFDE is obtained from the standard diffusion equation by replacing the second-order space derivative with the Riesz fractional derivative of order αset membership, variant(1,2]. The RFADE is obtained from the standard advection–dispersion equation by replacing the first-order and second-order space derivatives with the Riesz fractional derivatives of order βset membership, variant(0,1) and of order αset membership, variant(1,2], respectively. Firstly, analytic solutions of both the RFDE and RFADE are derived. Secondly, three numerical methods are provided to deal with the Riesz space fractional derivatives, namely, the L1/L2-approximation method, the standard/shifted Grünwald method, and the matrix transform method (MTM). Thirdly, the RFDE and RFADE are transformed into a system of ordinary differential equations, which is then solved by the method of lines. Finally, numerical results are given, which demonstrate the effectiveness and convergence of the three numerical methods.